grdmask: fix region '-Rg' not recognising longitudes more than 180 (#9001)

Fix #8430

Assisted-by: Claude Sonnet 4.6

Author: joa-quim

src/gmt_support.c

@@ -12313,21 +12313,21 @@ int gmt_get_format (struct GMT_CTRL *GMT, double interval, char *unit, char *pre
 	return (ndec);
 }
 
-void gmt_set_inside_mode (struct GMT_CTRL *GMT, struct GMT_DATASET *D, unsigned int mode) {
+void gmt_set_inside_mode(struct GMT_CTRL *GMT, struct GMT_DATASET *D, unsigned int mode) {
 	/* Determine if we use spherical or Cartesian function for in--on-out polygon tests */
 	static char *method[2] = {"Cartesian", "spherical"};
 	if (mode == GMT_IOO_SPHERICAL)	/* Force spherical */
 		GMT->current.proj.sph_inside = true;
 	else if (mode == GMT_IOO_CARTESIAN)	/* Force Cartesian */
 		GMT->current.proj.sph_inside = false;
-	else if (gmt_M_is_cartesian (GMT, GMT_IN))	/* If data is Cartesian then we do that */
+	else if (gmt_M_is_cartesian(GMT, GMT_IN))	/* If data is Cartesian then we do that */
 		GMT->current.proj.sph_inside = false;
 	else if (GMT->current.map.is_world)	/* Here we are dealing with geographic data that has 360 degree range */
 		GMT->current.proj.sph_inside = true;
 	else if (D) {	/* Geographic data less than 360 degree range in longitudes */
 		double lat[2];
 		lat[0] = D->min[GMT_Y]; lat[1] = D->max[GMT_Y];
-		if (doubleAlmostEqual (lat[0], -90.0) || doubleAlmostEqual (lat[1], +90.0))	/* Goes to a pole, must do spherical */
+		if (doubleAlmostEqual(lat[0], -90.0) || doubleAlmostEqual(lat[1], +90.0))	/* Goes to a pole, must do spherical */
 			GMT->current.proj.sph_inside = true;
 		else {	/* Limited in lon and lat, can do Cartesian but must ensure polygons do not jump within range */
 			uint64_t tbl, seg, row;
@@ -12347,13 +12347,22 @@ void gmt_set_inside_mode (struct GMT_CTRL *GMT, struct GMT_DATASET *D, unsigned
 					S = D->table[tbl]->segment[seg];	/* Shorthand */
 					for (row = 0; row < S->n_rows; row++)
 						gmt_lon_range_adjust (range, &S->data[GMT_X][row]);
+					/* Recompute segment x-range from adjusted vertices; S->min/max may be in a
+					 * different longitude convention (e.g. 0-360 from quad_finalize) than the
+					 * vertices just adjusted above, causing gmt_non_zero_winding to receive x
+					 * and polygon vertices in mismatched ranges. */
+					S->min[GMT_X] = S->max[GMT_X] = S->data[GMT_X][0];
+					for (row = 1; row < S->n_rows; row++) {
+						if (S->data[GMT_X][row] < S->min[GMT_X]) S->min[GMT_X] = S->data[GMT_X][row];
+						if (S->data[GMT_X][row] > S->max[GMT_X]) S->max[GMT_X] = S->data[GMT_X][row];
+					}
 				}
 			}
 		}
 	}
 	else
-		GMT_Report (GMT->parent, GMT_MSG_ERROR, "Not enough information given to gmt_set_inside_mode.\n");
-	GMT_Report (GMT->parent, GMT_MSG_INFORMATION, "A point's inside/outside status w.r.t. polygon(s) will be determined using a %s algorithm.\n", method[GMT->current.proj.sph_inside]);
+		GMT_Report(GMT->parent, GMT_MSG_ERROR, "Not enough information given to gmt_set_inside_mode.\n");
+	GMT_Report(GMT->parent, GMT_MSG_INFORMATION, "A point's inside/outside status w.r.t. polygon(s) will be determined using a %s algorithm.\n", method[GMT->current.proj.sph_inside]);
 }
 
 /*! . */
@@ -13238,22 +13247,7 @@ int gmt_grd_BC_set (struct GMT_CTRL *GMT, struct GMT_GRID *G, unsigned int direc
 }
 
 /* Clipper to ensure a byte stays in 0-255 range */
-#if 1
-GMT_LOCAL inline unsigned char gmtsupport_clip_to_byte (int byte) { if (byte < 0) return (0); else if (byte > 255) return (255); else return ((unsigned char)byte);}
-#else
-/* For debugging BC actions for images when the result is outside byte range  */
-GMT_LOCAL unsigned char gmtsupport_clip_to_byte (int byte) {
-	if (byte < 0) {
-		fprintf (stderr, "byte = %d\n", byte);
-		return (0);
-	}
-	if (byte > 255) {
-		fprintf (stderr, "byte = %d\n", byte);
-		return (255);
-	}
-	return ((unsigned char)byte);
-}
-#endif
+GMT_LOCAL inline unsigned char gmtsupport_clip_to_byte(int byte) { if (byte < 0) return (0); else if (byte > 255) return (255); else return ((unsigned char)byte);}
 
 /*! . */
 int gmtlib_image_BC_set (struct GMT_CTRL *GMT, struct GMT_IMAGE *I) {