Linked below is a polygon shapefile whose projection is NAD83 and whose values are LU/LC codes. Coordinates are meters of northing and easting. I would like to convert it to an ASCII raster with coordinates measured in latitude and longitude. I'm sure it is a simple process, but I can't figure out the re-projection. Tools I have are ArcGIS or SAGA or R. I can even open the ASCII text and modify the header if no recalculation is necessary.
http://hawaii.gov/dbedt/gis/lulc.htm
EDIT: Thanks to Dan Patterson for pointing this out: the projection is actually UTM 4, and NAD83 is the DATUM. The output raster should have the same datum. If I want output to have lat/long coordinates, does that mean it has no projection?
Answer
You have two tasks here. The first is to to coordinate transform (re-project) from UTM to Lat/Long (also called geographic or "no projection") and the second to convert from vector to raster.
In arcgis one approach is to use the Project tool to convert to geographic, then use one the To Raster geoprocessing tools located in the Conversion Tools toolbox, in this case Polygon to Raster. And now finally Raster to ASCII tool.
A second and perhaps more intutive vector to raster conversion is to load the shapefile into Arcmap, symbolize the polygons with "categories > unique values" as you want them to look, set the data frame coordinate system to Geographic NAD83 or WGS-1984, and then use File > Export map, save as TIFF with Format>"write geotiff tags" checked, bump the dpi and/or pixel dimensions up as you need them, then follow up with RasterToAscii.
No comments:
Post a Comment