Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1626 1099 (N10W36 L= 12)   488    190        1 (1587/1011-1651/1194) -59.0 W31-W37 N04-N14
      2  404 1133 (S10E41 L=295)   302     62        1 ( 385/1099- 418/1168) -76.0 E42-E39 S11-S07
      3 1266 1110 (N03W12 L=349)   228     64        0 (1236/1081-1289/1135) -33.7 W10-W14 N01-N03
      4  305 1063 (S15E48 L=288)  1147    311        3 ( 199/1017- 377/1118)   3.4 E59-E42 S17-S11
      5  129  916 (S25E65 L=271)  1564    332        2 (  97/ 836- 161/ 987)  71.6 E68-E61 S29-S20
      6 1789  919 (N03W53 L= 29)   727    173        1 (1767/ 869-1812/ 973) -45.0 W50-W56 N00-N06
      7 1334  722 (S18W25 L=  1)   456    171        0 (1246/ 695-1399/ 769)  36.9 W19-W28 S20-S14
      8 1699  675 (S12W52 L= 29)  1330    414        6 (1588/ 588-1792/ 728)  -2.4 W41-W65 S14-S09
      9 1697  583 (S17W56 L= 33)   516    132        2 (1670/ 540-1729/ 626) -29.7 W51-W59 S19-S15
     10 1276  602 (S26W25 L=  1)   287     76        0 (1235/ 580-1324/ 621)  31.0 W22-W27 S28-S24
     11  393  586 (S44E41 L=295)   422     98        0 ( 359/ 565- 428/ 607)  14.0 E45-E37 S46-S42
     12 1506  464 (S29W47 L= 23)  3721    734        8 (1276/ 421-1727/ 521)  16.6 W30-W67 S35-S20
  

    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