Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  847 1846 (N39E46 L=207)   413     79        1 ( 834/1805- 858/1882) -41.2 E51-E42 N35-N42
      2  704 1821 (N33E55 L=199)   572    123        1 ( 683/1766- 733/1867) -39.3 E60-E48 N28-N37
      3 1602 1669 (N56W25 L=279)   813     68        1 (1578/1634-1622/1711)  63.4 W21-W29 N51-N61
      4  664 1557 (N14E38 L=216)   560    208        1 ( 611/1444- 706/1614) -33.7 E43-E30 N08-N17
      5  495 1261 (S07E38 L=215)   726    140        0 ( 433/1226- 555/1289)  40.6 E42-E33 S10-S04
      6 1599 1209 (N21W29 L=282)   854    231        4 (1510/1126-1674/1272)  -8.1 W21-W35 N17-N23
      7  269 1084 (S22E53 L=200)  1077    318        0 ( 190/ 953- 360/1178)  48.6 E61-E46 S31-S14
      8  174 1105 (S22E65 L=188)   499     84        1 ( 147/1070- 197/1144)  45.0 E69-E62 S24-S19
      9  409  983 (S25E38 L=215)   284     67        1 ( 383/ 958- 445/1007)  45.0 E40-E35 S26-S22
     10  603  813 (S30E20 L=233)  1295    382        2 ( 455/ 748- 778/ 863)   9.1 E31-E06 S32-S27
     11 1325  472 (S27W38 L=291)   313     95        0 (1279/ 441-1361/ 504)  71.6 W36-W39 S30-S24
  

    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