Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1236 1399 (N27W15 L=196)   566    160        2 (1206/1319-1265/1461)  68.2 W13-W18 N21-N31
      2  415 1322 (N23E43 L=137)   401     71        1 ( 397/1281- 434/1359) -78.7 E45-E42 N20-N25
      3 1621 1284 (N18W42 L=223)  1129    162        3 (1562/1236-1686/1333) -41.2 W37-W47 N14-N21
      4 1796  840 (S11W56 L=237)  2137    333        5 (1732/ 757-1868/ 923)  34.7 W49-W64 S15-S05
      5 1703  785 (S14W47 L=228)   396     70        1 (1676/ 755-1721/ 820)  53.1 W45-W49 S15-S11
      6 1849  718 (S19W65 L=246)  1087    152        2 (1815/ 682-1882/ 751)  12.5 W60-W69 S21-S17
      7 1452  669 (S20W28 L=209)   482    163        0 (1416/ 607-1482/ 736) -76.0 W25-W30 S23-S15
      8  544  682 (S17E33 L=148)   401    110        1 ( 511/ 651- 592/ 715) -26.6 E35-E29 S19-S15
      9 1691  638 (S23W49 L=229)   881    118        3 (1633/ 599-1741/ 682)  24.0 W43-W52 S25-S20
     10  631  438 (S35E32 L=149)   371     79        0 ( 609/ 415- 658/ 459)  31.0 E34-E29 S36-S33
     11 1639  421 (S40W52 L=233)   903    159        2 (1585/ 406-1678/ 445)  -5.7 W46-W57 S41-S37
  

    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