Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  599 1718 (N48E44 L=345)   592     98        1 ( 576/1687- 619/1751) -45.0 E47-E40 N45-N51
      2  834 1498 (N33E16 L= 13)   694    163        1 ( 763/1444- 898/1543) -28.6 E21-E10 N29-N35
      3  431 1392 (N23E44 L=344)   777    134        1 ( 407/1332- 454/1468) -77.5 E46-E42 N19-N28
      4 1380 1319 (N22W23 L= 52)  1910    431        3 (1208/1210-1566/1421) -27.5 W11-W36 N15-N28
      5  889 1349 (N23E10 L= 19)   595    116        1 ( 854/1286- 916/1400) -69.4 E12-E07 N18-N26
      6 1047 1301 (N20W01 L= 29)   275     91        1 (1022/1259-1084/1339) -59.0 E01-W03 N17-N22
      7  914 1250 (N16E08 L= 21)   261     53        1 ( 880/1227- 945/1271) -21.8 E10-E05 N15-N17
      8  674 1111 (N07E22 L=  6)  1122    246        2 ( 622/1010- 753/1227) -56.3 E26-E16 N01-N13
      9  814 1020 (N02E13 L= 16)   382    130        1 ( 749/ 982- 871/1070)  35.5 E16-E09 N00-N05
     10  567  944 (S04E29 L=360)   332    111        1 ( 524/ 926- 625/ 964)   0.0 E32-E25 S04-S02
     11 1305  254 (S52W32 L= 61)   437     88        0 (1278/ 232-1331/ 279) -26.6 W28-W36 S54-S49
  

    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