Dapscoptyltd / QGIS

All things related to QGIS in here.
0 stars 1 forks source link

Technical Log: Raster Map Overlay Preparation #10

Open Dapscoptyltd opened 7 years ago

Dapscoptyltd commented 7 years ago

Intro

I needed to prepare a raster map to overlay onto my sample map project, to confirm it was possible and that I could do it. This meant preparing an image file from a map file.

Upon conferring with Dr. Adela Subotkova, I confirmed that the printed map should be converted into an image file.

Confirming my OSM Map was linked to the correct CRS

To ensure that any raster file imported to my map was correctly projected, I needed to confirm what my existing CRS was, and that it was correctly chosen for my map and the region. I had downloaded, as my first map a source from:

Properties
Title: OSM Standard
Credit: © OpenStreetMap contributors, CC-BY-SA
URL: http://tile.openstreetmap.org/{z}/{x}/{y}.png

The CRS chosen for the OSM street map was WGS 84 / World Mercator.

The project properties window: screen shot 2017-06-01 at 17 59 51

Testing the coordinates

To obtain a correct North and East coordinate for testing my OSM map and CRS choices in QGIS, I chose the intersection of Langemarkstraat, Ieperstraat, and Roeselarestraat in Zonnebeke. The exact point chosen was the rendered "centre" of the roundabout. This intersection is in front of the Onze-Lieve-Vrouwkerk, another prominent landmark on maps. screen shot 2017-06-01 at 17 53 47

Step One

Using the guide found here, I obtained northing and easting points from Google Maps: 50°52'21.8"N 2°59'15.7"E 50.872731, 2.987694. This is the N, E latitude and longitude coordinates format from the datum location. screen shot 2017-06-01 at 17 54 19 The results are shown in the figure above.

The latitude and longitude coordinates are of the most interest and confirm in my mind that this was an accurate (certainly for the purposes of this exercise). 2° is very close to Greenwich, and its positive number showed it as an easterly direction. The 50° bearing for latitude also is indicative that this bearing was accurate.

Step Two

  1. I confirmed my project settings in QGIS for map and CRS were as above.
  2. On my OSM map layer, I zoomed into the intersection of Langemarkstraat, Ieperstraat, and Roeselarestraat in Zonnebeke, to the rendered "centre" of the roundabout.
  3. I selected the OSM Tools button.
  4. I clicked the Map Coordinate. -This turned the cursor into a target.
  5. I selected the centre of the roundabout.
    • The From Map values from my selection with the cursor target were: -- 2.98770, 50.87272, which is an 'x', 'y' coordinate, where 'x' = North and 'y' = East.
  6. The result is shown below.

screen shot 2017-06-01 at 17 56 44

Step Three

  1. I went to the eps.io website.
  2. I performed a search on Belgium, Zonnebeke.
  3. I clicked through the results screens.
  4. The coordinate results were delivered to be:
    • 2.987481 'x' or East.
    • 50.872722 'y' or North.
    • This conformed to the results from above.

While this is far from an exhaustive test, I accept that for the purposes of this exercise, this is enough to confirm the CRS for this sample project's OSM map are correct.

Further Image creation Research

I made a few searches in Google to ascertain hidden problems:

This forum response proved helpful.

Prepare the Raster map tif file

This map needed to be, as confirmed by Dr. Sobotkova, converted into a raster graphic I chose a tif file without progressive compression. This map also needed information removed. The conversion was thus chosen to be performed in Adobe PhotoShop, but many graphic software packages could be used.

  1. In Adobe Photoshop I opened the file "1917-10-10_19-00_D11toK2EA_SitMap_RCDIG1014546-4.pdf" (for more information refer to my #9 Technical log). screen shot 2017-06-01 at 18 50 53
  2. I saved it as "1917-10-10_19-00_D11toK2EA_SitMap_RCDIG1014546-4.tif.
  3. I unlocked the background layer. screen shot 2017-06-01 at 18 53 49 screen shot 2017-06-01 at 18 54 01
  4. Next I needed to remove the extremities of the map and the legend.
    • This is done with the Lasso tool.
    • It has a number of options.
    • I chose the Polygonal Lasso tool, because it is best for straight edge trimming like this.
  5. I entered a cut line. screen shot 2017-06-01 at 18 57 38
  6. I selected Layer > Inverse.
  7. I then deleted the selected sections.
  8. Save as 1917-10-10_19-00_D11toK2EA_SitMap_RCDIG1014546-4.tif
    • This file had a large white space section in the top left as a result of removing the legend.
  9. I cropped the file and saved that as "1917-10-10_19-00_D11toK2EA_SitMap_RCDIG1014546-4_a.tif"
  10. I opened the original tif and cropped that image, and saved as "1917-10-10_19-00_D11toK2EA_SitMap_RCDIG1014546-4_b.tif"

These files are attached to #9 to explain file saving conventions.