Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1111 1864 (N55E09 L=342)   628     96        1 (1078/1836-1163/1882)  35.5 E12-E03 N51-N58
      2  454 1708 (N33E60 L=291)   542     57        0 ( 443/1670- 471/1741) -36.9 E64-E56 N29-N35
      3  434 1632 (N26E55 L=296)   436     68        2 ( 410/1600- 458/1658) -20.6 E58-E50 N24-N28
      4 1180 1508 (N26W04 L=355)   480    153        0 (1147/1442-1202/1577) -77.5 W00-W05 N21-N30
      5  888 1544 (N25E16 L=335)   217     95        0 ( 851/1514- 937/1558)  -9.5 E18-E12 N23-N25
      6 1786 1252 (N20W51 L= 42)  1036    234        1 (1749/1144-1814/1355) -76.0 W48-W53 N13-N25
      7  395  951 (S16E41 L=310)   324     64        0 ( 379/ 915- 419/ 982)  63.4 E42-E38 S18-S14
      8  623  787 (S24E23 L=328)   527    184        0 ( 587/ 718- 645/ 887)  79.7 E25-E21 S29-S18
      9  448  841 (S23E37 L=314)   295     97        1 ( 408/ 814- 465/ 878)  45.0 E39-E35 S24-S20
     10 1117  205 (S60W32 L= 23)   667     78        2 (1093/ 174-1149/ 242) -28.6 W27-W38 S62-S56
  

    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