Re: Creating an imcountour for a bpt [message #94455 is a reply to message #94451] |
Thu, 25 May 2017 07:24  |
wlandsman
Messages: 743 Registered: June 2000
|
Senior Member |
|
|
You show a full program and say you are having difficulty but you don't say what the difficulty is. Are you getting an error message? Does the plot not show up? Anyway here are some suggestions?
1. I am not certain where you are getting the program TVSCALE, but it is likely a very old program from David Fanning's library that was replaced by CGIMAGE (http://www.idlcoyote.com/programs/cgimage.pro). After 2011 all improvements, and bug fixes only occurred in CGIMAGE.
2. I also suspect that you getting IMCONTOUR from https://idlastro.gsfc.nasa.gov/ftp/pro/astro/imcontour.pro Make sure you have a recent version.
3. Why are you setting NLEVELS = 0? Do you want to show a CONTOUR plot with zero levels? If you don't want contours, and are just using imcontour.pro for astronomical labeling then it is better to set the /nodata keyword.
--Wayne
On Wednesday, May 24, 2017 at 3:55:12 PM UTC-4, Rafael Cirolini wrote:
> loadct, 12
> tvscale, mask, /KEEP_ASPECT_RATIO, POSITION=p, margin=0.001
> loadct, 6
> imcontour, mask, hdr, nlevels=0, /Noerase, min_value=500, max_value=501, TYPE=0, xminor=-1.5, yminor=-1.5, position=p
>
>
> device, /close
>
> end
>
>
> First the program makes a mask on the central pixel of the cube and does the operations of the emission lines that are the characteristic of the diagram. Right after the mascara I'm trying to create to first create the tvscale and then the imcountour, which is exactly where I'm having difficulty. I think this part is probably wrong because it's my first time creating an imcountour. The equations of the curves where the diagram is defined are these:
|
|
|