Function gmtplot_draw_dir_rose

src/gmt_plot.c:3584–3755

#define ROSE_OFFSET_SCL		0.4	/* Set label offset to 40% of font size */

GMT_LOCAL void gmtplot_draw_dir_rose (struct GMT_CTRL *GMT, struct PSL_CTRL *PSL, struct GMT_MAP_ROSE *mr) {
	bool adjust = (GMT->current.setting.map_embellishment_mode);	/* A bit shorter to use */
	unsigned int i, kind, form, just[4] = {PSL_TC, PSL_ML, PSL_BC, PSL_MR};
	int k;
	double angle, title_size, txt_offset, L[4], R[4], x[PSL_MAX_DIMS], y[PSL_MAX_DIMS], xp[PSL_MAX_DIMS], yp[PSL_MAX_DIMS], tx[3], ty[3];
	double lon, lat, s, c, rot[4] = {0.0, 45.0, 22.5, -22.5};
	struct GMT_FILL f;

	/* Set default font size in modern mode based on rose size - or honor setting in classic mode */
	title_size = (adjust) ? ROSE_LABEL_SCL * mr->size * 72.0 : GMT->current.setting.font_title.size;

	/* Initialize fill structure */
	gmt_init_fill (GMT, &f, GMT->current.setting.color_patch[GMT_BGD][0], GMT->current.setting.color_patch[GMT_BGD][1], GMT->current.setting.color_patch[GMT_BGD][2]);

	gmt_xy_to_geo (GMT, &lon, &lat, mr->refpoint->x, mr->refpoint->y);
	angle = gmt_azim_to_angle (GMT, lon, lat, DIST_TO_2ND_POINT, 90.0);	/* Get angle of E-W direction at this location */

#if 0
	/* This section can be activated if we need to test for angles of directional roses.
	 * Then, use export GMT_FAKE_ANGLE=value to change the orientation of the rose. */
#ifdef DEBUG
	{
		char *this_c = NULL;
		if ((this_c = getenv ("GMT_FAKE_ANGLE")) != NULL)
			angle = atof (this_c);
	}
#endif
#endif

	if (angle < 0.0) angle += 360.0;	if (angle >= 360.0) angle -= 360.0;
	gmt_setpen (GMT, &GMT->current.setting.map_tick_pen[GMT_PRIMARY]);

	PSL_command (PSL, "V\n");	/* Place the rose under gsave/grestore */

	if (mr->type == GMT_ROSE_DIR_FANCY) {	/* Fancy scale */
		PSL_comment (PSL, "Draw fancy directional rose of level %d\n", mr->kind);
		/* Set default offset size in modern mode based on rose size - or honor setting in classic mode */
		txt_offset = (adjust) ? ROSE_OFFSET_SCL * title_size / 72.0 : GMT->current.setting.map_title_offset;
		if (mr->do_label) {	/* Need to determine the shift in (x,y) due to the dimensions of label and offset */
			int hh = mr->justify % 4 - 2, vv = mr->justify / 4 - 1;	/* Horizontal and vertical shift indicators */
			char *T_label = NULL;
			if (hh) {	/* We need horizontal adjustments. Do label offset here and adjust for text size in PSL */
				T_label = (hh == -1) ? mr->label[3] : mr->label[1];	/* Get the W or E text label */
				if (T_label[0]) {	/* Get width of label and use it to shift origin left or right */
					mr->refpoint->x -= hh * txt_offset;	/* Any horizontal shifts we know exactly */
					PSL_deftextdim (PSL, "-w", title_size, T_label);	/* Get label width and place on PSL stack */
					PSL_command (PSL, "%d mul 0 T\n", -hh);	/* Multiply by +1 or -1 and transform x-origin by that amount */
				}
			}
			if (vv) {	/* We need vertical adjustments. Do label offset here and adjust for text size in PSL */
				T_label = (vv == -1) ? mr->label[0] : mr->label[2];	/* Get the S or N text label */
				if (T_label[0]) {	/* Get height of string and use it to shift origin up or down */
					mr->refpoint->y -= vv * txt_offset;	/* Any horizontal shifts we know exactly */
					PSL_deftextdim (PSL, "-H", title_size, T_label);	/* Get label height and place on PSL stack */
					PSL_command (PSL, "%d mul 0 exch T\n", -vv);	/* Multiply by +1 or -1 and transform y-origin by that amount */
				}
			}
		}