Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  912 1488 (N17E20 L=346)   227     58        0 ( 883/1470- 943/1506)  14.0 E21-E17 N15-N18
      2 1783 1396 (N41W44 L= 50)   479    110        1 (1736/1378-1810/1413)  18.4 W38-W47 N39-N43
      3 1631 1328 (N30W29 L= 36)  1001    245        2 (1550/1226-1715/1409) -21.8 W22-W37 N25-N33
      4  915 1346 (N09E15 L=351)   531     65        2 ( 875/1307- 947/1380)  78.7 E17-E13 N06-N11
      5  935 1265 (N05E11 L=355)   351     73        2 ( 904/1239- 961/1296) -33.7 E13-E09 N03-N07
      6  301  976 (S27E46 L=320)   633    210        2 ( 247/ 919- 371/1030) -13.0 E52-E39 S29-S24
      7 1860  882 (N12W60 L= 67)   470     36        0 (1844/ 858-1881/ 903)  71.6 W58-W62 N10-N13
      8 1114  784 (S16W12 L= 18)   705    178        2 (1055/ 739-1194/ 832)  -5.7 W07-W17 S18-S13
      9 1518  272 (S30W64 L= 70)   636    202        0 (1491/ 253-1551/ 299)   0.0 W60-W68 S31-S29
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y