Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  502 1265 (N24E30 L=208)   441    127        0 ( 456/1217- 540/1306) -49.4 E33-E27 N20-N27
      2  175 1212 (N27E63 L=174)   569     42        1 ( 150/1188- 203/1231)   7.1 E67-E58 N25-N28
      3  226 1140 (N21E55 L=182)   351     50        1 ( 197/1123- 256/1160)   0.0 E58-E51 N20-N22
      4  232 1075 (N17E55 L=183)   324     41        2 ( 212/1050- 251/1093) -45.0 E56-E52 N15-N18
      5 1897  961 (S18W64 L=301)   491     46        1 (1875/ 932-1916/ 985)  45.0 W61-W67 S20-S16
      6 1757  702 (S32W48 L=286)  1924    258        1 (1659/ 627-1834/ 802)  22.4 W38-W56 S36-S26
      7 1502  694 (S27W25 L=262)   309     94        0 (1470/ 653-1532/ 737)  45.0 W22-W27 S29-S24
      8  925  468 (S30E18 L=219)   379    119        0 ( 878/ 446- 987/ 487)   0.0 E21-E13 S31-S29
      9 1619  365 (S56W45 L=283)   886     99        0 (1590/ 345-1661/ 409) -24.0 W40-W58 S59-S51
     10  595  373 (S31E46 L=191)   412     55        1 ( 568/ 349- 623/ 399) -56.3 E48-E44 S32-S28
  

    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