Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  658 1743 (N27E48 L=115)   544     84        0 ( 637/1702- 685/1784) -60.3 E50-E44 N23-N30
      2 1014 1699 (N33E21 L=141)   476     76        1 ( 985/1669-1051/1729) -15.9 E24-E17 N31-N34
      3  378 1519 (N07E56 L=106)  1504    288        2 ( 270/1481- 498/1556)  24.0 E66-E47 N02-N10
      4  786 1374 (N06E22 L=140)   263    103        1 ( 753/1336- 817/1423)  90.0 E23-E21 N03-N09
      5  351 1340 (S04E51 L=111)   707    196        0 ( 286/1270- 403/1403) -18.4 E59-E44 S07-S01
      6 1617 1357 (N29W28 L=191)   419    106        0 (1582/1303-1646/1395) -33.7 W24-W31 N26-N30
      7  225 1270 (S10E61 L=101)   616    102        1 ( 209/1213- 254/1317)  90.0 E63-E59 S13-S07
      8 1386 1214 (N13W15 L=178)   201     80        0 (1363/1188-1423/1232)  45.0 W13-W17 N11-N14
      9  362 1172 (S14E45 L=117)   558    118        2 ( 324/1108- 402/1212)  66.8 E47-E42 S18-S11
     10  447 1149 (S13E38 L=124)   441    123        0 ( 417/1107- 504/1179)  45.0 E40-E34 S16-S11
     11  104 1119 (S20E76 L= 86)   758     64        1 (  84/1090- 129/1149)  18.4 E81-E71 S21-S18
     12  573  811 (S30E23 L=140)   214     47        1 ( 543/ 797- 600/ 834)  14.0 E24-E20 S31-S28
     13 1364  757 (S12W27 L=189)   938    319        1 (1220/ 725-1502/ 793)  32.0 W18-W34 S17-S07
  

    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