Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1404 1780 (N55W48 L=132)   426     44        0 (1385/1756-1430/1797)  14.0 W42-W54 N53-N56
      2 1231 1773 (N56W25 L=110)   750    147        0 (1171/1738-1301/1795)   9.5 W17-W35 N52-N57
      3 1031 1450 (N31W01 L= 86)  1018    309        1 ( 931/1367-1150/1526) -35.5 E05-W09 N24-N36
      4 1725 1475 (N29W60 L=145)  1898    221        3 (1638/1439-1822/1510)   4.6 W48-W85 N27-N30
      5  277 1349 (N23E59 L= 26)   872    104        2 ( 249/1310- 313/1392) -59.0 E63-E55 N21-N26
      6 1127 1281 (N19W07 L= 92)   257    130        0 (1065/1270-1187/1298)   7.1 W03-W11 N18-N20
      7 1424 1140 (N10W26 L=110)   246    119        0 (1398/1087-1466/1193)  54.5 W23-W29 N06-N13
      8  342 1033 (N04E47 L= 38)   490    121        2 ( 289/1010- 386/1068)  24.0 E51-E42 N02-N06
      9  744  939 (S01E17 L= 67)   233     47        0 ( 731/ 910- 758/ 977)  90.0 E18-E16 S02-N01
     10  692  787 (S11E21 L= 63)   275     98        0 ( 640/ 765- 742/ 804)  15.9 E25-E18 S11-S09
     11 1771  711 (S19W56 L=140)   418     53        0 (1745/ 685-1799/ 730)   9.5 W53-W59 S20-S17
     12 1394  413 (S38W28 L=113)   542     81        0 (1361/ 386-1424/ 447) -40.6 W24-W31 S40-S34
  

    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