Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1814 (N37E42 L=166)   607    144        1 ( 818/1768- 865/1853) -31.0 E46-E36 N33-N39
      2  528 1076 (S16E31 L=176)   797    228        1 ( 442/ 981- 599/1151)  61.7 E34-E27 S23-S10
      3 1829 1014 (N18W52 L=260)   580    116        0 (1802/ 960-1847/1064) -39.8 W48-W55 N15-N20
      4  999  961 (S10W00 L=208)   290    128        0 ( 945/ 918-1037/ 993)  11.3 E02-W03 S12-S08
      5  654  918 (S22E19 L=188)   275     95        0 ( 633/ 887- 693/ 972) -26.6 E21-E15 S23-S19
      6 1744  883 (N08W49 L=256)  1299    241        0 (1704/ 790-1807/ 971) -21.8 W42-W58 N04-N11
      7 1265  843 (S09W18 L=226)   443    132        2 (1243/ 784-1280/ 914) -69.4 W16-W20 S12-S04
      8 1205  791 (S14W17 L=224)   209     89        0 (1155/ 763-1235/ 812)   0.0 W13-W19 S14-S12
      9  448  687 (S42E30 L=177)   942    270        0 ( 346/ 665- 563/ 708)  19.3 E40-E20 S44-S37
     10 1205  512 (S29W28 L=235)   263     70        0 (1172/ 495-1254/ 534)  53.1 W25-W29 S30-S26
     11 1572  482 (S18W54 L=261)   837    356        2 (1456/ 446-1673/ 522)  12.8 W42-W64 S20-S15
  

    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