Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  422 1655 (N23E66 L= 40)   537     69        1 ( 390/1636- 455/1670)   0.0 E71-E60 N22-N24
      2 1085 1456 (N25E06 L=100)   892    331        6 (1032/1360-1132/1576) -50.2 E13-E01 N19-N31
      3  887 1391 (N16E17 L= 90)   462    138        4 ( 847/1339- 922/1432) -45.0 E19-E14 N13-N19
      4 1119 1322 (N17E01 L=105)   701    192        2 (1049/1286-1206/1349)  32.0 E05-W04 N14-N20
      5  180 1079 (S16E64 L= 42)   726     99        1 ( 158/1041- 198/1133)  90.0 E67-E62 S18-S12
      6 1754 1029 (N15W49 L=155)   387     92        1 (1739/ 983-1766/1075) -68.2 W47-W50 N11-N17
      7  496 1000 (S14E33 L= 74)   243     79        0 ( 464/ 974- 530/1016)  31.0 E34-E29 S16-S13
      8  800  693 (S27E07 L=100)  1341    390        7 ( 632/ 604- 928/ 794)  46.6 E17-W00 S36-S17
      9 1177  669 (S20W18 L=124)   224     94        1 (1139/ 637-1211/ 695)  63.4 W15-W19 S21-S17
     10 1616  436 (S23W62 L=168)  1233    102        2 (1577/ 401-1649/ 487) -71.6 W58-W68 S25-S19
     11  680  417 (S50E11 L= 96)   318     83        0 ( 651/ 401- 719/ 439)   0.0 E13-E06 S50-S48
     12  957  347 (S47W16 L=122)  1529    285        3 ( 835/ 326-1094/ 372)  26.6 W05-W25 S51-S41
     13 1469  317 (S34W59 L=165)   415     62        0 (1438/ 301-1500/ 338)  76.0 W56-W61 S36-S31
  

    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