Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1120 1886 (N54E34 L=293)   822     83        1 (1087/1863-1159/1911)  71.6 E39-E28 N50-N57
      2  738 1827 (N35E53 L=274)   631     47        2 ( 721/1801- 762/1856) -35.5 E57-E49 N32-N37
      3  703 1524 (N13E34 L=293)  1352    326        1 ( 671/1379- 737/1674) -58.6 E41-E29 N04-N22
      4 1216 1539 (N29E04 L=322)   424    133        0 (1164/1490-1248/1588) -24.0 E09-W00 N26-N30
      5  156 1051 (S26E68 L=259)   740    116        1 ( 138/ 996- 171/1100) -71.6 E70-E64 S28-S22
      6  937  709 (S26W04 L=331)   776    326        0 ( 885/ 593- 993/ 828) -63.4 W02-W10 S31-S17
      7  590  722 (S36E19 L=307)  1786    348        4 ( 421/ 675- 749/ 772)  31.0 E33-E08 S43-S28
      8 1266  697 (S17W24 L=350)   471    108        2 (1219/ 667-1320/ 734)  38.7 W20-W26 S18-S14
      9 1089  602 (S27W17 L=344)   944    184        2 ( 988/ 580-1173/ 635)  31.0 W11-W22 S30-S24
     10 1403  542 (S21W39 L=  5)  1213    218        3 (1322/ 495-1482/ 597)  63.4 W34-W42 S25-S15
     11  616  354 (S62E08 L=319)  2244    299        7 ( 507/ 312- 747/ 411)  26.6 E25-W03 S69-S51
     12 1496  375 (S26W58 L= 24)   819    201        2 (1393/ 358-1583/ 395)  29.7 W50-W65 S30-S22
  

    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