Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  373 1527 (N13E60 L=292)   651    110        1 ( 355/1470- 394/1574) -54.5 E63-E56 N09-N16
      2  866 1547 (N25E23 L=328)   710    230        1 ( 783/1504- 933/1573)   5.7 E28-E17 N21-N27
      3 1495 1386 (N30W22 L= 14)   698    218        2 (1454/1286-1533/1478) -61.4 W18-W26 N24-N35
      4  828 1262 (N07E17 L=335)   706    233        1 ( 762/1218- 891/1338) -16.7 E22-E12 N04-N10
      5  275  996 (S20E52 L=300)   453     84        1 ( 232/ 971- 314/1025)  24.0 E56-E47 S22-S17
      6 1191  635 (S21W20 L= 11)   261     70        2 (1155/ 617-1228/ 652)  36.9 W17-W21 S22-S19
      7 1518  597 (S16W43 L= 35)   277     61        1 (1483/ 582-1553/ 609)  21.8 W40-W45 S16-S14
      8 1446  321 (S34W56 L= 48)   631     78        1 (1422/ 297-1479/ 362) -16.7 W51-W61 S35-S31
  

    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