Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1646 (N23E62 L= 90)  1156    240        1 ( 392/1587- 432/1699) -41.6 E68-E58 N18-N26
      2  458 1439 (N11E47 L=105)   335    108        0 ( 430/1395- 482/1474) -29.7 E50-E43 N08-N12
      3  331 1271 (S02E52 L=100)   534    131        1 ( 308/1217- 365/1325) -39.8 E55-E47 S04-N01
      4 1573 1203 (N20W31 L=183)   636    122        0 (1533/1139-1602/1262) -49.4 W27-W33 N16-N23
      5 1343 1190 (N14W15 L=167)   210     88        0 (1305/1149-1368/1240) -38.7 W12-W17 N12-N16
      6 1807 1207 (N27W52 L=203)   367     41        0 (1787/1185-1828/1227) -33.7 W49-W53 N25-N28
      7  594  904 (S18E24 L=128)   235     65        0 ( 573/ 873- 625/ 928)  63.4 E25-E22 S20-S16
      8  358  875 (S25E42 L=110)  1098    331        1 ( 276/ 799- 470/ 924)  30.5 E49-E32 S31-S20
      9 1117  832 (S11W09 L=161)   303     91        0 (1071/ 814-1172/ 852)  38.7 W06-W12 S13-S09
     10  886  713 (S23E02 L=150)   482    196        0 ( 792/ 690- 953/ 775)  41.2 E07-W01 S25-S18
     11  959  486 (S36W09 L=161)   481    149        0 ( 884/ 471-1040/ 503)  26.6 W03-W14 S38-S33
     12 1414  446 (S28W42 L=194)   909    212        0 (1330/ 417-1480/ 495)  49.4 W35-W45 S31-S23
  

    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