I am trying to match rainfall data from USGS/FEWS to coordinate locations of survey households in Ethiopia, so I will need to be able to project the rainfall data. So far, I downloaded the dekadal RFE 2.0 files for East Africa (bulk by year), and for each dekadal, there are five files: .bil, .blw, .clr, .hdr, and .stx. I first added the data (the .bil file) to QGIS, using Africa Albers Equal Area Conic as Projected Coordinate Systems. And then, I opened CSV data of village latitude and longitude and saved it to Africa Albers Equal Area Conic.They looked fine and I could extract point values using point sampling tool.
However, there was wide dispersion from the precipitation data provided by survey. The survey uses the same USGS/FEWS rainfall data as I did, but uses household location rather than village location information. Since dispersion is so wide, I suspect that dispersion came from CRS I chose, rather than household/village level difference. From the USGS/FEWS website, I know following projection information:
Coordinates for corners: Lower left lat : -42.243 deg Lower left lon : -23.490 deg Upper left lat : 43.711 deg Upper left lon : -24.600 deg Lower right lat : -42.242 deg Lower right lon : 63.414 deg Upper right lat : 43.712 deg Upper right lon : 64.523 deg Image size : 1152 rows x 1152 cols Center lat,lon : 1.000000, 20.000000 Pixel size h x w : 8.000000 km x 8.000000 km Origin of latitudes : 1.000000 deg Central meridian : 20.000000 deg First std parallel : -19.000000 deg Second std parallel : 21.000000 deg Projection = ALBERS Conical Equal-area projection uses the clarke 1866 spheroid
Using this information, I am trying to custom projection as follows:
+proj=aea +lat_1=-19 +lat_2=21 +lat_0=1 +lon_0=20 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=clrk66
Can anybody confirm whether this projection is correct? I am beginner of QGIS user and only know the following basic parameters.
- lat_1 : First std parallel
- lat_2 : 2nd std paralell
- lat_0 : Origin of latitudes
- lon_0 : Central meridian
But I am not sure about +datum and +units.