Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1491 1531 (N24W45 L=106)   528     95        0 (1447/1495-1535/1561) -63.4 W43-W47 N20-N26
      2 1126 1494 (N32W19 L= 80)   404     68        1 (1093/1455-1153/1522) -78.7 W17-W20 N29-N34
      3 1497 1358 (N14W38 L= 98)   331    116        1 (1482/1297-1506/1415)  54.5 W35-W40 N10-N17
      4 1694 1319 (N07W51 L=111)  1333    222        6 (1614/1258-1777/1382) -30.3 W45-W58 N01-N11
      5 1523 1094 (S01W31 L= 92)   468    191        1 (1510/1008-1535/1193)  73.3 W29-W33 S05-N05
      6  942 1127 (N15E03 L= 58)   363    217        0 ( 903/1042-1005/1175)  45.0 E06-W02 N10-N16
      7 1914 1046 (S16W67 L=127)  3225    255        7 (1886/ 938-1944/1150)  63.4 W61-W74 S22-S09
      8 1012 1092 (N11W01 L= 61)   206    109        0 ( 997/1055-1046/1142) -90.0 W00-W02 N08-N14
      9 1855  997 (S16W57 L=118)   970    223        1 (1836/ 913-1876/1085) -90.0 W54-W59 S21-S10
     10  474 1023 (N18E35 L= 26)  1441    428        2 ( 320/ 967- 643/1070) -24.6 E47-E23 N11-N22
     11 1800  928 (S18W50 L=111)  1172    207        3 (1745/ 867-1844/ 989) -55.0 W44-W55 S23-S13
     12 1343  760 (S15W13 L= 73)   461    133        0 (1301/ 703-1378/ 818)  39.8 W09-W15 S17-S12
     13 1581  517 (S39W27 L= 88)  5464    825        4 (1220/ 312-1844/ 699)   5.5 E05-W58 S44-S34
     14 1121  644 (S17E03 L= 58)   296    100        0 (1087/ 617-1164/ 685) -63.4 E04-W00 S19-S13
     15  561  550 (S11E40 L= 20)   359     84        0 ( 525/ 529- 602/ 580) -36.9 E42-E37 S13-S09
     16 1436  224 (S64W14 L= 74)   709     59        2 (1414/ 202-1464/ 257) -54.0 W10-W20 S69-S58
     17 1391  204 (S64W07 L= 67)   656     77        0 (1369/ 178-1415/ 241) -72.9 W03-W11 S70-S57
  

    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