Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1375 1553 (N32W36 L= 32)   497    188        1 (1293/1536-1450/1579)  -4.4 W29-W42 N30-N33
      2 1078 1427 (N30W11 L=  7)   739    167        1 (1009/1374-1164/1460)  14.0 W05-W17 N27-N31
      3 1644 1345 (N14W48 L= 44)  1187    184        2 (1569/1282-1715/1383) -45.0 W42-W52 N08-N17
      4 1012 1152 (N14W01 L=358)   326    126        0 ( 937/1133-1067/1167) -23.2 E03-W04 N12-N15
      5 1109  919 (S01W03 L=360)  1605    432        6 (1052/ 765-1154/1117) -90.0 W00-W05 S10-N11
      6  498 1063 (N16E34 L=323)   469    132        0 ( 442/1035- 550/1094) -26.6 E37-E29 N13-N18
      7 1411  892 (S08W22 L= 18)   782    211        2 (1326/ 829-1501/ 944) -39.3 W16-W27 S13-S03
      8 1276  792 (S12W12 L=  8)   405    103        0 (1259/ 736-1300/ 844)  63.4 W10-W13 S15-S08
      9 1288  583 (S25W10 L=  6)   602    119        1 (1257/ 527-1321/ 642)  54.5 W07-W12 S28-S21
     10  946  350 (S37E19 L=338)  1243    312        2 ( 870/ 275- 992/ 462)  40.9 E27-E12 S42-S28
  

    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