[Qgis-user] 1928-1930 USGS maps (Nicolas Cadieux)
Nicolas Cadieux
nicolas.cadieux at archeotec.ca
Mon Jan 4 08:04:15 PST 2016
Hi,
I looked at your files. The borderless map are 8 byte file but they are
Boolean value files 0 and 1 only (black and white file. No gray scale).
In this case, 0 = data and 1 = no data??? (use the value plugin to see
this). So if you set your 1 values to no data (layer
properties/transparency/additional no data value), and set your layer
style from black to white (and not the default white to black) you will
see that your data is well georeferenced. (use an open layer layer like
osm or Google map as background.)
The files with borders (honolulu1928.tif) are Eight bit unsigned integer
but they are on a gray scale so the vary between 0 and 255. If you
look at the geotagg, (raster/miscellaneous/info (see below) that the
file has some georeferencing info like the crs but the coordinate are
wrong (Corner Coordinates:
Upper Left ( 0.0, 0.0)). So if you turn on the "CRS on the fly"
and then zoom to the layer extent, you will see the map in the sea.
It's not that the map does not appear... it's that it is in the wrong spot.
I also see this tag (TIFFTAG_SOFTWARE=Adobe Photoshop CS Windows). It is
possible that Photoshp was used in the past and that it striped away the
coordinates. Since you have the coordinates on the map, you can use the
gdal georeferencing plugin to place the map in the correct spot.
Good luck
Nicolas
==============
Driver: GTiff/GeoTIFF
Files: C:/Users/Nicolas/Downloads/1928 Honolulu/honolulu1928.tif
C:/Users/Nicolas/Downloads/1928 Honolulu/honolulu1928.tif.aux.xml
Size is 12000, 10519
Coordinate System is:
PROJCS["NAD_1983_HARN_UTM_Zone_4N",
GEOGCS["GCS_North_American_1983_HARN",
DATUM["NAD83_High_Accuracy_Reference_Network",
SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",500000.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",-159.0],
PARAMETER["Scale_Factor",0.9996],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Meter",1.0],
AUTHORITY["EPSG","3750"]]
Metadata:
TIFFTAG_DATETIME=2006:06:15 07:13:01
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
*_TIFFTAG_SOFTWARE=Adobe Photoshop CS Windows_*
TIFFTAG_XRESOLUTION=400
TIFFTAG_YRESOLUTION=400
Image Structure Metadata:
INTERLEAVE=BAND
_*Corner Coordinates:*__*
*__*Upper Left ( 0.0, 0.0)*__*
*__*Lower Left ( 0.0,10519.0)*__*
*__*Upper Right (12000.0, 0.0)*__*
*__*Lower Right (12000.0,10519.0)*_
Center ( 6000.0, 5259.5)
Band 1 Block=12000x1 Type=Byte, ColorInterp=Gray
Min=0.000 Max=255.000
Minimum=0.000, Maximum=255.000, Mean=193.512, StdDev=41.612
Metadata:
LAYER_TYPE=athematic
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=193.51169711538
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=41.611752735755
====================
Le 2016-01-04 03:04, Jan Becket [via OSGeo.org] a écrit :
> Yes, the min/max values are between 0 and 1. This is strange, because
> when I open the file in Photoshop (and then exit without saving), it
> is an 8-bit grayscale. I actually have two sets of 1928-30 USGS
> geotiffs for the island of Oahu. The set I was able to mount contains
> 16 geotiff files for each of the 16 quads that are 4-5 mb in size. The
> other set contains 16 geotiff files that are about 125 mb each, and
> clearly grayscale. These USGS quads, however, have collars and are
> less desirable for that reason. In any case, they are not visible when
> I attempt to import / mount them. In case anyone has a suggestion, I
> have placed the Honolulu quad for both sets of maps in this dropbox
> folder:
> https://www.dropbox.com/sh/te9et5yuery207d/AAAYnBBkBzUod53RauD7ArTda?dl=0
>
> As far as I can tell, the CRS for both sets of maps is Old Hawaiian. I
> think I mentioned earlier that I am reprojecting to CRS 3711: NAD 83,
> UTM zone 4N.
>
> Thanks, as always, for suggestions -
> Jan
>
>
>
>
>> Message: 1
>> Date: Sat, 2 Jan 2016 15:58:12 -0800 (PST)
>> From: Nicolas Cadieux <[hidden email]
>> </user/SendEmail.jtp?type=node&node=5243557&i=0>>
>> To: [hidden email] </user/SendEmail.jtp?type=node&node=5243557&i=1>
>> Subject: Re: [Qgis-user] 1928-1930 USGS maps
>> Message-ID: <[hidden email]
>> </user/SendEmail.jtp?type=node&node=5243557&i=2>>
>> Content-Type: text/plain; charset="utf-8"
>>
>> Good but weird.
>> If you look at the histogram, are all values between 0 and 1?
>> Nicolas Cadieux M.Sc.
>> Les Entreprises Archéotec inc.
>> 8548, rue Saint-Denis Montréal H2P 2H2
>> Téléphone: 514.381.5112 Fax: 514.381.4995
>> www.archeotec.ca <http://www.archeotec.ca>
>> On Dec 28, 2015 17:29, "William Kyngesburye [via OSGeo.org
>> <http://OSGeo.org>]" <[hidden email]
>> </user/SendEmail.jtp?type=node&node=5243557&i=3>> wrote:
>>
>> Max 1? Odd for a greyscale image. Or are they 1bit images?
>> At least it's working.
>> > On Dec 28, 2015, at 4:12 PM, Jan Becket < [hidden email] >
>> wrote:
>> >
>> > With your suggestions and those of Nicolas I was able to make
>> the maps appear. I trashed the copies I had altered in Photoshop,
>> loaded the original (black) images, set the min as 0 and the max as 1
>> and then stretched to min/max. That turned out to be easy. Thanks for
>> the help!
>> >
>> > - Jan
>> >
>> >
>> >> On Dec 28, 2015, at 4:42 AM, William Kyngesburye <
>> [hidden email] > wrote:
>> >>
>> >> Or, QGIS should be able to read 16bit greyscale images
>> without breaking them with Photoshop. You just need to set a min and
>> max to include the 16bit range and stretch the contrast to minmax.
>> Make sure you are not using QGIS 2.10 - that had a broken raster
>> styling on OS X.
>> >>
>> >> Alternately, use GDAL to convert to 8bit, not Photoshop.
>> "gdal_translate -scale -ot Byte" should do it (in Terminal,
>> I don't see the -scale or -ot options in QGIS GDAL convert format).
>> >>
>> >>> On Dec 28, 2015, at 12:53 AM, Nicolas Cadieux <
>> [hidden email] > wrote:
>> >>>
>> >>> Hi,
>> >>> You are correct, you have stripes the CRS info from the
>> map because Photoshop does not know what to do. If you load the map,
>> you can click on raster.. Something... Gdal-info. That should give
>> you all the tagged info left in the image.
>> >>>
>> >>> Load the map in QGIS and click on "zoom to layer
>> extent". If you still don't see it, make sure you are in a
>> CRS that is projected in meters. You should be able to see it. The
>> map should be at 0,0.
>> >>>
>> >>> All is not loss... You can reproject the map to the
>> correct CRS using the gdal georeferencing tool. Always back up
>> your data before manipulating it.
>> >>> Hope this helps!
>> >>>
>> >>> Nicolas Cadieux M.Sc.
>> >>> Les Entreprises Archéotec inc.
>> >>> 8548, rue Saint-Denis Montréal H2P 2H2
>> >>> Téléphone: 514.381.5112 Fax: 514.381.4995
>> >>> www.archeotec.ca <http://www.archeotec.ca>
>> >>>
>> >>> On Dec 28, 2015 01:16, "Jan Becket [via OSGeo.org
>> <http://OSGeo.org>]" <[hidden email]> wrote:
>> >>> After much effort, I obtained an old set of 1928-1930
>> USGS topo maps of Oʻahu Island, Hawaiʻi - for a project I am working
>> on. The maps are without collars and used by a local archaeological
>> firm. The geotiff maps at first showed as completely black on my Mac
>> system (10.10.5 - Yosemite). However I found this link with a
>> resultuion for the issue:
>> >>>
>> >>>
>> http://stackoverflow.com/questions/13925995/options-to-convert-16-bit-image
>> >>>
>> >>> I opened the 16 quads for Oʻahu in Photoshop and applied
>> auto levels, after which they became discernable quad maps. However,
>> when imported into QGIS, they are now invisible. I wonder if it is
>> possible that I stripped away the geotiff CRS information. Prior to
>> the Photoshop operation, the maps did appear correctly referenced to
>> Oʻahu Island (although they were totally opaque / black).
>> >>>
>> >>> Now, the maps in QGIS are not visible, although the
>> geotiffs are several megabytes in size and are visible when opened in
>> Photoshop, Preview, or whatever image viewer.
>> >>>
>> >>> I should mention that I have successfully imported other
>> sets of USGS maps of Oʻahu for this project - the 1902 series, the
>> 1916 series, the 1935 series - as well as the current series. The
>> early maps used the Old Hawaiian CRS (EPSG 3564) but more modern maps
>> use WGS83, Zone 4N (EPSG 3711). I have on the fly projection enabled
>> and have tried many CRS settings for the maps, but nothing results in
>> a visible quad appearing.
>> >>>
>> >>> Suggestions very welcome ...
>> -----
>> William Kyngesburye <kyngchaos*at*kyngchaos*dot*com>
>> http://www.kyngchaos.com/ "This is a question about the past, is
>> it? ... How can I tell that the past isn't a fiction designed to
>> account for the discrepancy between my immediate physical sensations
>> and my state of mind?"
>> - The Ruler of the Universe
>> _______________________________________________
>
>
> _______________________________________________
> Qgis-user mailing list
> [hidden email] </user/SendEmail.jtp?type=node&node=5243557&i=4>
> List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
> Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> ------------------------------------------------------------------------
> If you reply to this email, your message will be added to the
> discussion below:
> http://osgeo-org.1560.x6.nabble.com/Re-1928-1930-USGS-maps-Nicolas-Cadieux-tp5243557.html
>
> To start a new topic under Quantum GIS - User, email
> ml-node+s1560n4125267h38 at n6.nabble.com
> To unsubscribe from Quantum GIS - User, click here
> <http://osgeo-org.1560.x6.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4125267&code=bmljb2xhcy5jYWRpZXV4QGFyY2hlb3RlYy5jYXw0MTI1MjY3fDYzNDQ4MjQxNg==>.
> NAML
> <http://osgeo-org.1560.x6.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Re-1928-1930-USGS-maps-Nicolas-Cadieux-tp5243557p5243638.html
Sent from the Quantum GIS - User mailing list archive at Nabble.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-user/attachments/20160104/41c5ac69/attachment.html>
More information about the Qgis-user
mailing list