Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  722 1863 (N37E54 L= 71)   568     76        0 ( 698/1845- 743/1880)   0.0 E56-E51 N35-N39
      2  765 1797 (N33E45 L= 81)  1082    189        3 ( 733/1733- 801/1847) -45.0 E50-E40 N27-N36
      3  389 1617 (N12E63 L= 63)  1221    180        1 ( 344/1580- 468/1646)  29.1 E67-E55 N09-N14
      4 1330 1322 (N19W09 L=134)   699    193        1 (1297/1238-1386/1418) -38.7 W04-W14 N15-N23
      5  453 1003 (S22E34 L= 91)   353    101        2 ( 402/ 979- 497/1029)   7.1 E38-E30 S24-S20
      6 1313  325 (S36W46 L=171)   411     87        1 (1278/ 297-1349/ 348)   8.1 W42-W49 S36-S34
      7 1507  295 (S30W61 L=187)   936    255        1 (1478/ 248-1539/ 342) -16.7 W57-W67 S30-S27
  

    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