Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1278 1471 (N27W27 L=183)   521    110        2 (1237/1434-1329/1504)  14.0 W23-W31 N26-N28
      2  521 1230 (N29E29 L=127)  1797    584        4 ( 376/1056- 742/1406) -46.1 E42-E15 N15-N42
      3  214 1128 (N26E57 L= 99)  2066    390        4 ( 178/ 957- 260/1246) -80.5 E61-E52 N15-N33
      4 1684 1111 (S03W44 L=200)   935    139        1 (1628/1085-1742/1147) -18.4 W39-W48 S05-N00
      5 1544 1063 (S02W32 L=189)   635    157        1 (1512/ 986-1591/1124)  40.6 W29-W37 S06-N00
      6 1688 1021 (S08W42 L=198)   564    197        1 (1626/ 975-1748/1044) -41.2 W37-W46 S12-S05
      7  273  961 (N16E53 L=103)   437    125        1 ( 239/ 921- 300/1011)  35.5 E57-E49 N13-N18
      8  866  905 (N03E12 L=145)  1044    218        3 ( 804/ 819- 937/ 982) -73.3 E13-E07 S02-N08
      9 1554  838 (S16W29 L=185)   501    113        0 (1539/ 784-1571/ 893)  74.1 W27-W31 S19-S12
     10 1716  840 (S20W42 L=198)   460    113        1 (1670/ 812-1749/ 882)  18.4 W37-W44 S21-S17
     11  854  580 (S16E20 L=136)   768    167        0 ( 776/ 519- 899/ 613) -58.0 E23-E17 S20-S12
     12  698  281 (S35E44 L=112)  1140    204        2 ( 651/ 236- 742/ 319)  22.6 E50-E38 S38-S31
  

    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