PostGIS  2.1.10dev-r@@SVN_REVISION@@
static void test_gdal_warp ( )
static

Definition at line 584 of file cu_gdal.c.

References ovdump::band, cu_add_band(), cu_free_raster(), ES_NONE, PT_64BF, rtpixdump::rast, rtrowdump::raster, rt_band_get_hasnodata_flag(), rt_band_get_nodata(), rt_band_get_pixel(), rt_band_set_pixel(), rt_raster_gdal_warp(), rt_raster_get_band(), rt_raster_get_height(), rt_raster_get_num_bands(), rt_raster_get_width(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), genraster::value, pixval::x, and pixval::y.

Referenced by gdal_suite_setup().

584  {
585  rt_pixtype pixtype = PT_64BF;
586  rt_band band = NULL;
587 
589  rt_raster rast;
590  uint32_t x;
591  uint32_t width = 100;
592  uint32_t y;
593  uint32_t height = 100;
594  double value = 0;
595 
596  char src_srs[] = "PROJCS[\"unnamed\",GEOGCS[\"unnamed ellipse\",DATUM[\"unknown\",SPHEROID[\"unnamed\",6370997,0]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_center\",45],PARAMETER[\"longitude_of_center\",-100],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1],AUTHORITY[\"EPSG\",\"2163\"]]";
597 
598  char dst_srs[] = "PROJCS[\"NAD83 / California Albers\",GEOGCS[\"NAD83\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"6269\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4269\"]],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"standard_parallel_1\",34],PARAMETER[\"standard_parallel_2\",40.5],PARAMETER[\"latitude_of_center\",0],PARAMETER[\"longitude_of_center\",-120],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",-4000000],AUTHORITY[\"EPSG\",\"3310\"],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]]";
599 
600  raster = rt_raster_new(width, height);
601  CU_ASSERT(raster != NULL); /* or we're out of virtual memory */
602 
603  band = cu_add_band(raster, pixtype, 1, 0);
604  CU_ASSERT(band != NULL);
605 
606  rt_raster_set_offsets(raster, -500000, 600000);
607  rt_raster_set_scale(raster, 1000, -1000);
608 
609  for (x = 0; x < width; x++) {
610  for (y = 0; y < height; y++) {
611  rt_band_set_pixel(band, x, y, (((double) x * y) + (x + y) + (x + y * x)) / (x + y + 1), NULL);
612  }
613  }
614 
615  rast = rt_raster_gdal_warp(
616  raster,
617  src_srs, dst_srs,
618  NULL, NULL,
619  NULL, NULL,
620  NULL, NULL,
621  NULL, NULL,
622  NULL, NULL,
623  GRA_NearestNeighbour, -1
624  );
625  CU_ASSERT(rast != NULL);
626  CU_ASSERT_EQUAL(rt_raster_get_width(rast), 122);
627  CU_ASSERT_EQUAL(rt_raster_get_height(rast), 116);
628  CU_ASSERT_NOT_EQUAL(rt_raster_get_num_bands(rast), 0);
629 
630  band = rt_raster_get_band(rast, 0);
631  CU_ASSERT(band != NULL);
632 
633  CU_ASSERT(rt_band_get_hasnodata_flag(band));
634  rt_band_get_nodata(band, &value);
635  CU_ASSERT_DOUBLE_EQUAL(value, 0., DBL_EPSILON);
636 
637  CU_ASSERT_EQUAL(rt_band_get_pixel(band, 0, 0, &value, NULL), ES_NONE);
638  CU_ASSERT_DOUBLE_EQUAL(value, 0., DBL_EPSILON);
639 
640  cu_free_raster(rast);
641  cu_free_raster(raster);
642 }
int rt_raster_get_num_bands(rt_raster raster)
Definition: rt_api.c:5677
rt_raster rt_raster_gdal_warp(rt_raster raster, const char *src_srs, const char *dst_srs, double *scale_x, double *scale_y, int *width, int *height, double *ul_xw, double *ul_yw, double *grid_xw, double *grid_yw, double *skew_x, double *skew_y, GDALResampleAlg resample_alg, double max_err)
Return a warped raster using GDAL Warp API.
Definition: rt_api.c:9710
tuple band
Definition: ovdump.py:57
tuple rast
Definition: rtpixdump.py:62
tuple raster
Be careful!! Zeros function's input parameter can be a (height x width) array, not (width x height): ...
Definition: rtrowdump.py:123
void cu_free_raster(rt_raster raster)
rt_band cu_add_band(rt_raster raster, rt_pixtype pixtype, int hasnodata, double nodataval)
void rt_raster_set_scale(rt_raster raster, double scaleX, double scaleY)
Set scale in projection units.
Definition: rt_api.c:5442
rt_pixtype
Definition: rt_api.h:172
rt_errorstate rt_band_get_nodata(rt_band band, double *nodata)
Get NODATA value.
Definition: rt_api.c:3058
rt_errorstate rt_band_get_pixel(rt_band band, int x, int y, double *value, int *nodata)
Get pixel value.
Definition: rt_api.c:2549
tuple x
Definition: pixval.py:53
uint16_t rt_raster_get_height(rt_raster raster)
Definition: rt_api.c:5434
rt_band rt_raster_get_band(rt_raster raster, int n)
Return Nth band, or NULL if unavailable.
Definition: rt_api.c:5686
rt_raster rt_raster_new(uint32_t width, uint32_t height)
Construct a raster with given dimensions.
Definition: rt_api.c:5353
int rt_band_get_hasnodata_flag(rt_band band)
Get hasnodata flag value.
Definition: rt_api.c:2002
void rt_raster_set_offsets(rt_raster raster, double x, double y)
Set insertion points in projection units.
Definition: rt_api.c:5504
tuple y
Definition: pixval.py:54
uint16_t rt_raster_get_width(rt_raster raster)
Definition: rt_api.c:5426
rt_errorstate rt_band_set_pixel(rt_band band, int x, int y, double val, int *converted)
Set single pixel's value.
Definition: rt_api.c:2302

Here is the call graph for this function:

Here is the caller graph for this function: