Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  810 1715 (N57W08 L= 42)   534     69        3 ( 781/1682- 835/1747)  26.6 W04-W12 N53-N59
      2 1549 1638 (N27W59 L= 93)   420     92        2 (1517/1615-1587/1679) -63.4 W55-W63 N25-N29
      3 1393 1551 (N27W40 L= 75)   412    121        1 (1354/1524-1446/1575) -15.9 W36-W43 N24-N28
      4 1217 1374 (N22W21 L= 55)   361    159        0 (1136/1351-1271/1389) -29.7 W16-W24 N19-N24
      5  989 1182 (N17W02 L= 36)   900    260        1 ( 878/1101-1060/1253) -67.4 E03-W05 N11-N23
      6  709 1033 (N15E18 L= 16)   230    112        0 ( 695/ 988- 740/1072)  53.1 E20-E15 N12-N16
      7 1051  972 (N04W00 L= 35)   520    105        1 (1012/ 927-1077/1024)  68.2 E02-W02 N01-N06
      8 1789  867 (S23W48 L= 83)  2076    384        2 (1740/ 755-1837/1013)  60.3 W42-W55 S31-S15
      9 1877  896 (S25W60 L= 94)  1229    183        0 (1828/ 819-1924/ 963)  18.4 W52-W68 S28-S22
     10  962  821 (S03E08 L= 26)   317     80        1 ( 937/ 782- 981/ 862)  78.7 E09-E07 S05-N00
     11 1169  629 (S19E01 L= 33)   452    174        0 (1093/ 601-1231/ 652) -29.7 E06-W02 S22-S17
     12 1123  461 (S29E09 L= 25)   597    161        0 (1052/ 429-1181/ 506)   5.7 E13-E03 S30-S27
     13 1400  468 (S37W10 L= 44)   605     73        1 (1364/ 436-1429/ 505)  45.0 W07-W12 S39-S35
     14  444  383 (S20E62 L=332)   550     69        2 ( 419/ 357- 485/ 410)  10.3 E66-E55 S21-S18
     15 1296  367 (S42E00 L= 34)   341     97        1 (1253/ 336-1338/ 393)   0.0 E03-W03 S43-S40
     16 1172  305 (S43E12 L= 23)   395     75        1 (1147/ 278-1200/ 337) -90.0 E13-E10 S46-S39
  

    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