Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  906 1898 (N73W44 L= 76)   927     88        1 ( 884/1861- 931/1941)   8.1 W31-W80 N68-N76
      2 1722 1488 (N15W64 L= 97)   477     75        0 (1685/1468-1755/1504)  -5.7 W59-W69 N13-N16
      3 1775 1441 (N11W68 L=100)   546     42        0 (1758/1419-1793/1470)  18.4 W64-W70 N09-N12
      4 1097 1404 (N28W14 L= 46)   247     58        0 (1079/1373-1114/1433)  76.0 W12-W15 N25-N29
      5 1728 1408 (N11W60 L= 92)   483     43        0 (1704/1390-1754/1430) -36.9 W57-W62 N09-N12
      6 1032 1290 (N22W07 L= 39)   619    205        1 ( 963/1202-1075/1372)  36.0 W00-W11 N18-N26
      7  696 1284 (N29E15 L= 17)   760    227        2 ( 636/1184- 726/1367) -77.0 E18-E13 N22-N35
      8 1734 1201 (S00W51 L= 83)  1097    218        0 (1692/1104-1769/1294)  84.8 W49-W52 S06-N05
      9  213  940 (N16E62 L=331)   392     50        1 ( 186/ 914- 240/ 959)   7.1 E65-E57 N14-N17
     10 1611  866 (S16W33 L= 66)   588    155        2 (1559/ 827-1663/ 927)  24.0 W29-W38 S17-S13
     11 1505  625 (S29W22 L= 55)   893    187        2 (1477/ 547-1538/ 702)  79.7 W20-W24 S34-S23
     12  879  431 (S26E24 L=  8)  1082    269        3 ( 780/ 378-1004/ 502)   9.0 E32-E13 S28-S23
  

    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