Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  387 1642 (N34E67 L=306)   736     63        2 ( 374/1607- 408/1678) -22.6 E75-E62 N31-N36
      2  607 1594 (N33E39 L=335)  2262    393        6 ( 464/1456- 718/1666) -17.9 E56-E25 N25-N37
      3  494 1443 (N22E42 L=331)   477    225        0 ( 423/1379- 537/1544)  76.0 E46-E40 N17-N29
      4  425 1479 (N24E50 L=324)   369     54        0 ( 411/1447- 444/1506) -76.0 E51-E47 N21-N25
      5 1056 1396 (N25E02 L= 12)   350    109        0 (1034/1345-1091/1443)  81.9 E02-W00 N21-N28
      6  918 1338 (N20E10 L=  4)   411    168        0 ( 840/1307- 991/1361)  -5.7 E14-E04 N18-N20
      7  948 1057 (N03E05 L=  9)   217     55        1 ( 913/1033- 979/1071)  18.4 E06-E03 N01-N03
      8  793 1025 (S01E14 L=360)   245     98        0 ( 747/1011- 842/1046)   0.0 E17-E10 S01-N00
      9 1258 1012 (N03W14 L= 28)   631    193        1 (1198/ 951-1342/1046)  20.6 W11-W19 S01-N04
     10  153  933 (S13E68 L=306)  1667    164        6 ( 117/ 858- 180/ 994) -48.8 E74-E63 S17-S09
     11 1604  690 (S14W43 L= 56)  1718    304        2 (1498/ 606-1708/ 791)  45.0 W36-W49 S20-S07
     12  331  638 (S30E51 L=323)   818    117        1 ( 293/ 590- 387/ 672)  23.2 E57-E43 S33-S27
  

    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