Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1030 1633 (N36W03 L=324)   413     57        1 (1011/1594-1045/1671) -90.0 W01-W04 N33-N39
      2 1762 1374 (N17W54 L= 15)   371     64        0 (1735/1353-1789/1394) -36.9 W51-W56 N15-N18
      3 1482 1199 (N06W29 L=350)  1308    314        3 (1349/1062-1564/1291) -45.0 W21-W34 S01-N12
      4 1819 1165 (N04W56 L= 17)   526     91        2 (1786/1137-1844/1195) -90.0 W52-W59 N02-N05
      5  531  934 (S06E31 L=290)   342     95        0 ( 490/ 900- 575/ 966)  23.2 E34-E27 S07-S03
      6  413  755 (S16E42 L=279)  1732    408        4 ( 347/ 615- 519/ 859)  35.8 E50-E32 S23-S09
      7 1353  772 (S19W20 L=341)   260     60        0 (1334/ 737-1372/ 805)  76.0 W18-W21 S20-S16
      8 1510  400 (S45W39 L=359)   615    164        1 (1449/ 382-1582/ 418)   0.0 W32-W45 S46-S43
  

    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