comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » IDL ROI Objects

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: IDL ROI Objects [message #90131 is a reply to message #90130]

Thu, 29 January 2015 02:50

Fabzi
Messages: 305
Registered: July 2010

Senior Member

On 29.01.2015 11:05, David B wrote:
> IDl>Are you looking at the floor or the sky?

For the record, a better "tutorial" here:

PRO ROIPolygonAlternativeView

poly = [[0.8, 1.7, 2.7, 2.7, 0.8], [0.2, 2.2, 2.2, 0.2, 0.2]]

nn = 4
centersX = transpose(fltarr(nn)+1) ## (INDGEN(nn)+ 0.5)
centersY = transpose(INDGEN(nn)+ 0.5) ## (fltarr(nn)+1)

TVLCT, cgColor('red6', /Triple), 255
TVLCT, cgColor('white', /Triple), 1

cgDisplay, 700, 800
pos = cgLayout([2,2], OXMARGIN=[2, 2], OYMARGIN=[3, 4], XGAP=2, YGAP=6)

; Rule 0
p = OBJ_NEW('IDLanROI', poly[*,0], poly[*,1])
mask = p->ComputeMask(DIMENSIONS=[nn, nn], MASK_RULE=0)
anImage = BytArr(nn,nn)+1B
anImage[Where(mask EQ 255)] = 255
po = pos[*,0]
cgImage, anImage, /AXES, /SAVE, POS=po, TITLE='RULE 0', /NOERASE
cgPlotS, centersX, centersY, PSYM=16, /DATA
cgPlotS, poly[*, 0]+ 0.5, poly[*, 1]+ 0.5, COLOR='green', /DATA, THICK=3
cgPlot, [1], XRange=[0,4], YRange=[0,4], /NoData, $
YTickLen=1.0, XTickLen=1.0, POS=po, /NOERASE

; Rule 1
p = OBJ_NEW('IDLanROI', poly[*,0], poly[*,1])
mask = p->ComputeMask(DIMENSIONS=[nn, nn], MASK_RULE=1)
anImage = BytArr(nn,nn)+1B
anImage[Where(mask EQ 255)] = 255
po = pos[*,1]
cgImage, anImage, /AXES, /SAVE, POS=po, TITLE='RULE 1', /NOERASE
cgPlotS, centersX, centersY, PSYM=16, /DATA
cgPlotS, poly[*, 0]+ 0.5, poly[*, 1]+ 0.5, COLOR='green', /DATA, THICK=3
cgPlot, [1], XRange=[0,4], YRange=[0,4], /NoData, $
YTickLen=1.0, XTickLen=1.0, POS=po, /NOERASE

; Rule 2
p = OBJ_NEW('IDLanROI', poly[*,0], poly[*,1])
mask = p->ComputeMask(DIMENSIONS=[nn, nn], MASK_RULE=2)
anImage = BytArr(nn,nn)+1B
anImage[Where(mask EQ 255)] = 255
po = pos[*,2]
cgImage, anImage, /AXES, /SAVE, POS=po, TITLE='RULE 2', /NOERASE
cgPlotS, centersX, centersY, PSYM=16, /DATA
cgPlotS, poly[*, 0]+ 0.5, poly[*, 1]+ 0.5, COLOR='green', /DATA, THICK=3
cgPlot, [1], XRange=[0,4], YRange=[0,4], /NoData, $
YTickLen=1.0, XTickLen=1.0, POS=po, /NOERASE

; Contains points
p = OBJ_NEW('IDLanROI', poly[*,0], poly[*,1])
mask = p->ContainsPoints(centersX[*]-0.5, centersY[*]-0.5)
anImage = BytArr(nn,nn)+1B
p = where(mask eq 1)
po = pos[*,3]
cgImage, anImage, /AXES, /SAVE, POS=po, TITLE='CONTAINS POINTS', /NOERASE
cgPlotS, centersX, centersY, PSYM=16, /DATA
cgPlotS, centersX[p], centersY[p], PSYM=16, /DATA, COLOR='red'
cgPlotS, poly[*, 0]+ 0.5, poly[*, 1]+ 0.5, COLOR='green', /DATA, THICK=3

end

Report message to a moderator

[Message index]

		IDL ROI Objects By: David B on Tue, 27 January 2015 09:28
		Re: IDL ROI Objects By: David Fanning on Tue, 27 January 2015 10:53
		Re: IDL ROI Objects By: David Fanning on Tue, 27 January 2015 11:00
		Re: IDL ROI Objects By: Fabzi on Tue, 27 January 2015 15:10
		Re: IDL ROI Objects By: David B on Wed, 28 January 2015 02:00
		Re: IDL ROI Objects By: David B on Wed, 28 January 2015 02:06
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 06:12
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 06:19
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 09:17
		Re: IDL ROI Objects By: Fabzi on Wed, 28 January 2015 09:49
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 10:17
		Re: IDL ROI Objects By: Fabzi on Wed, 28 January 2015 10:30
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 10:42
		Re: IDL ROI Objects By: Fabzi on Wed, 28 January 2015 10:56
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 11:20
		Re: IDL ROI Objects By: Fabzi on Wed, 28 January 2015 11:24
		Re: IDL ROI Objects By: David Fanning on Wed, 28 January 2015 11:38
		Re: IDL ROI Objects By: David B on Thu, 29 January 2015 02:05
		Re: IDL ROI Objects By: Fabzi on Thu, 29 January 2015 02:50

Previous Topic:	Krig2d from point to image
Next Topic:	Re: Why did IDL just start crashing

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 17:30:06 PDT 2025

Total time taken to generate the page: 1.60012 seconds