Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  627 1683 (N23E47 L=113)  2000    356        1 ( 595/1563- 671/1808) -56.3 E53-E40 N15-N33
      2  389 1603 (N13E62 L= 98)  1492    242        1 ( 367/1512- 421/1684) -59.7 E67-E57 N06-N18
      3  901 1505 (N19E21 L=139)   224     71        0 ( 885/1470- 923/1542) -78.7 E21-E19 N16-N21
      4  410 1405 (N01E50 L=111)   434     79        2 ( 387/1366- 434/1449) -63.4 E51-E47 N00-N04
      5  258 1370 (S04E63 L= 97)   924    128        2 ( 229/1332- 289/1411)  78.7 E67-E59 S05-N00
      6  778 1360 (N07E23 L=137)   436    104        0 ( 727/1330- 826/1398)  56.3 E25-E21 N04-N10
      7 1176 1279 (N13W02 L=162)   582    152        1 (1120/1212-1225/1354) -32.0 E03-W05 N09-N16
      8  892 1297 (N07E15 L=146)   282    119        0 ( 862/1241- 941/1325)  71.6 E16-E12 N03-N09
      9 1776 1136 (N24W46 L=206)   653     81        0 (1755/1093-1798/1180) -68.2 W43-W47 N20-N25
     10  174  988 (S28E60 L=101)  1373    373        1 ( 115/ 889- 241/1069)  47.3 E65-E53 S34-S21
     11  636  830 (S26E18 L=142)  1038    270        1 ( 519/ 811- 759/ 851)  26.6 E26-E10 S29-S21
     12 1181  762 (S15W16 L=176)   511    225        0 (1081/ 720-1298/ 782)  36.0 W11-W22 S19-S11
     13 1859  720 (N04W70 L=230)   893    118        1 (1839/ 678-1871/ 768) -59.0 W65-W74 N01-N06
     14 1299  731 (S13W24 L=184)   228     62        2 (1266/ 714-1336/ 748)  36.9 W22-W26 S14-S11
     15 1205  659 (S20W21 L=181)   869    253        0 (1097/ 637-1339/ 682)  23.2 W14-W28 S22-S14
     16 1635  542 (S13W55 L=215)  1308    331        8 (1545/ 469-1697/ 631)  72.9 W51-W60 S19-S06
     17 1543  530 (S17W48 L=209)   524    102        3 (1502/ 503-1584/ 560)  45.0 W45-W51 S18-S13
  

    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