為空照圖進行空間對位

文章推薦指數: 80 %
投票人數:10人

你還會學到這些¶ · 下載超高解析度的公領域影像資料 · 在QGIS 中使用OpenLayers 附加元件 · 使用cs2cs 指令工具在不同的投影系統之間轉換坐標 · 使用手邊已經過空間對位的圖層 ... QGISTutorialsandTips v1 TutorialsList 簡介 製作地圖(QGIS3) WorkingwithAttributes(QGIS3) 匯入工作表或CSV檔(QGIS3) BasicVectorStyling(QGIS3) CalculatingLineLengthsandStatistics(QGIS3) BasicRasterStylingandAnalysis(QGIS3) RasterMosaicingandClipping(QGIS3) WorkingwithTerrainData(QGIS3) WorkingwithWMSData(QGIS3) WorkingwithProjections(QGIS3) 為紙本地圖進行空間對位(QGIS3) 為空照圖進行空間對位(QGIS3) DigitizingMapData(QGIS3) 資料表連接(QGIS3) 空間連接(QGIS3) PerformingSpatialQueries(QGIS3) NearestNeighborAnalysis(QGIS3) SamplingRasterDatausingPointsorPolygons(QGIS3) CalculatingRasterArea(QGIS3) 製作熱區圖(QGIS3) AnimatingTimeSeriesData(QGIS3) HandlingInvalidGeometries(QGIS3) AdvancedRasterAnalysis(QGIS3) InterpolatingPointData(QGIS3) CalculatingStreetIntersectionDensity(QGIS3) 使用處理框架進行批次處理(QGIS3) 利用處理建模自動執行複雜的工作流程(QGIS3) AutomatingMapCreationwithPrintLayoutAtlas(QGIS3) MultiCriteriaOverlayAnalysis(QGIS3) BasicNetworkVisualizationandRouting(QGIS3) LocatingNearestFacilitywithOrigin-DestinationMatrix(QGIS3) ServiceAreaAnalysisusingOpenrouteservice(QGIS3) TravelTimeAnalysiswithUberMovement(QGIS3) CalculatingArealMeanRainfall(QGIS3) Python程式設計的初步上手(QGIS3) RunningProcessingAlgorithmsviaPython(QGIS3) 使用自訂的Python表達式函數(QGIS3) 撰寫Python腳本並在處理框架下執行(QGIS3) 以Python製作附加元件(QGIS3) BuildingaProcessingPlugin(QGIS3) 執行QGIS工作排程 WebMappingwithQGIS2Web(QGIS3) 用QTiles製作底圖 UsingPlugins(QGIS3) 搜尋和下載開放街圖資料 QGIS學習資源 資料來源 使用處理框架進行批次處理(QGIS2) Author UjavalGandhi Follow@spatialthoughts WanttolearnQGISinastructuredway?CheckoutSpatialThoughtsformorelearningmaterialsandinstructor-ledonlineprogramswithQGIS.orgcertification. ChangeLanguage English Chinese(Taiwan) Dutch Finnish French German Greek Indonesian Italian Korean Persian Portuguese Romanian Russian Slovenian Spanish Thai Turkish Ukrainian Vietnamese 為空照圖進行空間對位¶ 警告 本教學已經過時,如需閱讀新的版本請前往為空照圖進行空間對位(QGIS3)。

在為紙本地圖進行空間對位教學中,我們介紹了一種QGIS中進行空間對位的方法。

這種方法可運用在掃描版地圖上,我們可以尋找地圖上的座標點當作控制點。

但有的時候,地圖上會沒有這種網格點,又或是你要為影像進行空間對位,這方法也行不通。

在這種情況下,我們可以使用其他的空間對位資料來源來當作輸入。

本教學中,我們會學習如何在空間對位的操作中使用這些開放的資料來源。

內容說明¶ 我們要使用高解析度的氣球影像(Balloon-imagery),以OpenStreetMap的參考座標進行空間對位。

你還會學到這些¶ 下載超高解析度的公領域影像資料 在QGIS中使用OpenLayers附加元件 使用cs2cs指令工具在不同的投影系統之間轉換坐標 使用手邊已經過空間對位的圖層,提供空間對位時的控制點座標資訊 為某圖層自訂「無資料值」(No-dataValue) 取得資料¶ 本教學中,我們要使用ThePublicLaboratory收集的一系列絕美的氣球、風箏空拍影像。

他們其實也有釋出空間對位過的版本,不過為了練習,我們在這裡要下載的是未經空間對位的JPG檔,然後再於QGIS中進行空間對位。

如果你喜歡他們提供的影像,可以在GoogleEarth中進一步探索。

下載紐約市WashingtonSquarePark的JPG影像。

你可以直接在JPG的按鈕上按右鍵然後選Savelinkas…。

為了方便起見,你也可以直接用下面的連結下載: UNTRANSLATEDSTRING newyorkcity-washingtonsquarepark.jpg 操作流程¶ 本教學會使用開放街圖來當作參考圖層。

請到附加元件‣管理與安裝附加元件中安裝OpenLayersPlugin。

有關在QGIS中如何使用附加元件,可參考使用附加元件。

安裝後,就可選擇:menuselection:`附加元件(或網路)-->OpenLayersplugin-->AddOpenStreetMaplayer`_,如此一來QGIS會載入OpenStreetMap預先處理好的資料塊。

現在開放街圖的資料已經讀到QGIS中了。

注意此圖層的座標參考系統(CRS)是EPSG3857-偽麥卡托投影(PseudoMercator),我們接下來在此圖層中使用的座標都會是這個CRS,所以值得在此一提。

下一件事是尋找我們要進行空間對位的地區。

開放街圖的圖層可以自由縮放尋找地點,不過這裏我們來用另外一種未來可能會有用的工具進行搜尋。

我們知道剛剛下載的圖片是位在紐約的WashingtonSquarePark,如果你把這個地點丟到網路上搜尋,就可以在維基百科的頁面中找到此條目,內容還包括公園的座標。

你會發現這裡的座標是經緯度,以度/分/秒的方式顯現,但我們的圖層是麥卡托投影,需要把它轉換成麥卡托直角坐標才能定位。

因此,我們要使用一個方便的指令叫做cs2cs來幫我們這個忙。

一般來說,在Windows安裝QGIS時,應該會把一個稱為OSGEO4W的套件也一併安裝,而在Linux和Mac上亦然。

在後兩個作業系統上,可以進入終端機,輸入cs2cs查看此程式是否存在,而Windows使用者可以去開始(所有程式)‣OSGeo4W‣MSYS啟動終端機(譯按:或是開始(所有程式)‣QGIS‣OSGeo4WShell。

) 備註 Mac的使用者可以在此找到cs2cs工具:/Library/Frameworks/PROJ.Framework/Programs/cs2cs 確定cs2cs指令有安裝在系統中後,就可以著手來把經緯度座標轉換成麥卡托座標。

使用這個工具時,必須要指定初始和目標的CRS,指定的方法可以為PROJ4字串或是EPSG代碼。

因為我們已經知道目標圖層的CRS,所以直接使用後者即可。

此指令最簡單的使用方法是直接在命令列中輸入座標,這個指令讀入座標的順序是XY,所以我們要輸入的是「經度緯度」。

輸入如下所示的指令並按下Enter鍵(指令中的反斜線()是雙引號(」)的逃脫字符。

在某些終端機版本中,不用加上逃脫字符和字串前後的雙引號)。

完成之後,可以看到螢幕上會顯示EPSG3857座標系統的XY座標。

echo"-73d59'51\"40d43'51\""|cs2cs+init=EPSG:4326+to+init=EPSG:3857 -8237364.024972720.340.00 如果以上的指令沒有用,可以試試看以下的不同版本 echo-73d59'5140d43'51|cs2cs+init=EPSG:4326+to+init=EPSG:3857 複製輸出的座標,然後在QGIS視窗最下端的坐標欄位中,把剛才複製的坐標貼上,記得格是要符合X,Y的樣子。

按下Enter鍵之後,可以看到地圖有稍微平移但沒有放大,我們必須要自己調整右側的比例選單,例如調到1:2500後再按下Enter(譯按:如果地圖空白一片,請試試1:5000。

) 太棒了,WashingtonSquarePark出現在畫面中了!總算可以開始來進行空間對位,請選擇影像‣空間對位‣空間對位,如果你找不到此選項的話,要到附加元件‣管理與安裝附加元件‣已安裝的分頁中啟用空間對位GDAL這個附加元件。

在空間對位視窗中選擇檔案‣Openraster,選擇剛下載的JPG檔,然後按開啟。

在選擇座標參考系統視窗中,選擇EPSG:3857PseudoMercator。

接下來按下在工具列上的AddPoint鈕,然後選擇幾個比較容易定位的特徵點,例如轉角、十字路口、雕像等等,這些點稍後要拿來當作控制點。

當你在圖上按下點位時,會有個視窗跳出來要求輸入地圖坐標。

選擇從地圖畫面取得坐標。

在我們用來當作參考圖層的開放街圖上找到相同的位置點一下,地圖上的坐標就會自動被使用為控制點坐標。

使用同樣的方法選擇並記錄至少4個點位的坐標。

接下來選擇設定‣Transformationsettings。

輸入如圖所示的參數設定,記得勾選回到處理完畢後載入QGIS中。

按下確定鈕回到空間對位的視窗中,點選檔案‣Startgeoreferencing,程式就會開始使用控制點轉換產生新影像。

處理程式結束後,已經過空間對位修正的圖層會被載入到QGIS中。

沒意外的話,你會看到新影像與開放街圖的圖層疊合得相當不錯。

如要讓它看起來更順眼,得要移除掉影像周圍的黑色和白色無資料值。

在這個圖層上按右鍵,選擇屬性。

切到透明度分頁,我們要把黑色和白色的像素通通指定為無資料值,就可以讓他們都變透明。

在附加無資料值中輸入0,接著在下方的自訂透明度選項中,按下+按鈕然後在每個波段中都輸入255,透明百分比則輸入100,最後按下確定。

現在,影像看起來就很漂亮的疊在地圖上了。

commentspoweredbyDisqus ThisworkislicensedunderaCreativeCommonsAttribution4.0InternationalLicense



請為這篇文章評分?