<div dir="ltr">Trying this again. Thanks, Mike Flanagan, for trying to help me the first time.<div><br></div><div>I start a new project in QGIS 3.12. I have a USGS map CA_Concord_297163_1915_62500_geo.tif, downloaded from <a href="https://ngmdb.usgs.gov/topoview/viewer/#11/37.8612/-122.0759" target="_blank">here</a>, the "Concord, CA 1915, 1915 edition" map. I add "Google Maps" from Browser > XYZ Tiles as a base map. My project properties show WGS/84 Pseudo-Mercator, EPSG:3857. I select Layer > Add Layer > Add Raster Layer... > I can select this geotiff. Clicking Add brings up this dialog. What do I choose?</div><div><br></div><div><div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:#f5f5f5;padding:5px;color:#222;font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid #ddd;line-height:1"><a href="https://drive.google.com/file/d/1FKVPrHLtXXKK8a-QTTz3IigT9dpj1Hou/view?usp=drive_web" target="_blank" style="display:inline-block;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:none;width:100%"><img style="vertical-align: bottom; border: none;" src="https://drive-thirdparty.googleusercontent.com/16/type/image/png"> <span dir="ltr" style="color:#15c;text-decoration:none;vertical-align:bottom">Screen Shot 2020-04-18 at 12.54.31 PM.png</span></a></div><br></div><div>What do I choose?</div><div><br></div><div>The reason why I ask, is that all my usgs maps are offset from the base map enough to be maddening. I know these things are imperfect, but even the corners labeled with known points are off. Below is the gdalinfo output for this file. Last question; if I georeference this tif file, versus trusting its metadata, the top corners left and right have two divergent borders. They start from a common point at bottom and at the top they are maybe a degree apart. Which of these two endpoints matches the stated coordinate?</div><div><br></div><div><font face="monospace">Morgans-MBP-2:CA_Concord_297163_1915_62500_geo_tif morganfletcher$ /Applications/QGIS3.12.app/Contents/MacOS/bin/gdalinfo CA_Concord_297163_1915_62500_geo.tif<br>Driver: GTiff/GeoTIFF<br>Files: CA_Concord_297163_1915_62500_geo.tif<br>Size is 4775, 5908<br>Coordinate System is:<br>PROJCS["unnamed",<br>    GEOGCS["NAD27",<br>        DATUM["North_American_Datum_1927",<br>            SPHEROID["Clarke 1866",6378206.4,294.9786982138982,<br>                AUTHORITY["EPSG","7008"]],<br>            AUTHORITY["EPSG","6267"]],<br>        PRIMEM["Greenwich",0],<br>        UNIT["degree",0.0174532925199433],<br>        AUTHORITY["EPSG","4267"]],<br>    PROJECTION["Polyconic"],<br>    PARAMETER["latitude_of_origin",0],<br>    PARAMETER["central_meridian",-122.125],<br>    PARAMETER["false_easting",0],<br>    PARAMETER["false_northing",0],<br>    UNIT["metre",1,<br>        AUTHORITY["EPSG","9001"]]]<br>Origin = (-12632.067326823029362,4208868.700336153618991)<br>Pixel Size = (5.291666666666666,-5.291666666666640)<br>Metadata:<br>  AREA_OR_POINT=Area<br>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br>  TIFFTAG_XRESOLUTION=300<br>  TIFFTAG_YRESOLUTION=300<br>Image Structure Metadata:<br>  COMPRESSION=YCbCr JPEG<br>  INTERLEAVE=PIXEL<br>  SOURCE_COLOR_SPACE=YCbCr<br>Corner Coordinates:<br>Upper Left  (  -12632.067, 4208868.700) (122d16' 7.84"W, 38d 0'50.81"N)<br>Lower Left  (  -12632.067, 4177605.534) (122d16' 5.87"W, 37d43'56.80"N)<br>Upper Right (   12635.641, 4208868.700) (121d58'52.01"W, 38d 0'50.81"N)<br>Lower Right (   12635.641, 4177605.534) (121d58'53.98"W, 37d43'56.80"N)<br>Center      (       1.787, 4193237.117) (122d 7'29.93"W, 37d52'24.13"N)<br>Band 1 Block=512x512 Type=Byte, ColorInterp=Red<br>  Overviews: 2388x2954, 1194x1477, 597x739, 299x370, 150x185, 75x93, 38x47<br>Band 2 Block=512x512 Type=Byte, ColorInterp=Green<br>  Overviews: 2388x2954, 1194x1477, 597x739, 299x370, 150x185, 75x93, 38x47<br>Band 3 Block=512x512 Type=Byte, ColorInterp=Blue<br>  Overviews: 2388x2954, 1194x1477, 597x739, 299x370, 150x185, 75x93, 38x47</font><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Apr 18, 2020 at 2:43 PM Morgan Fletcher <<a href="mailto:morgan@hahaha.org" target="_blank">morgan@hahaha.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">(Resending, with screenshots as google drive links, to cut down message size.)<br><br>I am running 3.12.1-București on a mac, os-x 10.14.6. My knowledge of GIS is rudimentary, I've learned what I need to, to satisfy my curiosity about my region, using old maps.<br> <br>I have downloaded all the USGS maps for my region - Oakland, CA, USA, which corresponds to "CA Concord" - from USGS topoView, in geotiff format. Opening them in QGIS gives fairly wild offsets; they are not all encoded the same way, yet have common corners. Features on the map vary with what modern maps show, but getting their corners right would help me a lot.<br><br>If I create a new project, set my CRS to WGS 84 / Pseudo-Mercator (EPSG:3857) as I think is correct for my base map, Google Maps (EPSG:3857 - WGS 84 / Pseudo-Mercator - Projected) then Layer > Add Layer > Add Raster Layer > Raster Dataset > CA_Concord_465520_1943_62500_geo.tif > Add, I get prompted with:<br><br><a href="https://drive.google.com/file/d/1FKVPrHLtXXKK8a-QTTz3IigT9dpj1Hou/view?usp=sharing" target="_blank">https://drive.google.com/file/d/1FKVPrHLtXXKK8a-QTTz3IigT9dpj1Hou/view?usp=sharing<br></a><br>I don't know what to choose, here. I chose the first one. I've tried others. The tif.prj file does say <font face="monospace">DATUM["D_North_American_1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]]</font>, which I am guessing conforms to NAD27. If I add an opacity slider to the new raster layer, set its opacity, and then use the Zoom to Coordinate panel to navigate to <font face="monospace">37º45'N 122º15'W</font>, I find that the point is off from what's on the map, see:<br><br><a href="https://drive.google.com/file/d/10vYtcM1FZ-vdS3z1FPW2EYbrNy2sj65L/view?usp=sharing" target="_blank">https://drive.google.com/file/d/10vYtcM1FZ-vdS3z1FPW2EYbrNy2sj65L/view?usp=sharing<br></a><br>This is similarly true at the top left:<br><br><a href="https://drive.google.com/file/d/18naBKuwk6rZ9DP0tFOZxnGKlIpRQgDKM/view?usp=sharing" target="_blank">https://drive.google.com/file/d/18naBKuwk6rZ9DP0tFOZxnGKlIpRQgDKM/view?usp=sharing<br></a><br>Also, note how the map shows two lines at top left, it's true at top right. This is something to do with projection, I think. I've tried to georeference these lines, to get them to line up right, but which intersection do I use, when there are two vertical lines intersecting with the top, horizontal border, to reference <font face="monospace">122º 15' 38º 00'</font> ? I've had mixed results, trying to use deg-min-ss coordinates to georeference these; I get a map off the west coast of Africa. Using points on the map is more successful. <br><br>All the maps have slightly variant values in their tfw files. Their prj files are more uniform, but also show variations. I'd like to normalize them all, so that I can view the USGS topo maps for my region over time. Advice?<br><br>Thanks,<br><br>Morgan</div>
</blockquote></div>