I have downloaded a raster, which has the following information.
Item Data Date 2015-01-03 00:00:00 Product cold partial sst Satellite NOAA-18 Direction Unknown Dimension 1596 x 1617 BPP 8 Latitude range 47 — 62.999108 Longditude range -15 — 13 Real → DN dn = (rw - -3) / 0.1 DN → real rw = 0.1 * dn + -3 DN range 0 — 255 Pixel dimensions 1.100149 x 1.100149
On Arcmap it shows the pixel X/Y co-ordinates. How can I convert this to Long and Lat?
From the FAQ of the site that I downloaded the raster from the following information is provided. link
How do I convert image coordinates to latitude/longitude? For converting from x/y pixel co-ordinates on a Mercator projection image to latitude/longitude (lat/lon) positions use:
rows - number of rows in image cols - number of columns in image minlon/maxlon - minimum/maximum longitude of image (in decimal degrees) minlat/maxlat - minimum/maximum latitude of image (in decimal degrees) Note that these are the coordinates of the CENTRE of the corner pixels. DEGTORAD - conversion from degrees to radians (PI/180.0) RADTODEG - conversion from radians to degrees to radians ln - natural log Calculate longitude. lonfract = x / (cols - 1) lon = minlon + (lonfract * (maxlon - minlon)) Calculate latitude. latfract = 1.0 - (y / (rows - 1)) Ymin = ln (tan (DEGTORAD * (45.0 + (minlat / 2.0)))) Ymax = ln (tan (DEGTORAD * (45.0 + (maxlat / 2.0)))) Yint = Ymin + (latfract * (Ymax - Ymin)) lat = 2.0 * (RADTODEG * (arctan (exp (Yint))) - 45.0)
The maximum/minimum latitude/longitude values and sizes are given when the images are viewed.
Would it possible to run the above code in ArcMAP or R (raster package)?