Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  394 1613 (N29E55 L=138)  1271    145        3 ( 323/1575- 444/1658)  33.7 E61-E48 N26-N32
      2 1153 1519 (N26W04 L=198)   802    235        1 (1059/1443-1237/1569) -28.6 E01-W10 N22-N29
      3  370 1448 (N17E49 L=144)  1175    154        4 ( 323/1388- 404/1516) -76.0 E52-E45 N13-N21
      4  704 1397 (N15E23 L=171)   408     64        1 ( 666/1363- 736/1433)  63.4 E25-E20 N13-N17
      5  630 1407 (N16E28 L=165)   240     40        1 ( 610/1383- 657/1429)   0.0 E29-E26 N14-N16
      6  540 1381 (N14E34 L=160)   313     75        1 ( 507/1348- 573/1408)  45.0 E35-E31 N11-N15
      7  483 1297 (N08E36 L=157)  1416    269        3 ( 384/1193- 565/1392) -38.2 E45-E30 N02-N13
      8 1682 1245 (N14W42 L=235)   414     67        1 (1644/1220-1716/1269)   0.0 W38-W44 N12-N15
      9  341 1231 (N04E47 L=147)   260     49        1 ( 313/1214- 364/1248)  14.0 E48-E44 N02-N04
     10  180 1019 (S09E61 L=133)  1648    248        7 ( 145/ 914- 212/1124)  85.6 E65-E57 S15-S02
     11  294  873 (S18E50 L=144)  2724    494        9 ( 223/ 635- 343/1022)  78.2 E55-E44 S32-S08
     12  823  941 (S11E11 L=182)   431    112        3 ( 785/ 898- 867/ 978) -21.8 E13-E08 S13-S09
     13  398  887 (S16E40 L=153)   503    116        0 ( 363/ 848- 424/ 935)  45.0 E43-E38 S18-S13
     14  572  884 (S16E28 L=166)  1180    357        2 ( 469/ 794- 664/ 934)   8.1 E35-E21 S21-S13
     15 1464  483 (S34W39 L=232)   509    137        0 (1412/ 448-1518/ 521) -18.4 W32-W44 S35-S31
     16  976  458 (S41W02 L=195)   354     55        1 ( 948/ 428-1003/ 485) -45.0 W00-W04 S42-S38
     17 1076  424 (S42W10 L=203)   270    113        0 (1042/ 407-1129/ 448)  26.6 W07-W13 S43-S40
     18 1444  213 (S54W64 L=257)  1045     99        1 (1395/ 193-1496/ 238)  -2.1 W52-W80 S54-S52
  

    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