Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1632 (N40E13 L=179)   709    180        0 ( 918/1598-1045/1689) -18.4 E18-E06 N38-N44
      2 1708 1480 (N38W54 L=247)  2539    326        4 (1645/1350-1747/1614)  69.8 W48-W64 N29-N48
      3  352 1460 (N20E58 L=134)   399     31        1 ( 335/1438- 366/1486) -56.3 E59-E55 N18-N21
      4 1108 1260 (N16W03 L=195)   486    235        0 (1059/1216-1152/1329) -56.3 W00-W05 N14-N20
      5 1441 1182 (N15W25 L=217)   261    106        0 (1416/1144-1478/1225) -45.0 W23-W28 N13-N17
      6  806 1075 (N02E14 L=178)   360     94        2 ( 785/1025- 827/1116) -78.7 E15-E12 S01-N04
      7  502  888 (S13E33 L=159)   352     39        1 ( 476/ 864- 526/ 913)  45.0 E34-E30 S14-S11
      8 1679  872 (S01W46 L=238)   588    122        1 (1633/ 832-1720/ 913) -18.4 W41-W50 S03-N00
      9 1216  806 (S10W14 L=207)   216     66        1 (1182/ 780-1241/ 829)  53.1 W12-W15 S12-S08
     10  821  799 (S15E10 L=182)   242     76        0 ( 779/ 787- 857/ 811)  -9.5 E13-E07 S15-S14
     11  453  783 (S21E37 L=156)   365    102        2 ( 412/ 762- 506/ 810)   0.0 E40-E32 S22-S19
     12  708  652 (S26E17 L=175)   226     23        0 ( 692/ 631- 723/ 677) -71.6 E18-E15 S27-S24
     13  874  650 (S24E06 L=187)   304     58        1 ( 841/ 628- 904/ 672)  26.6 E07-E03 S25-S22
     14 1372  550 (S25W30 L=222)   823    103        2 (1320/ 523-1423/ 579)   7.1 W26-W34 S26-S22
     15  891  498 (S35E03 L=189)   340     91        0 ( 845/ 488- 940/ 512)   0.0 E06-W00 S35-S33
     16  455  462 (S44E44 L=149)  2144    215        5 ( 369/ 427- 540/ 508)  -2.9 E54-E33 S46-S40
  

    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