Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  538 1655 (N26E46 L=232)   992    119        1 ( 493/1600- 593/1700)  69.4 E47-E42 N21-N29
      2 1167 1642 (N35E02 L=276)   762    132        1 (1099/1608-1246/1679)   5.7 E06-W04 N32-N36
      3 1127 1407 (N18E00 L=278)   702    131        1 (1100/1336-1163/1465)  90.0 E01-W01 N13-N21
      4 1386 1428 (N24W16 L=294)   248     46        1 (1364/1411-1411/1453) -26.6 W14-W18 N22-N25
      5  569 1377 (N08E34 L=244)   449     84        1 ( 539/1335- 598/1416) -45.0 E36-E31 N05-N10
      6 1362 1282 (N14W16 L=294)   477    147        2 (1318/1219-1409/1341)  66.0 W13-W18 N09-N18
      7  993 1224 (N05E05 L=273)   290     76        0 ( 973/1182-1014/1272)  78.7 E05-E03 N02-N07
      8 1436 1236 (N13W21 L=299)   433    108        0 (1411/1183-1478/1264)  63.4 W20-W24 N09-N15
      9 1721 1045 (N07W44 L=322)   368     49        0 (1708/1012-1738/1076)  90.0 W43-W46 N05-N09
     10 1815  974 (N06W54 L=332)   687    110        3 (1781/ 930-1852/1018)  45.0 W51-W57 N02-N08
     11  689  967 (S15E19 L=258)  2097    562        1 ( 439/ 910- 959/ 998)  18.4 E36-E03 S20-S09
     12  298  727 (S34E49 L=229)  2174    614        4 ( 173/ 596- 425/ 886)  35.1 E64-E37 S42-S22
     13 1353  666 (S22W27 L=305)   873    192        0 (1257/ 628-1456/ 703)  32.5 W21-W33 S25-S17
  

    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