Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1217 1473 (N32W09 L=287)   342    102        0 (1202/1422-1232/1521) -74.1 W07-W09 N28-N35
      2 1404 1263 (N20W23 L=301)   317    105        1 (1365/1236-1458/1281)  -8.1 W19-W26 N18-N20
      3 1797 1092 (N12W56 L=334)   851    210        1 (1754/1004-1828/1188) -54.0 W51-W59 N07-N18
      4  945 1122 (N06E06 L=273)   205     59        1 ( 919/1095- 968/1153)  63.4 E07-E04 N04-N08
      5  189 1008 (S09E63 L=215)   624    101        1 ( 138/ 987- 230/1023)   0.0 E71-E57 S09-S07
      6  149  922 (S14E70 L=208)  1016    105        1 ( 111/ 878- 190/ 954) -13.0 E79-E63 S16-S12
      7  453  862 (S15E37 L=242)   429     72        2 ( 414/ 831- 477/ 902)  78.7 E39-E34 S16-S11
      8 1018  856 (S09W01 L=280)   318     84        1 ( 988/ 815-1046/ 896)  78.7 W00-W02 S11-S06
      9  381  844 (S16E43 L=236)   427    137        2 ( 321/ 799- 421/ 893) -24.0 E48-E39 S19-S14
     10  302  781 (S22E51 L=227)  1244    214        3 ( 261/ 720- 367/ 867) -41.6 E54-E44 S25-S16
     11  511  808 (S17E32 L=246)   260     69        0 ( 502/ 767- 525/ 846)  78.7 E32-E30 S19-S14
     12 1700  576 (S20W58 L=336)   416     64        0 (1678/ 548-1722/ 606)  63.4 W56-W59 S22-S18
  

    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