This forum is shutting down! Please post new discussions at community.naturalcapitalproject.org

Water Yield Model Unprojected dataset.

amelineameline Member
edited June 2014 in Freshwater Models
Hello,
I am trying to run the Water Yield model of Invest (version 3.0.0). I am working in a site located in Costa Rica. For the Average Annual Reference Evapotranspiration, I am using the Global Aridity and PET Database from CGIAR (http://www.cgiar-csi.org/data/global-aridity-and-pet-database). This layer has the following attributes :
spatial references: GCS_WGS_1984
Linear Unit : /
angular unit : Degree (0,0174532925199433)
Datum : D_WGS_1984


I encountered the following error while running InVest:

DatasetUnprojected: These datasets are unprojected set(['C:/Documents and Settings/Administrateur/Bureau/Water Yield/INPUT/pet_he_yr_cr.tif'])

I projected the PET layer in WGS_1984_UTM_Zone_17N through ArcGis. Now, the layer has the following attributes : Spatial references : WGS_1984_UTM_Zone_17N
Linear Unit: Meter (1,000000)
Angular Unit : Degree(0,0174532925199433)
False_Easting : 500000
False_Northing : 0
Central_Meridian : 81
Scale_Factor : 0,9996
Latitude_Of_Origin :0
Datum : D_WGS_1984

When I use it in the Water Yield Model, the following error occures :

Exception: The datasets' intersection is empty (i.e., not all the datasets touch each other).
None


The log file indicates :

[...]
06/12/2014 15:30:41 hydropower_water_yield DEBUG Calculate PET from Ref Evap times Kc
06/12/2014 15:30:41 raster_utils WARNING These two datasets are not in the same projection. The different projections are:

'filename: C:/Documents and Settings/Administrateur/Bureau/Water Yield/INPUT/pet_project.tif'
PROJCS["WGS_1984_UTM_Zone_17N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-81],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32617"]]

and:

'filename:c:\documents and settings\administrateur\bureau\water yield\tmp\tmpy3ortb'
PROJCS["WGS_1984_Transverse_Mercator",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-84],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]


06/12/2014 15:30:41 root ERROR ---------------------------------------------------
06/12/2014 15:30:41 root ERROR ---------------------- ERROR ----------------------
06/12/2014 15:30:41 root ERROR ---------------------------------------------------
06/12/2014 15:30:41 root ERROR Error: exception found while running hydropower_water_yield
06/12/2014 15:30:41 root DEBUG

[...]

Exception: The datasets' intersection is empty (i.e., not all the datasets touch each other).
None

06/12/2014 15:30:41 root ERROR Exiting due to failures


I don't understand what I am doing wrong... If someone can help me, I would be very grateful to him !
Ameline,

Comments

  • same errors with the new (3.0.1) version of InVest
  • swolnyswolny Member, NatCap Staff
    edited June 2014
    Hi Ameline -

    Do your other input layers also have the same projection as the PET layer (WGS_1984_UTM_Zone_17N)? If not, choose one to standardize on (UTM 17N is good) and reproject them all to match. In ArcGIS it's easy to pick a specific layer's coordinate system in to match when using the Project Raster tool - in the Output Coordinate System input window, expand the Layers drop-down and you should see the ones used by the layers in your .mxd.

    If all of your layers are in the same projection already, if you look at each of them in turn in ArcGIS, do they overlap, or is one off on its own?

    ~ Stacie
    Post edited by swolny on
  • Hi Stacie,

    Thank you so much ! One of my layer was in WGS_1984_Transverse_Mercator. I projected it through ArcGis in UTM 17N, and now, the model is running :) It sounds like a beginner mistake...
This discussion has been closed.