PageRenderTime 66ms CodeModel.GetById 28ms RepoModel.GetById 1ms app.codeStats 0ms

/raster/rt_pg/rtpg_spatial_relationship.c

https://github.com/therealtakeshi/postgis
C | 1315 lines | 1051 code | 131 blank | 133 comment | 245 complexity | d83e8082e62fd78e04040c137fb14af0 MD5 | raw file
Possible License(s): GPL-2.0, LGPL-2.1, BSD-3-Clause 
    

  1. /*
    2. 

  2.  * $Id$
    3. 

  3.  *
    4. 

  4.  * WKTRaster - Raster Types for PostGIS
    5. 

  5.  * http://www.postgis.org/support/wiki/index.php?WKTRasterHomePage
    6. 

  6.  *
    7. 

  7.  * Copyright (C) 2011-2013 Regents of the University of California
    8. 

  8.  *   <bkpark@ucdavis.edu>
    9. 

  9.  * Copyright (C) 2010-2011 Jorge Arevalo <jorge.arevalo@deimos-space.com>
    10. 

  10.  * Copyright (C) 2010-2011 David Zwarg <dzwarg@azavea.com>
    11. 

  11.  * Copyright (C) 2009-2011 Pierre Racine <pierre.racine@sbf.ulaval.ca>
    12. 

  12.  * Copyright (C) 2009-2011 Mateusz Loskot <mateusz@loskot.net>
    13. 

  13.  * Copyright (C) 2008-2009 Sandro Santilli <strk@keybit.net>
    14. 

  14.  *
    15. 

  15.  * This program is free software; you can redistribute it and/or
    16. 

  16.  * modify it under the terms of the GNU General Public License
    17. 

  17.  * as published by the Free Software Foundation; either version 2
    18. 

  18.  * of the License, or (at your option) any later version.
    19. 

  19.  *
    20. 

  20.  * This program is distributed in the hope that it will be useful,
    21. 

  21.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
    22. 

  22.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    23. 

  23.  * GNU General Public License for more details.
    24. 

  24.  *
    25. 

  25.  * You should have received a copy of the GNU General Public License
    26. 

  26.  * along with this program; if not, write to the Free Software Foundation,
    27. 

  27.  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
    28. 

  28.  *
    29. 

  29.  */
    30. 

  30. 

  31. #include <postgres.h> /* for palloc */
    32. 

  32. #include <fmgr.h>
    33. 

  33. 

  34. #include "../../postgis_config.h"
    35. 

  35. 

  36. #include "lwgeom_pg.h"
    37. 

  37. #include "rtpostgis.h"
    38. 

  38. 

  39. /* determine if two rasters intersect */
    40. 

  40. Datum RASTER_intersects(PG_FUNCTION_ARGS);
    41. 

  41. 

  42. /* determine if two rasters overlap */
    43. 

  43. Datum RASTER_overlaps(PG_FUNCTION_ARGS);
    44. 

  44. 

  45. /* determine if two rasters touch */
    46. 

  46. Datum RASTER_touches(PG_FUNCTION_ARGS);
    47. 

  47. 

  48. /* determine if the first raster contains the second raster */
    49. 

  49. Datum RASTER_contains(PG_FUNCTION_ARGS);
    50. 

  50. 

  51. /* determine if the first raster contains properly the second raster */
    52. 

  52. Datum RASTER_containsProperly(PG_FUNCTION_ARGS);
    53. 

  53. 

  54. /* determine if the first raster covers the second raster */
    55. 

  55. Datum RASTER_covers(PG_FUNCTION_ARGS);
    56. 

  56. 

  57. /* determine if the first raster is covered by the second raster */
    58. 

  58. Datum RASTER_coveredby(PG_FUNCTION_ARGS);
    59. 

  59. 

  60. /* determine if the two rasters are within the specified distance of each other */
    61. 

  61. Datum RASTER_dwithin(PG_FUNCTION_ARGS);
    62. 

  62. 

  63. /* determine if the two rasters are fully within the specified distance of each other */
    64. 

  64. Datum RASTER_dfullywithin(PG_FUNCTION_ARGS);
    65. 

  65. 

  66. /* determine if two rasters are aligned */
    67. 

  67. Datum RASTER_sameAlignment(PG_FUNCTION_ARGS);
    68. 

  68. Datum RASTER_notSameAlignmentReason(PG_FUNCTION_ARGS);
    69. 

  69. 

  70. /**
    71. 

  71.  * See if two rasters intersect
    72. 

  72.  */
    73. 

  73. PG_FUNCTION_INFO_V1(RASTER_intersects);
    74. 

  74. Datum RASTER_intersects(PG_FUNCTION_ARGS)
    75. 

  75. {
    76. 

  76. 	const int set_count = 2;
    77. 

  77. 	rt_pgraster *pgrast[2];
    78. 

  78. 	int pgrastpos[2] = {-1, -1};
    79. 

  79. 	rt_raster rast[2] = {NULL};
    80. 

  80. 	uint32_t bandindex[2] = {0};
    81. 

  81. 	uint32_t hasbandindex[2] = {0};
    82. 

  82. 

  83. 	uint32_t i;
    84. 

  84. 	uint32_t j;
    85. 

  85. 	uint32_t k;
    86. 

  86. 	uint32_t numBands;
    87. 

  87. 	int rtn;
    88. 

  88. 	int result;
    89. 

  89. 

  90. 	for (i = 0, j = 0; i < set_count; i++) {
    91. 

  91. 		/* pgrast is null, return null */
    92. 

  92. 		if (PG_ARGISNULL(j)) {
    93. 

  93. 			for (k = 0; k < i; k++) {
    94. 

  94. 				rt_raster_destroy(rast[k]);
    95. 

  95. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    96. 

  96. 			}
    97. 

  97. 			PG_RETURN_NULL();
    98. 

  98. 		}
    99. 

  99. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    100. 

  100. 		pgrastpos[i] = j;
    101. 

  101. 		j++;
    102. 

  102. 

  103. 		/* raster */
    104. 

  104. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    105. 

  105. 		if (!rast[i]) {
    106. 

  106. 			for (k = 0; k <= i; k++) {
    107. 

  107. 				if (k < i)
    108. 

  108. 					rt_raster_destroy(rast[k]);
    109. 

  109. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    110. 

  110. 			}
    111. 

  111. 			elog(ERROR, "RASTER_intersects: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    112. 

  112. 			PG_RETURN_NULL();
    113. 

  113. 		}
    114. 

  114. 

  115. 		/* numbands */
    116. 

  116. 		numBands = rt_raster_get_num_bands(rast[i]);
    117. 

  117. 		if (numBands < 1) {
    118. 

  118. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    119. 

  119. 			if (i > 0) i++;
    120. 

  120. 			for (k = 0; k < i; k++) {
    121. 

  121. 				rt_raster_destroy(rast[k]);
    122. 

  122. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    123. 

  123. 			}
    124. 

  124. 			PG_RETURN_NULL();
    125. 

  125. 		}
    126. 

  126. 

  127. 		/* band index */
    128. 

  128. 		if (!PG_ARGISNULL(j)) {
    129. 

  129. 			bandindex[i] = PG_GETARG_INT32(j);
    130. 

  130. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    131. 

  131. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    132. 

  132. 				if (i > 0) i++;
    133. 

  133. 				for (k = 0; k < i; k++) {
    134. 

  134. 					rt_raster_destroy(rast[k]);
    135. 

  135. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    136. 

  136. 				}
    137. 

  137. 				PG_RETURN_NULL();
    138. 

  138. 			}
    139. 

  139. 			hasbandindex[i] = 1;
    140. 

  140. 		}
    141. 

  141. 		else
    142. 

  142. 			hasbandindex[i] = 0;
    143. 

  143. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    144. 

  144. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    145. 

  145. 		j++;
    146. 

  146. 	}
    147. 

  147. 

  148. 	/* hasbandindex must be balanced */
    149. 

  149. 	if (
    150. 

  150. 		(hasbandindex[0] && !hasbandindex[1]) ||
    151. 

  151. 		(!hasbandindex[0] && hasbandindex[1])
    152. 

  152. 	) {
    153. 

  153. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    154. 

  154. 		for (k = 0; k < set_count; k++) {
    155. 

  155. 			rt_raster_destroy(rast[k]);
    156. 

  156. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    157. 

  157. 		}
    158. 

  158. 		PG_RETURN_NULL();
    159. 

  159. 	}
    160. 

  160. 

  161. 	/* SRID must match */
    162. 

  162. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    163. 

  163. 		for (k = 0; k < set_count; k++) {
    164. 

  164. 			rt_raster_destroy(rast[k]);
    165. 

  165. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    166. 

  166. 		}
    167. 

  167. 		elog(ERROR, "The two rasters provided have different SRIDs");
    168. 

  168. 		PG_RETURN_NULL();
    169. 

  169. 	}
    170. 

  170. 

  171. 	rtn = rt_raster_intersects(
    172. 

  172. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    173. 

  173. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    174. 

  174. 		&result
    175. 

  175. 	);
    176. 

  176. 	for (k = 0; k < set_count; k++) {
    177. 

  177. 		rt_raster_destroy(rast[k]);
    178. 

  178. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    179. 

  179. 	}
    180. 

  180. 

  181. 	if (rtn != ES_NONE) {
    182. 

  182. 		elog(ERROR, "RASTER_intersects: Could not test for intersection on the two rasters");
    183. 

  183. 		PG_RETURN_NULL();
    184. 

  184. 	}
    185. 

  185. 

  186. 	PG_RETURN_BOOL(result);
    187. 

  187. }
    188. 

  188. 

  189. /**
    190. 

  190.  * See if two rasters overlap
    191. 

  191.  */
    192. 

  192. PG_FUNCTION_INFO_V1(RASTER_overlaps);
    193. 

  193. Datum RASTER_overlaps(PG_FUNCTION_ARGS)
    194. 

  194. {
    195. 

  195. 	const int set_count = 2;
    196. 

  196. 	rt_pgraster *pgrast[2];
    197. 

  197. 	int pgrastpos[2] = {-1, -1};
    198. 

  198. 	rt_raster rast[2] = {NULL};
    199. 

  199. 	uint32_t bandindex[2] = {0};
    200. 

  200. 	uint32_t hasbandindex[2] = {0};
    201. 

  201. 

  202. 	uint32_t i;
    203. 

  203. 	uint32_t j;
    204. 

  204. 	uint32_t k;
    205. 

  205. 	uint32_t numBands;
    206. 

  206. 	int rtn;
    207. 

  207. 	int result;
    208. 

  208. 

  209. 	for (i = 0, j = 0; i < set_count; i++) {
    210. 

  210. 		/* pgrast is null, return null */
    211. 

  211. 		if (PG_ARGISNULL(j)) {
    212. 

  212. 			for (k = 0; k < i; k++) {
    213. 

  213. 				rt_raster_destroy(rast[k]);
    214. 

  214. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    215. 

  215. 			}
    216. 

  216. 			PG_RETURN_NULL();
    217. 

  217. 		}
    218. 

  218. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    219. 

  219. 		pgrastpos[i] = j;
    220. 

  220. 		j++;
    221. 

  221. 

  222. 		/* raster */
    223. 

  223. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    224. 

  224. 		if (!rast[i]) {
    225. 

  225. 			for (k = 0; k <= i; k++) {
    226. 

  226. 				if (k < i)
    227. 

  227. 					rt_raster_destroy(rast[k]);
    228. 

  228. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    229. 

  229. 			}
    230. 

  230. 			elog(ERROR, "RASTER_overlaps: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    231. 

  231. 			PG_RETURN_NULL();
    232. 

  232. 		}
    233. 

  233. 

  234. 		/* numbands */
    235. 

  235. 		numBands = rt_raster_get_num_bands(rast[i]);
    236. 

  236. 		if (numBands < 1) {
    237. 

  237. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    238. 

  238. 			if (i > 0) i++;
    239. 

  239. 			for (k = 0; k < i; k++) {
    240. 

  240. 				rt_raster_destroy(rast[k]);
    241. 

  241. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    242. 

  242. 			}
    243. 

  243. 			PG_RETURN_NULL();
    244. 

  244. 		}
    245. 

  245. 

  246. 		/* band index */
    247. 

  247. 		if (!PG_ARGISNULL(j)) {
    248. 

  248. 			bandindex[i] = PG_GETARG_INT32(j);
    249. 

  249. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    250. 

  250. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    251. 

  251. 				if (i > 0) i++;
    252. 

  252. 				for (k = 0; k < i; k++) {
    253. 

  253. 					rt_raster_destroy(rast[k]);
    254. 

  254. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    255. 

  255. 				}
    256. 

  256. 				PG_RETURN_NULL();
    257. 

  257. 			}
    258. 

  258. 			hasbandindex[i] = 1;
    259. 

  259. 		}
    260. 

  260. 		else
    261. 

  261. 			hasbandindex[i] = 0;
    262. 

  262. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    263. 

  263. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    264. 

  264. 		j++;
    265. 

  265. 	}
    266. 

  266. 

  267. 	/* hasbandindex must be balanced */
    268. 

  268. 	if (
    269. 

  269. 		(hasbandindex[0] && !hasbandindex[1]) ||
    270. 

  270. 		(!hasbandindex[0] && hasbandindex[1])
    271. 

  271. 	) {
    272. 

  272. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    273. 

  273. 		for (k = 0; k < set_count; k++) {
    274. 

  274. 			rt_raster_destroy(rast[k]);
    275. 

  275. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    276. 

  276. 		}
    277. 

  277. 		PG_RETURN_NULL();
    278. 

  278. 	}
    279. 

  279. 

  280. 	/* SRID must match */
    281. 

  281. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    282. 

  282. 		for (k = 0; k < set_count; k++) {
    283. 

  283. 			rt_raster_destroy(rast[k]);
    284. 

  284. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    285. 

  285. 		}
    286. 

  286. 		elog(ERROR, "The two rasters provided have different SRIDs");
    287. 

  287. 		PG_RETURN_NULL();
    288. 

  288. 	}
    289. 

  289. 

  290. 	rtn = rt_raster_overlaps(
    291. 

  291. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    292. 

  292. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    293. 

  293. 		&result
    294. 

  294. 	);
    295. 

  295. 	for (k = 0; k < set_count; k++) {
    296. 

  296. 		rt_raster_destroy(rast[k]);
    297. 

  297. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    298. 

  298. 	}
    299. 

  299. 

  300. 	if (rtn != ES_NONE) {
    301. 

  301. 		elog(ERROR, "RASTER_overlaps: Could not test for overlap on the two rasters");
    302. 

  302. 		PG_RETURN_NULL();
    303. 

  303. 	}
    304. 

  304. 

  305. 	PG_RETURN_BOOL(result);
    306. 

  306. }
    307. 

  307. 

  308. /**
    309. 

  309.  * See if two rasters touch
    310. 

  310.  */
    311. 

  311. PG_FUNCTION_INFO_V1(RASTER_touches);
    312. 

  312. Datum RASTER_touches(PG_FUNCTION_ARGS)
    313. 

  313. {
    314. 

  314. 	const int set_count = 2;
    315. 

  315. 	rt_pgraster *pgrast[2];
    316. 

  316. 	int pgrastpos[2] = {-1, -1};
    317. 

  317. 	rt_raster rast[2] = {NULL};
    318. 

  318. 	uint32_t bandindex[2] = {0};
    319. 

  319. 	uint32_t hasbandindex[2] = {0};
    320. 

  320. 

  321. 	uint32_t i;
    322. 

  322. 	uint32_t j;
    323. 

  323. 	uint32_t k;
    324. 

  324. 	uint32_t numBands;
    325. 

  325. 	int rtn;
    326. 

  326. 	int result;
    327. 

  327. 

  328. 	for (i = 0, j = 0; i < set_count; i++) {
    329. 

  329. 		/* pgrast is null, return null */
    330. 

  330. 		if (PG_ARGISNULL(j)) {
    331. 

  331. 			for (k = 0; k < i; k++) {
    332. 

  332. 				rt_raster_destroy(rast[k]);
    333. 

  333. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    334. 

  334. 			}
    335. 

  335. 			PG_RETURN_NULL();
    336. 

  336. 		}
    337. 

  337. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    338. 

  338. 		pgrastpos[i] = j;
    339. 

  339. 		j++;
    340. 

  340. 

  341. 		/* raster */
    342. 

  342. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    343. 

  343. 		if (!rast[i]) {
    344. 

  344. 			for (k = 0; k <= i; k++) {
    345. 

  345. 				if (k < i)
    346. 

  346. 					rt_raster_destroy(rast[k]);
    347. 

  347. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    348. 

  348. 			}
    349. 

  349. 			elog(ERROR, "RASTER_touches: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    350. 

  350. 			PG_RETURN_NULL();
    351. 

  351. 		}
    352. 

  352. 

  353. 		/* numbands */
    354. 

  354. 		numBands = rt_raster_get_num_bands(rast[i]);
    355. 

  355. 		if (numBands < 1) {
    356. 

  356. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    357. 

  357. 			if (i > 0) i++;
    358. 

  358. 			for (k = 0; k < i; k++) {
    359. 

  359. 				rt_raster_destroy(rast[k]);
    360. 

  360. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    361. 

  361. 			}
    362. 

  362. 			PG_RETURN_NULL();
    363. 

  363. 		}
    364. 

  364. 

  365. 		/* band index */
    366. 

  366. 		if (!PG_ARGISNULL(j)) {
    367. 

  367. 			bandindex[i] = PG_GETARG_INT32(j);
    368. 

  368. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    369. 

  369. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    370. 

  370. 				if (i > 0) i++;
    371. 

  371. 				for (k = 0; k < i; k++) {
    372. 

  372. 					rt_raster_destroy(rast[k]);
    373. 

  373. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    374. 

  374. 				}
    375. 

  375. 				PG_RETURN_NULL();
    376. 

  376. 			}
    377. 

  377. 			hasbandindex[i] = 1;
    378. 

  378. 		}
    379. 

  379. 		else
    380. 

  380. 			hasbandindex[i] = 0;
    381. 

  381. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    382. 

  382. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    383. 

  383. 		j++;
    384. 

  384. 	}
    385. 

  385. 

  386. 	/* hasbandindex must be balanced */
    387. 

  387. 	if (
    388. 

  388. 		(hasbandindex[0] && !hasbandindex[1]) ||
    389. 

  389. 		(!hasbandindex[0] && hasbandindex[1])
    390. 

  390. 	) {
    391. 

  391. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    392. 

  392. 		for (k = 0; k < set_count; k++) {
    393. 

  393. 			rt_raster_destroy(rast[k]);
    394. 

  394. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    395. 

  395. 		}
    396. 

  396. 		PG_RETURN_NULL();
    397. 

  397. 	}
    398. 

  398. 

  399. 	/* SRID must match */
    400. 

  400. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    401. 

  401. 		for (k = 0; k < set_count; k++) {
    402. 

  402. 			rt_raster_destroy(rast[k]);
    403. 

  403. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    404. 

  404. 		}
    405. 

  405. 		elog(ERROR, "The two rasters provided have different SRIDs");
    406. 

  406. 		PG_RETURN_NULL();
    407. 

  407. 	}
    408. 

  408. 

  409. 	rtn = rt_raster_touches(
    410. 

  410. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    411. 

  411. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    412. 

  412. 		&result
    413. 

  413. 	);
    414. 

  414. 	for (k = 0; k < set_count; k++) {
    415. 

  415. 		rt_raster_destroy(rast[k]);
    416. 

  416. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    417. 

  417. 	}
    418. 

  418. 

  419. 	if (rtn != ES_NONE) {
    420. 

  420. 		elog(ERROR, "RASTER_touches: Could not test for touch on the two rasters");
    421. 

  421. 		PG_RETURN_NULL();
    422. 

  422. 	}
    423. 

  423. 

  424. 	PG_RETURN_BOOL(result);
    425. 

  425. }
    426. 

  426. 

  427. /**
    428. 

  428.  * See if the first raster contains the second raster
    429. 

  429.  */
    430. 

  430. PG_FUNCTION_INFO_V1(RASTER_contains);
    431. 

  431. Datum RASTER_contains(PG_FUNCTION_ARGS)
    432. 

  432. {
    433. 

  433. 	const int set_count = 2;
    434. 

  434. 	rt_pgraster *pgrast[2];
    435. 

  435. 	int pgrastpos[2] = {-1, -1};
    436. 

  436. 	rt_raster rast[2] = {NULL};
    437. 

  437. 	uint32_t bandindex[2] = {0};
    438. 

  438. 	uint32_t hasbandindex[2] = {0};
    439. 

  439. 

  440. 	uint32_t i;
    441. 

  441. 	uint32_t j;
    442. 

  442. 	uint32_t k;
    443. 

  443. 	uint32_t numBands;
    444. 

  444. 	int rtn;
    445. 

  445. 	int result;
    446. 

  446. 

  447. 	for (i = 0, j = 0; i < set_count; i++) {
    448. 

  448. 		/* pgrast is null, return null */
    449. 

  449. 		if (PG_ARGISNULL(j)) {
    450. 

  450. 			for (k = 0; k < i; k++) {
    451. 

  451. 				rt_raster_destroy(rast[k]);
    452. 

  452. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    453. 

  453. 			}
    454. 

  454. 			PG_RETURN_NULL();
    455. 

  455. 		}
    456. 

  456. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    457. 

  457. 		pgrastpos[i] = j;
    458. 

  458. 		j++;
    459. 

  459. 

  460. 		/* raster */
    461. 

  461. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    462. 

  462. 		if (!rast[i]) {
    463. 

  463. 			for (k = 0; k <= i; k++) {
    464. 

  464. 				if (k < i)
    465. 

  465. 					rt_raster_destroy(rast[k]);
    466. 

  466. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    467. 

  467. 			}
    468. 

  468. 			elog(ERROR, "RASTER_contains: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    469. 

  469. 			PG_RETURN_NULL();
    470. 

  470. 		}
    471. 

  471. 

  472. 		/* numbands */
    473. 

  473. 		numBands = rt_raster_get_num_bands(rast[i]);
    474. 

  474. 		if (numBands < 1) {
    475. 

  475. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    476. 

  476. 			if (i > 0) i++;
    477. 

  477. 			for (k = 0; k < i; k++) {
    478. 

  478. 				rt_raster_destroy(rast[k]);
    479. 

  479. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    480. 

  480. 			}
    481. 

  481. 			PG_RETURN_NULL();
    482. 

  482. 		}
    483. 

  483. 

  484. 		/* band index */
    485. 

  485. 		if (!PG_ARGISNULL(j)) {
    486. 

  486. 			bandindex[i] = PG_GETARG_INT32(j);
    487. 

  487. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    488. 

  488. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    489. 

  489. 				if (i > 0) i++;
    490. 

  490. 				for (k = 0; k < i; k++) {
    491. 

  491. 					rt_raster_destroy(rast[k]);
    492. 

  492. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    493. 

  493. 				}
    494. 

  494. 				PG_RETURN_NULL();
    495. 

  495. 			}
    496. 

  496. 			hasbandindex[i] = 1;
    497. 

  497. 		}
    498. 

  498. 		else
    499. 

  499. 			hasbandindex[i] = 0;
    500. 

  500. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    501. 

  501. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    502. 

  502. 		j++;
    503. 

  503. 	}
    504. 

  504. 

  505. 	/* hasbandindex must be balanced */
    506. 

  506. 	if (
    507. 

  507. 		(hasbandindex[0] && !hasbandindex[1]) ||
    508. 

  508. 		(!hasbandindex[0] && hasbandindex[1])
    509. 

  509. 	) {
    510. 

  510. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    511. 

  511. 		for (k = 0; k < set_count; k++) {
    512. 

  512. 			rt_raster_destroy(rast[k]);
    513. 

  513. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    514. 

  514. 		}
    515. 

  515. 		PG_RETURN_NULL();
    516. 

  516. 	}
    517. 

  517. 

  518. 	/* SRID must match */
    519. 

  519. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    520. 

  520. 		for (k = 0; k < set_count; k++) {
    521. 

  521. 			rt_raster_destroy(rast[k]);
    522. 

  522. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    523. 

  523. 		}
    524. 

  524. 		elog(ERROR, "The two rasters provided have different SRIDs");
    525. 

  525. 		PG_RETURN_NULL();
    526. 

  526. 	}
    527. 

  527. 

  528. 	rtn = rt_raster_contains(
    529. 

  529. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    530. 

  530. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    531. 

  531. 		&result
    532. 

  532. 	);
    533. 

  533. 	for (k = 0; k < set_count; k++) {
    534. 

  534. 		rt_raster_destroy(rast[k]);
    535. 

  535. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    536. 

  536. 	}
    537. 

  537. 

  538. 	if (rtn != ES_NONE) {
    539. 

  539. 		elog(ERROR, "RASTER_contains: Could not test that the first raster contains the second raster");
    540. 

  540. 		PG_RETURN_NULL();
    541. 

  541. 	}
    542. 

  542. 

  543. 	PG_RETURN_BOOL(result);
    544. 

  544. }
    545. 

  545. 

  546. /**
    547. 

  547.  * See if the first raster contains properly the second raster
    548. 

  548.  */
    549. 

  549. PG_FUNCTION_INFO_V1(RASTER_containsProperly);
    550. 

  550. Datum RASTER_containsProperly(PG_FUNCTION_ARGS)
    551. 

  551. {
    552. 

  552. 	const int set_count = 2;
    553. 

  553. 	rt_pgraster *pgrast[2];
    554. 

  554. 	int pgrastpos[2] = {-1, -1};
    555. 

  555. 	rt_raster rast[2] = {NULL};
    556. 

  556. 	uint32_t bandindex[2] = {0};
    557. 

  557. 	uint32_t hasbandindex[2] = {0};
    558. 

  558. 

  559. 	uint32_t i;
    560. 

  560. 	uint32_t j;
    561. 

  561. 	uint32_t k;
    562. 

  562. 	uint32_t numBands;
    563. 

  563. 	int rtn;
    564. 

  564. 	int result;
    565. 

  565. 

  566. 	for (i = 0, j = 0; i < set_count; i++) {
    567. 

  567. 		/* pgrast is null, return null */
    568. 

  568. 		if (PG_ARGISNULL(j)) {
    569. 

  569. 			for (k = 0; k < i; k++) {
    570. 

  570. 				rt_raster_destroy(rast[k]);
    571. 

  571. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    572. 

  572. 			}
    573. 

  573. 			PG_RETURN_NULL();
    574. 

  574. 		}
    575. 

  575. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    576. 

  576. 		pgrastpos[i] = j;
    577. 

  577. 		j++;
    578. 

  578. 

  579. 		/* raster */
    580. 

  580. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    581. 

  581. 		if (!rast[i]) {
    582. 

  582. 			for (k = 0; k <= i; k++) {
    583. 

  583. 				if (k < i)
    584. 

  584. 					rt_raster_destroy(rast[k]);
    585. 

  585. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    586. 

  586. 			}
    587. 

  587. 			elog(ERROR, "RASTER_containsProperly: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    588. 

  588. 			PG_RETURN_NULL();
    589. 

  589. 		}
    590. 

  590. 

  591. 		/* numbands */
    592. 

  592. 		numBands = rt_raster_get_num_bands(rast[i]);
    593. 

  593. 		if (numBands < 1) {
    594. 

  594. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    595. 

  595. 			if (i > 0) i++;
    596. 

  596. 			for (k = 0; k < i; k++) {
    597. 

  597. 				rt_raster_destroy(rast[k]);
    598. 

  598. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    599. 

  599. 			}
    600. 

  600. 			PG_RETURN_NULL();
    601. 

  601. 		}
    602. 

  602. 

  603. 		/* band index */
    604. 

  604. 		if (!PG_ARGISNULL(j)) {
    605. 

  605. 			bandindex[i] = PG_GETARG_INT32(j);
    606. 

  606. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    607. 

  607. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    608. 

  608. 				if (i > 0) i++;
    609. 

  609. 				for (k = 0; k < i; k++) {
    610. 

  610. 					rt_raster_destroy(rast[k]);
    611. 

  611. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    612. 

  612. 				}
    613. 

  613. 				PG_RETURN_NULL();
    614. 

  614. 			}
    615. 

  615. 			hasbandindex[i] = 1;
    616. 

  616. 		}
    617. 

  617. 		else
    618. 

  618. 			hasbandindex[i] = 0;
    619. 

  619. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    620. 

  620. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    621. 

  621. 		j++;
    622. 

  622. 	}
    623. 

  623. 

  624. 	/* hasbandindex must be balanced */
    625. 

  625. 	if (
    626. 

  626. 		(hasbandindex[0] && !hasbandindex[1]) ||
    627. 

  627. 		(!hasbandindex[0] && hasbandindex[1])
    628. 

  628. 	) {
    629. 

  629. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    630. 

  630. 		for (k = 0; k < set_count; k++) {
    631. 

  631. 			rt_raster_destroy(rast[k]);
    632. 

  632. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    633. 

  633. 		}
    634. 

  634. 		PG_RETURN_NULL();
    635. 

  635. 	}
    636. 

  636. 

  637. 	/* SRID must match */
    638. 

  638. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    639. 

  639. 		for (k = 0; k < set_count; k++) {
    640. 

  640. 			rt_raster_destroy(rast[k]);
    641. 

  641. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    642. 

  642. 		}
    643. 

  643. 		elog(ERROR, "The two rasters provided have different SRIDs");
    644. 

  644. 		PG_RETURN_NULL();
    645. 

  645. 	}
    646. 

  646. 

  647. 	rtn = rt_raster_contains_properly(
    648. 

  648. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    649. 

  649. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    650. 

  650. 		&result
    651. 

  651. 	);
    652. 

  652. 	for (k = 0; k < set_count; k++) {
    653. 

  653. 		rt_raster_destroy(rast[k]);
    654. 

  654. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    655. 

  655. 	}
    656. 

  656. 

  657. 	if (rtn != ES_NONE) {
    658. 

  658. 		elog(ERROR, "RASTER_containsProperly: Could not test that the first raster contains properly the second raster");
    659. 

  659. 		PG_RETURN_NULL();
    660. 

  660. 	}
    661. 

  661. 

  662. 	PG_RETURN_BOOL(result);
    663. 

  663. }
    664. 

  664. 

  665. /**
    666. 

  666.  * See if the first raster covers the second raster
    667. 

  667.  */
    668. 

  668. PG_FUNCTION_INFO_V1(RASTER_covers);
    669. 

  669. Datum RASTER_covers(PG_FUNCTION_ARGS)
    670. 

  670. {
    671. 

  671. 	const int set_count = 2;
    672. 

  672. 	rt_pgraster *pgrast[2];
    673. 

  673. 	int pgrastpos[2] = {-1, -1};
    674. 

  674. 	rt_raster rast[2] = {NULL};
    675. 

  675. 	uint32_t bandindex[2] = {0};
    676. 

  676. 	uint32_t hasbandindex[2] = {0};
    677. 

  677. 

  678. 	uint32_t i;
    679. 

  679. 	uint32_t j;
    680. 

  680. 	uint32_t k;
    681. 

  681. 	uint32_t numBands;
    682. 

  682. 	int rtn;
    683. 

  683. 	int result;
    684. 

  684. 

  685. 	for (i = 0, j = 0; i < set_count; i++) {
    686. 

  686. 		/* pgrast is null, return null */
    687. 

  687. 		if (PG_ARGISNULL(j)) {
    688. 

  688. 			for (k = 0; k < i; k++) {
    689. 

  689. 				rt_raster_destroy(rast[k]);
    690. 

  690. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    691. 

  691. 			}
    692. 

  692. 			PG_RETURN_NULL();
    693. 

  693. 		}
    694. 

  694. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    695. 

  695. 		pgrastpos[i] = j;
    696. 

  696. 		j++;
    697. 

  697. 

  698. 		/* raster */
    699. 

  699. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    700. 

  700. 		if (!rast[i]) {
    701. 

  701. 			for (k = 0; k <= i; k++) {
    702. 

  702. 				if (k < i)
    703. 

  703. 					rt_raster_destroy(rast[k]);
    704. 

  704. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    705. 

  705. 			}
    706. 

  706. 			elog(ERROR, "RASTER_covers: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    707. 

  707. 			PG_RETURN_NULL();
    708. 

  708. 		}
    709. 

  709. 

  710. 		/* numbands */
    711. 

  711. 		numBands = rt_raster_get_num_bands(rast[i]);
    712. 

  712. 		if (numBands < 1) {
    713. 

  713. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    714. 

  714. 			if (i > 0) i++;
    715. 

  715. 			for (k = 0; k < i; k++) {
    716. 

  716. 				rt_raster_destroy(rast[k]);
    717. 

  717. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    718. 

  718. 			}
    719. 

  719. 			PG_RETURN_NULL();
    720. 

  720. 		}
    721. 

  721. 

  722. 		/* band index */
    723. 

  723. 		if (!PG_ARGISNULL(j)) {
    724. 

  724. 			bandindex[i] = PG_GETARG_INT32(j);
    725. 

  725. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    726. 

  726. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    727. 

  727. 				if (i > 0) i++;
    728. 

  728. 				for (k = 0; k < i; k++) {
    729. 

  729. 					rt_raster_destroy(rast[k]);
    730. 

  730. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    731. 

  731. 				}
    732. 

  732. 				PG_RETURN_NULL();
    733. 

  733. 			}
    734. 

  734. 			hasbandindex[i] = 1;
    735. 

  735. 		}
    736. 

  736. 		else
    737. 

  737. 			hasbandindex[i] = 0;
    738. 

  738. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    739. 

  739. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    740. 

  740. 		j++;
    741. 

  741. 	}
    742. 

  742. 

  743. 	/* hasbandindex must be balanced */
    744. 

  744. 	if (
    745. 

  745. 		(hasbandindex[0] && !hasbandindex[1]) ||
    746. 

  746. 		(!hasbandindex[0] && hasbandindex[1])
    747. 

  747. 	) {
    748. 

  748. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    749. 

  749. 		for (k = 0; k < set_count; k++) {
    750. 

  750. 			rt_raster_destroy(rast[k]);
    751. 

  751. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    752. 

  752. 		}
    753. 

  753. 		PG_RETURN_NULL();
    754. 

  754. 	}
    755. 

  755. 

  756. 	/* SRID must match */
    757. 

  757. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    758. 

  758. 		for (k = 0; k < set_count; k++) {
    759. 

  759. 			rt_raster_destroy(rast[k]);
    760. 

  760. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    761. 

  761. 		}
    762. 

  762. 		elog(ERROR, "The two rasters provided have different SRIDs");
    763. 

  763. 		PG_RETURN_NULL();
    764. 

  764. 	}
    765. 

  765. 

  766. 	rtn = rt_raster_covers(
    767. 

  767. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    768. 

  768. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    769. 

  769. 		&result
    770. 

  770. 	);
    771. 

  771. 	for (k = 0; k < set_count; k++) {
    772. 

  772. 		rt_raster_destroy(rast[k]);
    773. 

  773. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    774. 

  774. 	}
    775. 

  775. 

  776. 	if (rtn != ES_NONE) {
    777. 

  777. 		elog(ERROR, "RASTER_covers: Could not test that the first raster covers the second raster");
    778. 

  778. 		PG_RETURN_NULL();
    779. 

  779. 	}
    780. 

  780. 

  781. 	PG_RETURN_BOOL(result);
    782. 

  782. }
    783. 

  783. 

  784. /**
    785. 

  785.  * See if the first raster is covered by the second raster
    786. 

  786.  */
    787. 

  787. PG_FUNCTION_INFO_V1(RASTER_coveredby);
    788. 

  788. Datum RASTER_coveredby(PG_FUNCTION_ARGS)
    789. 

  789. {
    790. 

  790. 	const int set_count = 2;
    791. 

  791. 	rt_pgraster *pgrast[2];
    792. 

  792. 	int pgrastpos[2] = {-1, -1};
    793. 

  793. 	rt_raster rast[2] = {NULL};
    794. 

  794. 	uint32_t bandindex[2] = {0};
    795. 

  795. 	uint32_t hasbandindex[2] = {0};
    796. 

  796. 

  797. 	uint32_t i;
    798. 

  798. 	uint32_t j;
    799. 

  799. 	uint32_t k;
    800. 

  800. 	uint32_t numBands;
    801. 

  801. 	int rtn;
    802. 

  802. 	int result;
    803. 

  803. 

  804. 	for (i = 0, j = 0; i < set_count; i++) {
    805. 

  805. 		/* pgrast is null, return null */
    806. 

  806. 		if (PG_ARGISNULL(j)) {
    807. 

  807. 			for (k = 0; k < i; k++) {
    808. 

  808. 				rt_raster_destroy(rast[k]);
    809. 

  809. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    810. 

  810. 			}
    811. 

  811. 			PG_RETURN_NULL();
    812. 

  812. 		}
    813. 

  813. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    814. 

  814. 		pgrastpos[i] = j;
    815. 

  815. 		j++;
    816. 

  816. 

  817. 		/* raster */
    818. 

  818. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    819. 

  819. 		if (!rast[i]) {
    820. 

  820. 			for (k = 0; k <= i; k++) {
    821. 

  821. 				if (k < i)
    822. 

  822. 					rt_raster_destroy(rast[k]);
    823. 

  823. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    824. 

  824. 			}
    825. 

  825. 			elog(ERROR, "RASTER_coveredby: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    826. 

  826. 			PG_RETURN_NULL();
    827. 

  827. 		}
    828. 

  828. 

  829. 		/* numbands */
    830. 

  830. 		numBands = rt_raster_get_num_bands(rast[i]);
    831. 

  831. 		if (numBands < 1) {
    832. 

  832. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    833. 

  833. 			if (i > 0) i++;
    834. 

  834. 			for (k = 0; k < i; k++) {
    835. 

  835. 				rt_raster_destroy(rast[k]);
    836. 

  836. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    837. 

  837. 			}
    838. 

  838. 			PG_RETURN_NULL();
    839. 

  839. 		}
    840. 

  840. 

  841. 		/* band index */
    842. 

  842. 		if (!PG_ARGISNULL(j)) {
    843. 

  843. 			bandindex[i] = PG_GETARG_INT32(j);
    844. 

  844. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    845. 

  845. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    846. 

  846. 				if (i > 0) i++;
    847. 

  847. 				for (k = 0; k < i; k++) {
    848. 

  848. 					rt_raster_destroy(rast[k]);
    849. 

  849. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    850. 

  850. 				}
    851. 

  851. 				PG_RETURN_NULL();
    852. 

  852. 			}
    853. 

  853. 			hasbandindex[i] = 1;
    854. 

  854. 		}
    855. 

  855. 		else
    856. 

  856. 			hasbandindex[i] = 0;
    857. 

  857. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    858. 

  858. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    859. 

  859. 		j++;
    860. 

  860. 	}
    861. 

  861. 

  862. 	/* hasbandindex must be balanced */
    863. 

  863. 	if (
    864. 

  864. 		(hasbandindex[0] && !hasbandindex[1]) ||
    865. 

  865. 		(!hasbandindex[0] && hasbandindex[1])
    866. 

  866. 	) {
    867. 

  867. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    868. 

  868. 		for (k = 0; k < set_count; k++) {
    869. 

  869. 			rt_raster_destroy(rast[k]);
    870. 

  870. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    871. 

  871. 		}
    872. 

  872. 		PG_RETURN_NULL();
    873. 

  873. 	}
    874. 

  874. 

  875. 	/* SRID must match */
    876. 

  876. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    877. 

  877. 		for (k = 0; k < set_count; k++) {
    878. 

  878. 			rt_raster_destroy(rast[k]);
    879. 

  879. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    880. 

  880. 		}
    881. 

  881. 		elog(ERROR, "The two rasters provided have different SRIDs");
    882. 

  882. 		PG_RETURN_NULL();
    883. 

  883. 	}
    884. 

  884. 

  885. 	rtn = rt_raster_coveredby(
    886. 

  886. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    887. 

  887. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    888. 

  888. 		&result
    889. 

  889. 	);
    890. 

  890. 	for (k = 0; k < set_count; k++) {
    891. 

  891. 		rt_raster_destroy(rast[k]);
    892. 

  892. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    893. 

  893. 	}
    894. 

  894. 

  895. 	if (rtn != ES_NONE) {
    896. 

  896. 		elog(ERROR, "RASTER_coveredby: Could not test that the first raster is covered by the second raster");
    897. 

  897. 		PG_RETURN_NULL();
    898. 

  898. 	}
    899. 

  899. 

  900. 	PG_RETURN_BOOL(result);
    901. 

  901. }
    902. 

  902. 

  903. /**
    904. 

  904.  * See if the two rasters are within the specified distance of each other
    905. 

  905.  */
    906. 

  906. PG_FUNCTION_INFO_V1(RASTER_dwithin);
    907. 

  907. Datum RASTER_dwithin(PG_FUNCTION_ARGS)
    908. 

  908. {
    909. 

  909. 	const int set_count = 2;
    910. 

  910. 	rt_pgraster *pgrast[2];
    911. 

  911. 	int pgrastpos[2] = {-1, -1};
    912. 

  912. 	rt_raster rast[2] = {NULL};
    913. 

  913. 	uint32_t bandindex[2] = {0};
    914. 

  914. 	uint32_t hasbandindex[2] = {0};
    915. 

  915. 	double distance = 0;
    916. 

  916. 

  917. 	uint32_t i;
    918. 

  918. 	uint32_t j;
    919. 

  919. 	uint32_t k;
    920. 

  920. 	uint32_t numBands;
    921. 

  921. 	int rtn;
    922. 

  922. 	int result;
    923. 

  923. 

  924. 	for (i = 0, j = 0; i < set_count; i++) {
    925. 

  925. 		/* pgrast is null, return null */
    926. 

  926. 		if (PG_ARGISNULL(j)) {
    927. 

  927. 			for (k = 0; k < i; k++) {
    928. 

  928. 				rt_raster_destroy(rast[k]);
    929. 

  929. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    930. 

  930. 			}
    931. 

  931. 			PG_RETURN_NULL();
    932. 

  932. 		}
    933. 

  933. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    934. 

  934. 		pgrastpos[i] = j;
    935. 

  935. 		j++;
    936. 

  936. 

  937. 		/* raster */
    938. 

  938. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    939. 

  939. 		if (!rast[i]) {
    940. 

  940. 			for (k = 0; k <= i; k++) {
    941. 

  941. 				if (k < i)
    942. 

  942. 					rt_raster_destroy(rast[k]);
    943. 

  943. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    944. 

  944. 			}
    945. 

  945. 			elog(ERROR, "RASTER_dwithin: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    946. 

  946. 			PG_RETURN_NULL();
    947. 

  947. 		}
    948. 

  948. 

  949. 		/* numbands */
    950. 

  950. 		numBands = rt_raster_get_num_bands(rast[i]);
    951. 

  951. 		if (numBands < 1) {
    952. 

  952. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    953. 

  953. 			if (i > 0) i++;
    954. 

  954. 			for (k = 0; k < i; k++) {
    955. 

  955. 				rt_raster_destroy(rast[k]);
    956. 

  956. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    957. 

  957. 			}
    958. 

  958. 			PG_RETURN_NULL();
    959. 

  959. 		}
    960. 

  960. 

  961. 		/* band index */
    962. 

  962. 		if (!PG_ARGISNULL(j)) {
    963. 

  963. 			bandindex[i] = PG_GETARG_INT32(j);
    964. 

  964. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    965. 

  965. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    966. 

  966. 				if (i > 0) i++;
    967. 

  967. 				for (k = 0; k < i; k++) {
    968. 

  968. 					rt_raster_destroy(rast[k]);
    969. 

  969. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    970. 

  970. 				}
    971. 

  971. 				PG_RETURN_NULL();
    972. 

  972. 			}
    973. 

  973. 			hasbandindex[i] = 1;
    974. 

  974. 		}
    975. 

  975. 		else
    976. 

  976. 			hasbandindex[i] = 0;
    977. 

  977. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    978. 

  978. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    979. 

  979. 		j++;
    980. 

  980. 	}
    981. 

  981. 

  982. 	/* distance */
    983. 

  983. 	if (PG_ARGISNULL(4)) {
    984. 

  984. 		elog(NOTICE, "Distance cannot be NULL.  Returning NULL");
    985. 

  985. 		for (k = 0; k < set_count; k++) {
    986. 

  986. 			rt_raster_destroy(rast[k]);
    987. 

  987. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    988. 

  988. 		}
    989. 

  989. 		PG_RETURN_NULL();
    990. 

  990. 	}
    991. 

  991. 

  992. 	distance = PG_GETARG_FLOAT8(4);
    993. 

  993. 	if (distance < 0) {
    994. 

  994. 		elog(NOTICE, "Distance cannot be less than zero.  Returning NULL");
    995. 

  995. 		for (k = 0; k < set_count; k++) {
    996. 

  996. 			rt_raster_destroy(rast[k]);
    997. 

  997. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    998. 

  998. 		}
    999. 

  999. 		PG_RETURN_NULL();
    1000. 

  1000. 	}
    1001. 

  1001. 

  1002. 	/* hasbandindex must be balanced */
    1003. 

  1003. 	if (
    1004. 

  1004. 		(hasbandindex[0] && !hasbandindex[1]) ||
    1005. 

  1005. 		(!hasbandindex[0] && hasbandindex[1])
    1006. 

  1006. 	) {
    1007. 

  1007. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    1008. 

  1008. 		for (k = 0; k < set_count; k++) {
    1009. 

  1009. 			rt_raster_destroy(rast[k]);
    1010. 

  1010. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1011. 

  1011. 		}
    1012. 

  1012. 		PG_RETURN_NULL();
    1013. 

  1013. 	}
    1014. 

  1014. 

  1015. 	/* SRID must match */
    1016. 

  1016. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    1017. 

  1017. 		for (k = 0; k < set_count; k++) {
    1018. 

  1018. 			rt_raster_destroy(rast[k]);
    1019. 

  1019. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1020. 

  1020. 		}
    1021. 

  1021. 		elog(ERROR, "The two rasters provided have different SRIDs");
    1022. 

  1022. 		PG_RETURN_NULL();
    1023. 

  1023. 	}
    1024. 

  1024. 

  1025. 	rtn = rt_raster_within_distance(
    1026. 

  1026. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    1027. 

  1027. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    1028. 

  1028. 		distance,
    1029. 

  1029. 		&result
    1030. 

  1030. 	);
    1031. 

  1031. 	for (k = 0; k < set_count; k++) {
    1032. 

  1032. 		rt_raster_destroy(rast[k]);
    1033. 

  1033. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1034. 

  1034. 	}
    1035. 

  1035. 

  1036. 	if (rtn != ES_NONE) {
    1037. 

  1037. 		elog(ERROR, "RASTER_dwithin: Could not test that the two rasters are within the specified distance of each other");
    1038. 

  1038. 		PG_RETURN_NULL();
    1039. 

  1039. 	}
    1040. 

  1040. 

  1041. 	PG_RETURN_BOOL(result);
    1042. 

  1042. }
    1043. 

  1043. 

  1044. /**
    1045. 

  1045.  * See if the two rasters are fully within the specified distance of each other
    1046. 

  1046.  */
    1047. 

  1047. PG_FUNCTION_INFO_V1(RASTER_dfullywithin);
    1048. 

  1048. Datum RASTER_dfullywithin(PG_FUNCTION_ARGS)
    1049. 

  1049. {
    1050. 

  1050. 	const int set_count = 2;
    1051. 

  1051. 	rt_pgraster *pgrast[2];
    1052. 

  1052. 	int pgrastpos[2] = {-1, -1};
    1053. 

  1053. 	rt_raster rast[2] = {NULL};
    1054. 

  1054. 	uint32_t bandindex[2] = {0};
    1055. 

  1055. 	uint32_t hasbandindex[2] = {0};
    1056. 

  1056. 	double distance = 0;
    1057. 

  1057. 

  1058. 	uint32_t i;
    1059. 

  1059. 	uint32_t j;
    1060. 

  1060. 	uint32_t k;
    1061. 

  1061. 	uint32_t numBands;
    1062. 

  1062. 	int rtn;
    1063. 

  1063. 	int result;
    1064. 

  1064. 

  1065. 	for (i = 0, j = 0; i < set_count; i++) {
    1066. 

  1066. 		/* pgrast is null, return null */
    1067. 

  1067. 		if (PG_ARGISNULL(j)) {
    1068. 

  1068. 			for (k = 0; k < i; k++) {
    1069. 

  1069. 				rt_raster_destroy(rast[k]);
    1070. 

  1070. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1071. 

  1071. 			}
    1072. 

  1072. 			PG_RETURN_NULL();
    1073. 

  1073. 		}
    1074. 

  1074. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(j));
    1075. 

  1075. 		pgrastpos[i] = j;
    1076. 

  1076. 		j++;
    1077. 

  1077. 

  1078. 		/* raster */
    1079. 

  1079. 		rast[i] = rt_raster_deserialize(pgrast[i], FALSE);
    1080. 

  1080. 		if (!rast[i]) {
    1081. 

  1081. 			for (k = 0; k <= i; k++) {
    1082. 

  1082. 				if (k < i)
    1083. 

  1083. 					rt_raster_destroy(rast[k]);
    1084. 

  1084. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1085. 

  1085. 			}
    1086. 

  1086. 			elog(ERROR, "RASTER_dfullywithin: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    1087. 

  1087. 			PG_RETURN_NULL();
    1088. 

  1088. 		}
    1089. 

  1089. 

  1090. 		/* numbands */
    1091. 

  1091. 		numBands = rt_raster_get_num_bands(rast[i]);
    1092. 

  1092. 		if (numBands < 1) {
    1093. 

  1093. 			elog(NOTICE, "The %s raster provided has no bands", i < 1 ? "first" : "second");
    1094. 

  1094. 			if (i > 0) i++;
    1095. 

  1095. 			for (k = 0; k < i; k++) {
    1096. 

  1096. 				rt_raster_destroy(rast[k]);
    1097. 

  1097. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1098. 

  1098. 			}
    1099. 

  1099. 			PG_RETURN_NULL();
    1100. 

  1100. 		}
    1101. 

  1101. 

  1102. 		/* band index */
    1103. 

  1103. 		if (!PG_ARGISNULL(j)) {
    1104. 

  1104. 			bandindex[i] = PG_GETARG_INT32(j);
    1105. 

  1105. 			if (bandindex[i] < 1 || bandindex[i] > numBands) {
    1106. 

  1106. 				elog(NOTICE, "Invalid band index (must use 1-based) for the %s raster. Returning NULL", i < 1 ? "first" : "second");
    1107. 

  1107. 				if (i > 0) i++;
    1108. 

  1108. 				for (k = 0; k < i; k++) {
    1109. 

  1109. 					rt_raster_destroy(rast[k]);
    1110. 

  1110. 					PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1111. 

  1111. 				}
    1112. 

  1112. 				PG_RETURN_NULL();
    1113. 

  1113. 			}
    1114. 

  1114. 			hasbandindex[i] = 1;
    1115. 

  1115. 		}
    1116. 

  1116. 		else
    1117. 

  1117. 			hasbandindex[i] = 0;
    1118. 

  1118. 		POSTGIS_RT_DEBUGF(4, "hasbandindex[%d] = %d", i, hasbandindex[i]);
    1119. 

  1119. 		POSTGIS_RT_DEBUGF(4, "bandindex[%d] = %d", i, bandindex[i]);
    1120. 

  1120. 		j++;
    1121. 

  1121. 	}
    1122. 

  1122. 

  1123. 	/* distance */
    1124. 

  1124. 	if (PG_ARGISNULL(4)) {
    1125. 

  1125. 		elog(NOTICE, "Distance cannot be NULL.  Returning NULL");
    1126. 

  1126. 		for (k = 0; k < set_count; k++) {
    1127. 

  1127. 			rt_raster_destroy(rast[k]);
    1128. 

  1128. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1129. 

  1129. 		}
    1130. 

  1130. 		PG_RETURN_NULL();
    1131. 

  1131. 	}
    1132. 

  1132. 

  1133. 	distance = PG_GETARG_FLOAT8(4);
    1134. 

  1134. 	if (distance < 0) {
    1135. 

  1135. 		elog(NOTICE, "Distance cannot be less than zero.  Returning NULL");
    1136. 

  1136. 		for (k = 0; k < set_count; k++) {
    1137. 

  1137. 			rt_raster_destroy(rast[k]);
    1138. 

  1138. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1139. 

  1139. 		}
    1140. 

  1140. 		PG_RETURN_NULL();
    1141. 

  1141. 	}
    1142. 

  1142. 

  1143. 	/* hasbandindex must be balanced */
    1144. 

  1144. 	if (
    1145. 

  1145. 		(hasbandindex[0] && !hasbandindex[1]) ||
    1146. 

  1146. 		(!hasbandindex[0] && hasbandindex[1])
    1147. 

  1147. 	) {
    1148. 

  1148. 		elog(NOTICE, "Missing band index.  Band indices must be provided for both rasters if any one is provided");
    1149. 

  1149. 		for (k = 0; k < set_count; k++) {
    1150. 

  1150. 			rt_raster_destroy(rast[k]);
    1151. 

  1151. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1152. 

  1152. 		}
    1153. 

  1153. 		PG_RETURN_NULL();
    1154. 

  1154. 	}
    1155. 

  1155. 

  1156. 	/* SRID must match */
    1157. 

  1157. 	if (rt_raster_get_srid(rast[0]) != rt_raster_get_srid(rast[1])) {
    1158. 

  1158. 		for (k = 0; k < set_count; k++) {
    1159. 

  1159. 			rt_raster_destroy(rast[k]);
    1160. 

  1160. 			PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1161. 

  1161. 		}
    1162. 

  1162. 		elog(ERROR, "The two rasters provided have different SRIDs");
    1163. 

  1163. 		PG_RETURN_NULL();
    1164. 

  1164. 	}
    1165. 

  1165. 

  1166. 	rtn = rt_raster_fully_within_distance(
    1167. 

  1167. 		rast[0], (hasbandindex[0] ? bandindex[0] - 1 : -1),
    1168. 

  1168. 		rast[1], (hasbandindex[1] ? bandindex[1] - 1 : -1),
    1169. 

  1169. 		distance,
    1170. 

  1170. 		&result
    1171. 

  1171. 	);
    1172. 

  1172. 	for (k = 0; k < set_count; k++) {
    1173. 

  1173. 		rt_raster_destroy(rast[k]);
    1174. 

  1174. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1175. 

  1175. 	}
    1176. 

  1176. 

  1177. 	if (rtn != ES_NONE) {
    1178. 

  1178. 		elog(ERROR, "RASTER_dfullywithin: Could not test that the two rasters are fully within the specified distance of each other");
    1179. 

  1179. 		PG_RETURN_NULL();
    1180. 

  1180. 	}
    1181. 

  1181. 

  1182. 	PG_RETURN_BOOL(result);
    1183. 

  1183. }
    1184. 

  1184. 

  1185. /**
    1186. 

  1186.  * See if two rasters are aligned
    1187. 

  1187.  */
    1188. 

  1188. PG_FUNCTION_INFO_V1(RASTER_sameAlignment);
    1189. 

  1189. Datum RASTER_sameAlignment(PG_FUNCTION_ARGS)
    1190. 

  1190. {
    1191. 

  1191. 	const int set_count = 2;
    1192. 

  1192. 	rt_pgraster *pgrast[2];
    1193. 

  1193. 	int pgrastpos[2] = {-1, -1};
    1194. 

  1194. 	rt_raster rast[2] = {NULL};
    1195. 

  1195. 

  1196. 	uint32_t i;
    1197. 

  1197. 	uint32_t j;
    1198. 

  1198. 	uint32_t k;
    1199. 

  1199. 	int rtn;
    1200. 

  1200. 	int aligned = 0;
    1201. 

  1201. 	char *reason = NULL;
    1202. 

  1202. 

  1203. 	for (i = 0, j = 0; i < set_count; i++) {
    1204. 

  1204. 		/* pgrast is null, return null */
    1205. 

  1205. 		if (PG_ARGISNULL(j)) {
    1206. 

  1206. 			for (k = 0; k < i; k++) {
    1207. 

  1207. 				rt_raster_destroy(rast[k]);
    1208. 

  1208. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1209. 

  1209. 			}
    1210. 

  1210. 			PG_RETURN_NULL();
    1211. 

  1211. 		}
    1212. 

  1212. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM_SLICE(PG_GETARG_DATUM(j), 0, sizeof(struct rt_raster_serialized_t));
    1213. 

  1213. 		pgrastpos[i] = j;
    1214. 

  1214. 		j++;
    1215. 

  1215. 

  1216. 		/* raster */
    1217. 

  1217. 		rast[i] = rt_raster_deserialize(pgrast[i], TRUE);
    1218. 

  1218. 		if (!rast[i]) {
    1219. 

  1219. 			for (k = 0; k <= i; k++) {
    1220. 

  1220. 				if (k < i)
    1221. 

  1221. 					rt_raster_destroy(rast[k]);
    1222. 

  1222. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1223. 

  1223. 			}
    1224. 

  1224. 			elog(ERROR, "RASTER_sameAlignment: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    1225. 

  1225. 			PG_RETURN_NULL();
    1226. 

  1226. 		}
    1227. 

  1227. 	}
    1228. 

  1228. 

  1229. 	rtn = rt_raster_same_alignment(
    1230. 

  1230. 		rast[0],
    1231. 

  1231. 		rast[1],
    1232. 

  1232. 		&aligned,
    1233. 

  1233. 		&reason
    1234. 

  1234. 	);
    1235. 

  1235. 	for (k = 0; k < set_count; k++) {
    1236. 

  1236. 		rt_raster_destroy(rast[k]);
    1237. 

  1237. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1238. 

  1238. 	}
    1239. 

  1239. 

  1240. 	if (rtn != ES_NONE) {
    1241. 

  1241. 		elog(ERROR, "RASTER_sameAlignment: Could not test for alignment on the two rasters");
    1242. 

  1242. 		PG_RETURN_NULL();
    1243. 

  1243. 	}
    1244. 

  1244. 

  1245. 	/* only output reason if not aligned */
    1246. 

  1246. 	if (reason != NULL && !aligned)
    1247. 

  1247. 		elog(NOTICE, "%s", reason);
    1248. 

  1248. 

  1249. 	PG_RETURN_BOOL(aligned);
    1250. 

  1250. }
    1251. 

  1251. 

  1252. /**
    1253. 

  1253.  * Return a reason why two rasters are not aligned
    1254. 

  1254.  */
    1255. 

  1255. PG_FUNCTION_INFO_V1(RASTER_notSameAlignmentReason);
    1256. 

  1256. Datum RASTER_notSameAlignmentReason(PG_FUNCTION_ARGS)
    1257. 

  1257. {
    1258. 

  1258. 	const int set_count = 2;
    1259. 

  1259. 	rt_pgraster *pgrast[2];
    1260. 

  1260. 	int pgrastpos[2] = {-1, -1};
    1261. 

  1261. 	rt_raster rast[2] = {NULL};
    1262. 

  1262. 

  1263. 	uint32_t i;
    1264. 

  1264. 	uint32_t j;
    1265. 

  1265. 	uint32_t k;
    1266. 

  1266. 	int rtn;
    1267. 

  1267. 	int aligned = 0;
    1268. 

  1268. 	char *reason = NULL;
    1269. 

  1269. 	text *result = NULL;
    1270. 

  1270. 

  1271. 	for (i = 0, j = 0; i < set_count; i++) {
    1272. 

  1272. 		/* pgrast is null, return null */
    1273. 

  1273. 		if (PG_ARGISNULL(j)) {
    1274. 

  1274. 			for (k = 0; k < i; k++) {
    1275. 

  1275. 				rt_raster_destroy(rast[k]);
    1276. 

  1276. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1277. 

  1277. 			}
    1278. 

  1278. 			PG_RETURN_NULL();
    1279. 

  1279. 		}
    1280. 

  1280. 		pgrast[i] = (rt_pgraster *) PG_DETOAST_DATUM_SLICE(PG_GETARG_DATUM(j), 0, sizeof(struct rt_raster_serialized_t));
    1281. 

  1281. 		pgrastpos[i] = j;
    1282. 

  1282. 		j++;
    1283. 

  1283. 

  1284. 		/* raster */
    1285. 

  1285. 		rast[i] = rt_raster_deserialize(pgrast[i], TRUE);
    1286. 

  1286. 		if (!rast[i]) {
    1287. 

  1287. 			for (k = 0; k <= i; k++) {
    1288. 

  1288. 				if (k < i)
    1289. 

  1289. 					rt_raster_destroy(rast[k]);
    1290. 

  1290. 				PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1291. 

  1291. 			}
    1292. 

  1292. 			elog(ERROR, "RASTER_notSameAlignmentReason: Could not deserialize the %s raster", i < 1 ? "first" : "second");
    1293. 

  1293. 			PG_RETURN_NULL();
    1294. 

  1294. 		}
    1295. 

  1295. 	}
    1296. 

  1296. 

  1297. 	rtn = rt_raster_same_alignment(
    1298. 

  1298. 		rast[0],
    1299. 

  1299. 		rast[1],
    1300. 

  1300. 		&aligned,
    1301. 

  1301. 		&reason
    1302. 

  1302. 	);
    1303. 

  1303. 	for (k = 0; k < set_count; k++) {
    1304. 

  1304. 		rt_raster_destroy(rast[k]);
    1305. 

  1305. 		PG_FREE_IF_COPY(pgrast[k], pgrastpos[k]);
    1306. 

  1306. 	}
    1307. 

  1307. 

  1308. 	if (rtn != ES_NONE) {
    1309. 

  1309. 		elog(ERROR, "RASTER_notSameAlignmentReason: Could not test for alignment on the two rasters");
    1310. 

  1310. 		PG_RETURN_NULL();
    1311. 

  1311. 	}
    1312. 

  1312. 

  1313. 	result = cstring2text(reason);
    1314. 

  1314. 	PG_RETURN_TEXT_P(result);
    1315. 

  1315. }
    1316. 

  1316.