Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  727 1793 (N63W04 L= 20)  2558    379        5 ( 566/1691- 877/1880)   7.3 E18-W29 N59-N66
      2  684 1485 (N40E10 L=  6)  1118    288        0 ( 559/1443- 826/1523) -26.6 E20-W00 N33-N44
      3 1647 1413 (N09W50 L= 65)  1151    192        1 (1598/1341-1691/1497)  78.7 W47-W52 N04-N14
      4 1448 1415 (N15W36 L= 51)   292     82        0 (1411/1396-1487/1435)  -9.5 W32-W38 N14-N15
      5  872 1345 (N26E01 L= 15)   498     94        0 ( 823/1306- 922/1381)   9.5 E04-W02 N24-N27
      6  192 1339 (N42E63 L=312)   561     50        1 ( 165/1322- 226/1353)  -8.7 E70-E57 N41-N43
      7 1577 1307 (N06W40 L= 56)   266     76        1 (1542/1291-1621/1327) -45.0 W37-W43 N03-N07
      8 1832 1269 (S04W61 L= 77)   387     63        0 (1811/1248-1853/1288) -33.7 W58-W63 S05-S02
      9 1596 1184 (S01W38 L= 53)   258     53        0 (1574/1159-1619/1215)  21.8 W35-W40 S02-N00
     10 1504 1074 (S05W29 L= 44)   223     44        0 (1486/1048-1526/1091)  18.4 W27-W30 S06-S03
     11  799  886 (N02E16 L=  0)   401    111        0 ( 760/ 830- 827/ 930) -90.0 E17-E14 S01-N05
     12  164  825 (N11E67 L=309)  1453    187        1 ( 117/ 767- 223/ 908) -45.0 E73-E60 N07-N16
     13 1167  761 (S14W02 L= 17)   207     66        1 (1151/ 727-1178/ 802)  63.4 W00-W02 S15-S11
     14  813  636 (S12E22 L=354)   346     99        0 ( 793/ 581- 834/ 683)  50.2 E24-E19 S15-S09
     15 1747  565 (S44W44 L= 60)   377     54        0 (1716/ 549-1773/ 580) -12.5 W39-W48 S46-S42
     16 1450  452 (S42W13 L= 29)  1629    398        2 (1284/ 359-1621/ 515)  -1.9 E01-W29 S45-S38
     17  326  467 (S12E67 L=309)   670     60        0 ( 312/ 432- 344/ 502)  21.8 E71-E61 S14-S10
     18  472  462 (S15E52 L=324)   390     74        0 ( 430/ 449- 508/ 474)  -8.1 E55-E48 S16-S14
  

    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