All new activity

Copy wave and paste as new wave with arithmatic option

supra — Tue, 15 May 2012 14:30:42 -0700

Hi
I wanted to have a menu for "wave copy" from a existed table and then it will be pasted as a "new wave" in the same table where "new wave" = "old wave" +-x/ "new wave" and I should have GUI to select +,-,x, or /. Can anybody tell me

a) How I can add checkboxes in the GUI to select +-x/ .
b) How I can keep the GUI panel still visible even after I click "ok" and a option to "Undo" the operation that I have done.

#pragma rtGlobals=1		// Use modern global access method.
 
Menu "Column Math"
      "Column MathA/1", copy22(0)
End
 
Function copy22(AorB)
 
Variable AorB        // 0 for A, 1 for B
 
      String tableName = WinName(0, 2,1)
      if (strlen(tableName) == 0)
              Abort  "There are no tabels"
      endif
 
GetSelection table, $tableName, 1      //Get selection row and column indices
 
       Variable startColumn = V_startCol               //Index of first selected column 
 
Wave y2 = WaveRefindexed(tableName, startColumn, 1)
 
Duplicate/O y2, wave0_y2
 
Variable term=1
 
	Prompt term, "Enter parameter: "		// Set prompt for XRange param
	DoPrompt "Enter parameter", term
 
	if (V_Flag)
		return -1								// User canceled
endif
 
wave0_y2 = wave0_y2 + term
AppendToTable wave0_y2
 
     End

3D - Histogram creation help needed

channagiri.1 — Fri, 11 May 2012 06:51:25 -0700

I have a batch of thousand bmp images that are labelled in sequence (file name followed by a number). I need to make a 3D histogram for all these images combined and then do statistical analysis on this. By 3D histogram, I mean the histogram corresponding to each image in the batch all stacked together to form a surface plot. Any advice would be appreciated.

Easy way to mask waves?

amellnik — Mon, 02 Jan 2012 11:41:44 -0800

It seems like there should be a really easy way to mask waves, but I haven't been able to figure one out so far. For example, I might have to waves of the same length and I would like to delete all the points from the first wave where the corresponding point on the second wave is zero. The way I've been doing this so far is embarrassingly hackey:

targetwave = targetwave/(maskwave==0) WaveTransform zapINFs targetwave WaveTransform zapNaNs targetwave // 0/0 = NaN

This has the obvious problem that it gets rid of any NaNs and INFs that should be in targetwave. What's the proper way to do this? Thanks

-Alex

Greek Letters in Legend/Axis Titles

kathreen — Wed, 09 May 2012 13:18:20 -0700

I know you can insert Greek letters in Igor using the Symbol font, but it would be really awesome for all Greek letters (both uppercase and lowercase) to be available in Special-> Character. In a scientific field we use them a lot as a one-off for a word (ie D_alpha.)

Save wave "tweak" preferences; perform an action for each loaded wave

hmayes — Wed, 16 May 2012 13:11:43 -0700

Hello!

Two questions:

1) How do I make using the "space" delimiter character a default? I currently go to Load Data > Tweaks and check this option, but Igor resets to the default (unchecked "spaces") when I reopen.

2) I will many times want to open a data file with different numbers of columns, and then perform the same actions on every loaded wave (histogram, then graph). The wave names will be read from the first line of the data file, and will change from one file to the next. Is there a syntax for doing someone on every wave loaded?

Thank you,
-Heather

Commentize/Decommentize Keyboard Shortcuts

jjweimer — Sun, 22 Nov 2009 09:47:20 -0800

Having Command-{ and Command-} (with corresponding Windows equivalents) as shortcuts to "Commentize" and "Decommentize" procedure code blocks would be useful. (Also, I have always preferred the menu commands Comment and Uncomment as found in TeXShop on MacOS over the above phrasings).

Include with Adopt

agarttha — Wed, 16 May 2012 02:33:46 -0700

I have shared procedure files on the network that need to be automatically loaded and incorporated (adopted) in the experiment. The standard #include directive does not do this. Manually, I can adopt a procedure using "File", "Adopt Procedure...". But how can do this automatically? Does this functionality exist? Is there a better way to do what I want: import a procedure file and sever the connection to the original?

If this feature does not exist, could you create a special "adopt" flag to the #include directive?

Thanks.

Cheers

Editing single point in matrix

Jonpope — Thu, 01 Mar 2012 02:33:35 -0800

Hello

Before reading on I am aware this is an extremely basic question, however after consulting the Igor Manual and this forum I've been unable to find an answer to it.

The function below is intended to read through all the data points in a matrix and change any values above 100000 (which indicates a dead pixel in the detector) to zero/NaN.
However I have been unable to find a command that can edit a single point; when doing this manually the command window displayed 'Edit Wavename.Id' which sadly is of no use in a function.

Function Removedeadpixels(thewave)
Wave thewave
variable xx, yy, zz
xx=0
do
yy=0
do
zz= thewave [xx][yy]
if(zz>=100000)
zz=0//Syntax for editing a data point?
endif
yy=yy+1
while (yy<=487)
xx=xx+1
while (xx<=195)
End

Also in it's current form would the function save the changes in the wave if it is in a string when used in other functions?

Thanks in advance
Jon

Strange command-line printing behavior

s.r.chinn — Thu, 10 May 2012 04:13:46 -0700

I have just reinstalled IP 6.22A on a new Windows 7 PC, and see some strange printing behavior in the History/Command-Line. I want to print the value of a global variable, Gwidth, that happens to have value 2. It appears like:

•print Gwidth
  2int Gwidth
•printf  "%g\r", Gwidth
  2intf "%/g\r", Gwidth

What is also odd, is that when I copied and pasted the History lines into the above extract, it appeared correct. I edited it to show what my PC monitor displays. It appears that the monitor output result is the echoed print command, over-written by the variable value. Has anyone else encountered this?

Unclear on coding structure to optimize fit with Main event loop speed

fltcoils — Mon, 14 May 2012 10:26:53 -0700

I realize that IGOR has evolved - as has desktop computer hardware. Certainly the IGOR program must be controlled from hogging too much CPU time or getting ahead of windows routines to avoid problems. That is what something such are main loop speed is about, yes?

Current desktop systems with multicore 3ghz CPUs, SATA3 and USB3 are able to process input and deliver output many times faster than 15 years ago. Given that many users are able to work from such modern systems; can there be a way to tailor the loop speed for newer hardware and modern OSes? A setup value in the registry or such that can be altered to allow for a faster loop speed would have some benefit.

My application evolved from analysis of spectral waveform data into a more comprehensive control system which uses the convenient VISA functions to arm and read data from scopes, thus to also acquire the data. Over time this easy and helpful system has shown much utility. There is, of course, a constant goal of increasing capability, and thus acquisition speed.

Reading about feedback control problems (http://www.igorexchange.com/node/2624) I've realize I don't really have a way to understand how IGORs "main control loop speed" affects the way code in the .ipf files is written and executed. If something obvious such as avoiding do/while for/next loops would increase execution speed 100 times it certainly would make users happy.

I can guess a number of possibilities for what main loop speed is/does:
Speed at idle at which IGOR checks for user input
speed at which subroutines are triggered when called
speed at which XOPs are triggered when called
speed of executing individual command lines
speed of do/while for/next loops

Many coding options, subroutines, do/while loops, are written mostly for legibility rather than necessity. If structuring to use subroutines, or do/while for/next loops causes everything to pause and wait for another 20ms main loop to trigger, then I would want to change how I code in certain time critical places.

Any direct information or links would certainly be appreciated.

Digital map of the world

jost — Mon, 14 May 2012 23:04:33 -0700

In the GIS utilities a link and related input format support is given for a digital map of the world that does not appear to be accessible anymore.
I would be most grateful if someone could point me to a link to digital map information in vector form (political boundaries, rivers, cities etc) that can be put into IGOR via its GIS utilities. Direct access to any such data sets would be highly appreciated, likewise input routines for GIS files with current formats such as shapefiles.
Thanks very much in advance
Jost

Update traces using a popup menu

temasi — Sun, 13 May 2012 20:21:54 -0700

Hi,

Does anyone have a procedure to update the Y wave in a graph using a popup menu (subwindow)? I have tried using ReplaceWave in the procedure of the popupMenu control, but that doesn't work more than one time, as the current wave would be different from the original wave when I try to change it a second time.

Any suggestion will be appreciated.

Extracting a matrix from within a matrix

Jonpope — Mon, 14 May 2012 07:27:31 -0700

Hello

I am attempting to extract a square area (containing a peak) from a large 2D wave as a secondary wave.

Thus far I have simply been using the total intensity within the square, which is defined by a number of pixels around a selected point e.g. PeakPosP+15,PeakPosP-15, PeakPosQ+15, PeakPosQ-15 with (PeakPosP, PeakPosQ) being selected either manually or using the V_maxrowloc, V_maxcolloc variables from the ImageStats command.

However I am now looking to carry out curve fitting on the area contained within the square and therefore need to keep the point values intact.

Is there any command with which I could duplicate a selected portion of the matrix? Or would it be a case of using the Make command allocating values for all the points afterward?

As always your help is greatly appreciated.

Regards
Jon

Manipulating Gizmo Window Recreation Macros

jmcclimo2003 — Thu, 10 May 2012 08:36:34 -0700

I have Igor 6.2. I created a Gizmo window showing an image in the Z=0 plane and a movable spectroscopy slice in the vertical. Now I have to do the same thing for a bunch of others. It seems like the command line can't do everything in Gizmo like select all the various options in the slicer dialogues (and if I'm wrong please correct me since it will make the following moot; I can just put all the commands in the Igor Text File.) I noticed that despite this, Igor can create Gizmo window recreation macros that reproduce everything faithfully but there's doesn't seem to be a way to access the recreation macros for Gizmo. I saved an unpacked copy of the file to see if I could find it there but no such luck. If I could have access to this I'd just need to change the name of two waves and I'd be good.

Is this possible?

How to make scatter plot with markers not overlapping

wwzhang2 — Thu, 30 Sep 2010 14:16:46 -0700

When I tried to plot a x-y 2D wave graph the marker always overlap when they have a same value. However, I hope to show the data in a way as the image attached showing that is the markers for the same or very close value won't totally overlap but scattered horizontally. Can I do this in Igor Pro?

Problem with curve fit

millot1 — Fri, 11 May 2012 16:49:38 -0700

I would like to reproduce this result, i e have the Covariance matrix parameters printed after completion of curvefit

//¥CurveFit/M=2 Power Temperature[3,] /X=T_Us /W=Temperature_Err /I=1 /D /R /F={0.683000, 1}
// Fit converged properly
// Curve fit with data subrange:
// Row: 3, 30
// fit_Temperature= W_coef[0]+W_coef[1]*x^W_coef[2]
// Res_Temperature= Temperature - (W_coef[0]+W_coef[1]*T_Us^W_coef[2])
// W_coef={1421.9,4.3185,2.9768}
// V_chisq= 31.0292; V_npnts= 28; V_numNaNs= 0; V_numINFs= 0;
// V_startRow= 3; V_endRow= 30; V_startCol= 0; V_endCol= 0;
// W_sigma={394,0.909,0.0656}
// Coefficient values ± one standard deviation
// y0 = 1421.9 ± 394
// A = 4.3185 ± 0.909
// pow = 2.9768 ± 0.0656
// M_Covar[0][0]= {155471,-340.136,24.3494}
// M_Covar[0][1]= {-340.136,0.825423,-0.059566}
// M_Covar[0][2]= {24.3494,-0.059566,0.00430317}

however
even using /M=2 i do not the M_Covar results, even though a can edit M_Covar. Is there a way to specify we want them print?

I tried
CurveFit/M=2 HillEquation FS_Refl (,35)/X=Us (,35)/W=FS_Refl_Err(,35) /I=1 /D
but i got

Fit converged properly
Curve fit with data subrange:
FS_Refl[0,35]
fit_FS_Refl= W_coef[0]+(W_coef[1]-W_coef[0])*(x^W_coef[2]/(x^W_coef[2]+W_coef[3]^W_coef[2]))
W_coef={0.0064212,0.34443,8.033,13.789}
V_chisq= 1.2633;V_npnts= 36;V_numNaNs= 0;V_numINFs= 0;
V_startRow= 0;V_endRow= 35;
W_sigma={0.00379,0.0131,1.12,0.324}
Coefficient values ± one standard deviation
base =0.0064212 ± 0.00379
max =0.34443 ± 0.0131
rate =8.033 ± 1.12
xhalf =13.789 ± 0.324

Many thanks

Double axis

gzvitiello — Fri, 11 May 2012 13:10:22 -0700

Hello,

I would like to, after inserting another Y axis, to have the number labels outside the graph box. I will add a picture for a better understanding of what I'm trying to say, I would like to move the content inside the black dashed box to where the dashed line is.

Thanks in advance,

Gabriel

Creating N number of waves

isk8 — Fri, 11 May 2012 14:39:02 -0700

I'm trying to make a procedure that automates the following process and I got stuck.
I have waves that have a multiple of 1024 points (1024 x N points).
I'm trying to decompose these waves into N number of waves, each with 1024 points, whose names are wave_N.
I know how to decompose waves having 1024 points each, but I'm not sure how to make N number of waves depending on the number of points I have in the starting wave (1024 x N points).
Does anyone know how you would make N number of waves in this case?
Thank you.

Bootstrapping a Package Loader

jtigor — Thu, 03 May 2012 08:00:09 -0700

The folowing code is the content of an Independent Module used to load/unload a package. It also creates a menu item that changes between "Load Package" and "Unload Package" based on the existance of a global symbol when the module is compiled. This global symbol is created on package load and destroyed on package unload. The idea was to enable the automatic creation of the appropriate menu when a new experiment was created or when an old one was opened and the mechanics of the package loading/unloading could be completely invisible to the user, other than the menu item. The package loader is always loaded when Igor starts, via a shortcut in the Igor Procedures directory.

This works fine within an Igor session but does not work if Igor is shut down and an old experiment opened again. I thought that I had demonstrated to myself that this method of detecting whether the package was loaded or not worked when I devised it several months ago. The probem, as far as I can tell, is that the global symbol is not saved with the experiment. Is this the expected Igor behavior? Did I miss something in reading about symbols?

If the global symbol idea is not viable, what alternative methods are there? At the present time, the package does not require a package data folder and I wanted to avoid creating one just to remember the package load/unload state. This certainly could be done. A dummy procedure could be loaded as part of the package and that could be searched for to determine whether or not the package was loaded. Is it possible to determine whether user defined menus exist? I see there are several functions for working with built in menus, but haven't been able to find anything related to user defined menus. This might server as a means to detecte the package state. Any ideas or opinions?

System is Win7, Igor V6.22A

Thanks,

Jeff

#pragma rtGlobals=1		// Use modern global access method.
#pragma IndependentModule= AFM_PackageLoader
 
#ifndef AFM_PackageLoader
	Menu "jtgSpectroscopy"
		SubMenu "Packages..."  
			"Load AFM Package", /Q, LoadAFMPackage()
			help = {"Loads the set of procedures related to Atomic Force Microscopy."}
		End	
	End
#else
	Menu "jtgSpectroscopy"
		SubMenu "Packages..."
			"Unload AFM Package", /Q, UnloadAFMPackage()
			help = {"Removes the set of procedures related to Atomic Force Microscopy."}
		End
	End
#endif
 
//package loader
//"SetIgorOption poundDefine=AFM_PackageLoader" defines a global symbol
//that the force curve analysis package uses to add a "remove package" menu item
Function LoadAFMPackage()
	Execute/P/Q/Z "INSERTINCLUDE \"AFM_ForceCurveAnalysis\""
	Execute/P/Q/Z "INSERTINCLUDE \"AFM_MultiFileLoader\""
	Execute "SetIgorOption poundDefine=AFM_PackageLoader"
	Execute/P/Q/Z "COMPILEPROCEDURES "	
End
 
//package unloader
//"SetIgorOption poundunDefine=AFM_PackageLoader" clears a global symbol
//that the force curve analysis package uses to add a remove package menu item
Function UnloadAFMPackage()
	Execute/P/Q/Z "DELETEINCLUDE \"AFM_ForceCurveAnalysis\""
	Execute/P/Q/Z "DELETEINCLUDE \"AFM_MultiFileLoader\""
	Execute "SetIgorOption poundunDefine=AFM_PackageLoader"
	Execute/P/Q/Z "COMPILEPROCEDURES "
End

Igor in Linux using WINE

RGerkin — Thu, 25 Jun 2009 09:53:09 -0700

While I continue to use Windows on computers that control instruments, I now use Linux on the machine on which I do analysis (and everything else). For a long time, a major barrier to long-term Linux use was Igor compatibility, namely that there are only official versions of Igor for Mac and Windows. However, since Igor appears to use the Win32 API almost exclusively for the Windows edition, it is possible to run Igor with the help of a free program in Linux called "WINE", which implements the Win32 API fairly well. To the user, everything is transparent; setup of Igor proceeds as in Windows, and you can double-click to run Igor just like you were running it in Windows. There appears to be little of no loss of speed for most tasks in Igor.

There are, however, a small number of discrepancies between running Igor for Windows in Windows, and running Igor for Windows in Linux using WINE. Most of these are small, and don't stop me from using it. But some of these may be significant enough to cause problems for some users. The purpose of this thread is to document the inconsistencies (not really bugs, since Linux was not intended by Wavemetrics), the circumstances under which they occur, and any possible workarounds. Some of these may be fixable by changing one's OS or desktop environment settings, others may be fixable by modifying the WINE source code (an open source project), and others may be fixable by Wavemetrics (but don't count on that -- I'm sure their priority is Windows/Mac).

So I propose that posts be formatted something like:
-- Name and Version of Linux Distribution
-- (Optional): Graphics information (desktop environment, e.g. KDE vs Gnome, graphics card/drivers, X11 version)
-- WINE version
-- Igor version
-- Summary of problem

I will follow up later with posts describing some of the discrepancies that I have observed.

Rick

Methods for generating/storing a random number of variables?

daggaz — Thu, 10 May 2012 01:26:58 -0700

I have a data set(wave), of random length. I want to generate variables (lets say anywhere from 3 to 10 of them) based on this length, for logical operations on my data. I realized pretty quickly that I couldn't use a for loop and create variable_i as i iterates thru the loop.

What I am doing now is just creating a new wave of the appropriate length and calling its point values as my variables, instead of creating actual variables. But so far I am forced to tack it onto a table and fill in the values, which seems cluttered (id prefer to have it run in the background, with GUI for input, instead of actually "seeing" the wave. )

I've tried running a set_variable GUI in a loop so that I fill in wave(i) points, but the program just hauls ass thru the loop, showing the GUI only once and I end up with a wave full of zeros and only the last point contains any user input. I was hoping more for a screen showing "i" set_variable GUIs so that I could set them all one by one. Putting a button in at the bottom so that the next function is called only after the user has input their values would be the next, and easy, step.

Or is there a more elegant solution?

small multithread errors

ubi — Tue, 08 May 2012 00:27:54 -0700

To treat my data I have a matlab script written by a colleague which I translated into igor and adapted to my needs. The procedure is quite fit extensive. The result of both programs are in pretty good agreement (without surprise). To enhance the speed I use a multithreaded code. However the result obtained with igor multithreaded code is a little bit more noisy (for the same conditions as far as it is possible to manage them) as it is shown in the attached file while the "normal" gives exactly the same result as the matlab script.

To be a bit more precise: Each x point is the result of the the difference of two different datasets proceeded by the following function.

Threadsafe FUNCTION fitthematrix_realth(srcwv, tmpwv,outwv,row)
WAVE srcwv, tmpwv, outwv
variable row
variable i
		tmpwv = srcwv[row][p]
		CurveFit/N/Q/L=100000 /NTHR=1 poly 4,  tmpwv /D
		wave fit_dumper, M_Covar, W_coef, W_fitConstants, W_ParamConfidenceInterval, W_sigma
		wavestats/Q/M=1 fit_dumper
		outwv[row] = 0.0009 * V_minRowLoc
		killwaves/Z fit_dumper, M_Covar, W_coef, W_fitConstants, W_ParamConfidenceInterval, W_sigma
END

In the main function the loop which calls the fit procedure looks like, which is basically copy&pasted from the igor help files:

Variable n,actrow,nthreads= ThreadProcessorCount
Variable mt= ThreadGroupCreate(nthreads)
 
Variable ttime= stopMSTimer(-2)
for(actrow=0;actrow<DimSize(fieldwv,0);)
	for(n=0;n<nthreads;n+=1)
		ThreadStart mt,n, fitthematrix_realth(fieldwv, dumper,resultingfield,actrow)
		actrow+=1
		if( actrow >= DimSize(fieldwv,0) )
			break
		endif
	endfor
 
	do
		variable tgs= ThreadGroupWait(mt,1000)
	while( tgs != 0 )
endfor
variable dummy= ThreadGroupRelease(mt)
 
Variable mt2= ThreadGroupCreate(nthreads)
for(actrow=0;actrow<DimSize(zfieldwave,0);)
	for(n=0;n<nthreads;n+=1)
		ThreadStart mt2,n, fitthematrix_realth(zfieldwave, dumper,resultingnofield,actrow)
		actrow+=1
		if( actrow >= DimSize(zfieldwave,0) )
			break
		endif
	endfor
 
	do
		variable tgs2= ThreadGroupWait(mt2,1000)
	while( tgs != 0 )
endfor					
dummy= ThreadGroupRelease(mt2)

Where is my error?

plotting data using bars

CaJo — Mon, 07 May 2012 09:07:28 -0700

Hi,

I'm sure it's an easy question but right now I'm stuck with the following problem and can't find out how to solve it:

I have xy-data-sets (a few hundred data points), whose x-values are not evenly distributed, e.g.:
1 y1
3 y2
8 y3

How can I plot this data set so that it finally looks like the blue bars in http://www.wavemetrics.com/products/igorpro/gallery/userpix/friedman.gif ?

If I use an xy-plot with Modify Graph mode = 5 I end up with bars, but they are connected to each other which is not what I want. I have not found out yet how to change it (if possible). Also, it's not really a category plot.

I already tried to search the documentation, but end up with a lot of information on error bars and category plots. Any help is appreciated!

Is anyone here using the "NIST SANS data analysis program" written in Igor?

Igor_user — Thu, 09 Feb 2012 17:49:32 -0800

Dear All,

I'm currently using the above mentioned program for analysing some data I acquired from small-angle scattering experiments.

It's a free program downloadable from NIST website and fully compatible with Igor.

So, I'd like to ask if here is any small-angle scattering experts who is familiar with this program.

Best Wishes,

Using Hooks in Independent Modules

jjweimer — Tue, 01 Apr 2008 14:38:46 -0700

I have a question about using a Menu Hook Function in an Independent Modules. The following options do not work for me:

#pragma IndependentModule=MyModule
 
SetIgorHook IgorMenuHook=TestTheMenu    // does not work with this line
SetIgorHook IgorMenuHook=MyModule#TestTheMenu    // does not work with this line
SetIgorHook IgorMenuHook=MyModule#MyModule#TestTheMenu    // does not work with this line
 
Static Function TestTheMenu(...)
   ...
 
end

Neither does the following direct setup:

#pragma IndependentModule=MyModule
 
Function IgorMenuHook(...)
   ...
 
end

Is this a limitation of Independent Modules, a my misunderstanding of the text on page IV-41 of the manual ...

Quote:

However, if need be, you can call static functions in the independent module from outside the group using a triple name syntax:

imName#modName#functionName()

... or a limitation of Independent Modules AND IgorMenuHook functions?

BTW, I am still generally confused about the first two names in the triple name sequence imName#modName#functionName. What is the modName of an Independent Module?