Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1103 1777 (N54E08 L=163)   852    168        0 (1015/1751-1178/1801)  26.6 E15-W00 N50-N58
      2 1339 1750 (N57W18 L=189)   422     57        1 (1318/1720-1364/1770)  71.6 W15-W21 N53-N59
      3 1438 1724 (N56W30 L=201)   650     70        2 (1410/1698-1473/1748)   6.3 W25-W35 N53-N58
      4  718 1567 (N30E31 L=140)   365     78        2 ( 698/1535- 741/1603) -68.2 E32-E28 N28-N33
      5  316 1484 (N19E64 L=108)  1754    209        7 ( 274/1404- 364/1581) -56.3 E73-E58 N13-N25
      6 1383 1456 (N33W19 L=191)   292     60        0 (1351/1441-1417/1471)   0.0 W16-W22 N31-N33
      7 1778 1345 (N31W57 L=228)   550     67        0 (1742/1327-1813/1364)  15.3 W51-W62 N29-N32
      8  864 1325 (N17E14 L=157)   212     74        0 ( 830/1307- 904/1340)   0.0 E16-E11 N16-N17
      9  678  975 (S07E21 L=150)   345     66        2 ( 649/ 937- 707/1006) -53.1 E22-E18 S09-S04
  

    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