Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1226 1752 (N45E12 L= 14)  2610    421        5 (1104/1600-1377/1887) -22.4 E31-W04 N37-N54
      2  601 1742 (N25E54 L=332)   933    121        1 ( 551/1710- 655/1774)  50.2 E57-E50 N21-N27
      3  990 1695 (N32E24 L=  2)   549    161        1 ( 927/1635-1017/1757) -28.6 E31-E20 N28-N34
      4  612 1675 (N20E48 L=339)  1573    293        4 ( 505/1600- 722/1747)  65.0 E51-E43 N13-N28
      5  854 1635 (N24E29 L=357)   289     84        0 ( 840/1590- 869/1675) -56.3 E31-E27 N20-N26
      6 1315 1433 (N24W05 L= 32)   409    168        0 (1279/1384-1363/1515) -33.7 W00-W09 N21-N28
      7  178 1026 (S26E64 L=323)   657    125        2 ( 142/ 981- 214/1088)  35.0 E69-E59 S29-S22
      8  619  952 (S21E22 L=  4)   388     59        0 ( 572/ 935- 653/ 971)  18.4 E25-E19 S23-S19
      9 1299  521 (S26W33 L= 60)   239     80        0 (1266/ 500-1336/ 542)  63.4 W32-W34 S28-S24
     10 1302  380 (S34W43 L= 69)   301     75        1 (1267/ 363-1343/ 404)  53.1 W40-W44 S35-S31
     11 1302  202 (S44W64 L= 90)   625     74        1 (1269/ 171-1335/ 233)   3.8 W57-W73 S44-S42
  

    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