Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1531 (N38E08 L=302)   513    157        0 ( 807/1492- 952/1559)   9.5 E13-E01 N35-N39
      2  554 1529 (N39E36 L=274)   401     88        2 ( 510/1503- 600/1552) -11.3 E40-E30 N37-N41
      3 1777 1420 (N22W65 L= 14)  2467    246        9 (1728/1340-1847/1531) -76.0 W60-W72 N17-N29
      4 1692 1416 (N23W54 L=  3)   704    148        0 (1644/1365-1747/1470) -58.0 W50-W58 N19-N27
      5  250 1210 (N18E58 L=252)  1056    142        1 ( 224/1147- 271/1286)  90.0 E60-E55 N14-N23
      6 1819 1262 (N12W62 L= 12)   401     44        0 (1801/1245-1836/1285)  26.6 W59-W64 N11-N13
      7  190 1152 (N14E64 L=245)   419     72        3 ( 171/1125- 216/1183) -23.2 E67-E60 N12-N16
      8  957 1149 (N13E03 L=306)   263    109        1 ( 913/1126-1008/1165)   9.5 E06-W00 N11-N13
      9  839 1092 (N10E11 L=298)   187    116        1 ( 810/1033- 861/1119)  68.2 E12-E09 N06-N11
     10  388 1057 (N09E43 L=267)   515    141        1 ( 349/1012- 448/1085)  23.2 E45-E37 N06-N10
     11  662  899 (S01E23 L=286)  1747    440        2 ( 594/ 764- 800/1073) -58.6 E27-E15 S09-N09
     12 1521  940 (S04W31 L=341)   688    205        0 (1431/ 890-1595/1010)  18.4 W24-W36 S06-N00
     13  348  656 (S16E51 L=258)  1129    351        1 ( 291/ 517- 411/ 813)  87.1 E54-E47 S25-S05
     14 1678  728 (S19W44 L=354)   325     46        1 (1655/ 704-1695/ 756)  45.0 W42-W45 S19-S16
     15  288  644 (S17E58 L=251)  1798    180        4 ( 247/ 586- 339/ 732) -66.0 E63-E54 S20-S11
     16  949  563 (S24E08 L=302)   498    209        1 ( 854/ 541-1049/ 593) -14.0 E13-E01 S25-S21
     17 1056  333 (S42E03 L=307)  2138    341        6 ( 953/ 248-1152/ 465)  45.0 E12-W05 S49-S32
     18  627  262 (S47E47 L=263)  1857    171        2 ( 536/ 235- 704/ 302)   7.4 E66-E35 S50-S43
  

    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