Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1774 1446 (N45W46 L=138)  1037    115        0 (1719/1424-1833/1463)  11.3 W37-W57 N41-N47
      2  538 1207 (S09E34 L= 58)   237     36        0 ( 520/1183- 554/1230)  90.0 E34-E32 S10-S07
      3  369 1172 (S15E46 L= 46)   268     90        0 ( 326/1157- 415/1184)  23.2 E49-E41 S16-S13
      4 1730 1087 (N19W42 L=134)   542    152        1 (1694/1015-1782/1133) -16.7 W38-W48 N16-N21
      5  579 1007 (S19E26 L= 66)   537    172        0 ( 520/ 956- 645/1063)  56.3 E29-E22 S23-S14
      6  417 1000 (S24E38 L= 54)   470    185        0 ( 335/ 963- 483/1021)  13.0 E44-E31 S26-S21
      7  291  762 (S41E49 L= 43)  1182    151        0 ( 248/ 701- 338/ 834)  40.2 E56-E43 S46-S35
      8 1623  759 (S03W44 L=136)   268     91        1 (1576/ 743-1666/ 772)  23.2 W40-W47 S04-N00
      9 1308  539 (S24W33 L=125)  1004    199        2 (1229/ 509-1386/ 575)  53.1 W28-W37 S28-S20
     10 1137  482 (S33W26 L=118)   268     65        0 (1112/ 452-1159/ 512) -78.7 W24-W27 S35-S30
     11 1006  348 (S46W26 L=118)  1318    228        3 ( 971/ 290-1083/ 458) -21.3 W17-W35 S50-S39
     12 1171  391 (S37W34 L=126)   523     71        0 (1126/ 361-1202/ 424)   0.0 W28-W37 S38-S35
  

    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