Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1363 1783 (N56E04 L= 51)   407     40        0 (1342/1760-1383/1802) -80.5 E06-E01 N52-N58
      2  768 1649 (N23E36 L= 19)   724    163        1 ( 693/1619- 841/1677)  10.3 E41-E30 N21-N25
      3 1409 1592 (N40W07 L= 62)   345     34        2 (1385/1572-1435/1617)  36.9 W05-W09 N38-N42
      4 1613 1369 (N32W28 L= 83)   425    109        0 (1568/1324-1649/1402) -14.0 W23-W31 N30-N33
      5  482 1074 (S17E34 L= 21)   629    219        1 ( 394/1026- 558/1145)  45.0 E40-E30 S21-S11
      6 1800 1012 (N17W51 L=106)   433     85        0 (1779/ 981-1818/1053)  78.7 W49-W52 N14-N19
      7 1719  965 (N12W45 L= 99)   313     70        0 (1697/ 921-1737/ 999)  90.0 W43-W46 N08-N13
      8  308  930 (S29E47 L=  8)   395     94        1 ( 266/ 893- 346/ 963)  40.6 E51-E43 S32-S26
      9 1530  543 (S16W47 L=102)   379     77        0 (1500/ 509-1561/ 575)  90.0 W46-W48 S19-S13
     10  650  403 (S56E06 L= 48)   288     33        0 ( 633/ 377- 665/ 427) -36.9 E08-E04 S57-S53
     11 1420  340 (S31W54 L=109)   523     69        0 (1398/ 316-1450/ 372) -14.0 W50-W58 S32-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