Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1727 1500 (N13W64 L=300)  1546    193        2 (1627/1459-1810/1547) -36.3 W55-W73 N07-N18
      2  340 1447 (N49E47 L=189)   356     42        1 ( 320/1420- 361/1473)  26.6 E49-E44 N47-N50
      3 1050 1375 (N26W11 L=247)   655    149        2 ( 989/1327-1111/1418) -54.5 W08-W13 N22-N30
      4 1094 1275 (N19W10 L=247)   448    168        0 (1012/1248-1171/1295) -24.0 W05-W14 N16-N20
      5 1638 1159 (S02W40 L=276)   608    155        0 (1612/1077-1665/1220)  60.3 W37-W43 S06-N01
      6 1110 1096 (N09W06 L=243)   246    138        0 (1076/1039-1144/1155)  35.5 W03-W10 N06-N11
      7  280  915 (N16E53 L=184)   690    131        1 ( 237/ 851- 317/ 955) -45.0 E56-E50 N12-N19
      8  243  819 (N11E59 L=177)  1405    217        2 ( 192/ 738- 310/ 908) -57.5 E64-E52 N05-N17
      9 1878  869 (S27W58 L=294)   855     93        0 (1831/ 848-1920/ 900) -23.2 W51-W66 S29-S23
     10  492  860 (N10E36 L=200)   356     62        0 ( 466/ 822- 516/ 887) -63.4 E37-E34 N07-N11
     11  460  556 (S08E48 L=188)   420    109        0 ( 419/ 512- 487/ 598) -80.5 E50-E46 S10-S04
     12 1421  410 (S42W10 L=246)  3248    652        3 (1159/ 251-1657/ 582)   9.5 E16-W32 S47-S36
     13 1219  533 (S26E01 L=235)   355     76        0 (1197/ 495-1238/ 580) -90.0 E02-W00 S29-S22
     14  667  493 (S15E36 L=200)  1025    244        2 ( 600/ 416- 727/ 558) -71.6 E39-E32 S21-S11
     15  811  440 (S22E29 L=207)   443    109        0 ( 770/ 412- 843/ 484)  23.2 E32-E25 S23-S19
     16  576  424 (S18E46 L=190)   292     45        0 ( 557/ 400- 597/ 450) -90.0 E47-E44 S19-S16
     17  389  363 (S20E71 L=165)  1203     68        1 ( 366/ 327- 409/ 396)  16.7 E84-E64 S23-S17
     18  778  256 (S34E43 L=193)   582    143        0 ( 759/ 196- 796/ 329)  40.2 E50-E37 S39-S28
  

    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