Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  881 1694 (N48E15 L=122)  2851    474        5 ( 788/1563- 989/1760)  51.8 E22-E05 N37-N54
      2 1451 1663 (N47W38 L=174)   727    124        1 (1411/1636-1493/1683)  16.7 W32-W42 N44-N48
      3 1313 1607 (N42W23 L=159)   726    142        0 (1286/1551-1347/1660)  56.3 W20-W27 N37-N46
      4 1625 1580 (N40W50 L=187)  1010    159        2 (1563/1536-1696/1629) -37.9 W45-W56 N36-N43
      5  463 1424 (N27E42 L= 94)   502    192        1 ( 448/1343- 486/1496) -59.0 E44-E38 N21-N31
      6 1328 1430 (N29W21 L=157)   263     50        1 (1308/1403-1351/1465) -63.4 W19-W22 N27-N31
      7 1012 1306 (N21E01 L=135)   425     80        1 ( 973/1279-1069/1338) -31.0 E03-W02 N19-N22
      8  943 1227 (N16E06 L=131)   653    261        2 ( 807/1201-1071/1259)   3.6 E14-W02 N13-N17
      9 1604 1071 (N06W38 L=175)   393     88        1 (1573/1037-1644/1112) -45.0 W36-W41 N04-N08
     10 1303  895 (S04W18 L=154)   336    103        0 (1290/ 843-1329/ 948)  80.5 W16-W19 S07-S01
     11  350  763 (S15E47 L= 89)   566    112        4 ( 308/ 736- 382/ 787)  -8.1 E51-E44 S17-S13
     12 1030  728 (S15W01 L=138)   273     78        0 (1007/ 692-1047/ 762)  90.0 W00-W02 S17-S13
     13 1008  646 (S21E00 L=136)   304     38        1 ( 983/ 622-1032/ 675) -45.0 E01-W01 S22-S18
     14 1668  446 (S35W55 L=192)   526    145        1 (1634/ 437-1709/ 459)  -7.1 W51-W59 S36-S34
     15  501  394 (S42E44 L= 93)  1387    271        1 ( 440/ 353- 548/ 458)  29.4 E52-E36 S45-S35
     16  448  425 (S39E48 L= 89)   486    148        1 ( 420/ 399- 484/ 451)  24.0 E51-E42 S41-S36
     17 1327  373 (S41W27 L=163)  1700    341        4 (1184/ 328-1503/ 431)  13.5 W15-W40 S44-S35
  

    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