Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1005 1873 (N60W32 L=225)  1783    246        3 ( 932/1806-1113/1935)   8.7 W18-W57 N55-N62
      2 1220 1741 (N42W36 L=228)   309     45        0 (1195/1715-1247/1763)  15.9 W32-W39 N40-N42
      3 1325 1728 (N38W43 L=235)   390     57        0 (1289/1715-1362/1740) -31.0 W39-W46 N36-N39
      4  247 1162 (N25E53 L=139)   788    160        0 ( 232/1085- 269/1247) -73.3 E54-E50 N20-N30
      5 1668 1016 (S11W40 L=232)   308     76        0 (1646/ 989-1697/1051) -63.4 W38-W42 S13-S09
      6 1878  720 (S35W65 L=258)  1065     87        2 (1848/ 690-1913/ 767) -39.8 W60-W77 S37-S32
      7 1126  667 (S21E01 L=191)   561    259        2 (1027/ 596-1232/ 733)  12.1 E08-W06 S23-S17
      8  240  657 (S05E65 L=127)   505     52        0 ( 221/ 628- 261/ 682)  26.6 E67-E60 S06-S03
      9  876  353 (S36E28 L=165)  1247    337        1 ( 723/ 335-1039/ 382) -20.0 E37-E15 S40-S30
  

    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