Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  982 1877 (N48E42 L=220)   656     75        1 ( 955/1841-1016/1897) -56.3 E45-E37 N44-N50
      2  886 1815 (N39E41 L=221)   802    104        0 ( 866/1760- 914/1863) -45.0 E45-E35 N34-N43
      3 1220 1583 (N33E06 L=256)   776    146        2 (1198/1512-1246/1655) -65.6 E07-E02 N27-N38
      4 1450 1508 (N36W12 L=274)  1282    338        0 (1353/1369-1523/1615) -31.4 W02-W20 N28-N41
      5 1613 1084 (N15W33 L=295)   781    125        2 (1567/1042-1656/1132)  -8.1 W29-W36 N12-N17
      6 1572  905 (N04W35 L=297)   947    247        1 (1506/ 781-1625/ 989)  85.6 W32-W37 S04-N09
      7  462  769 (S36E30 L=231)   255     89        0 ( 418/ 756- 510/ 786)  14.0 E34-E26 S36-S34
      8  871  636 (S32W03 L=264)   240    105        0 ( 853/ 592- 885/ 689) -74.1 W01-W03 S35-S28
      9 1381  353 (S31W51 L=312)   813    104        1 (1333/ 332-1429/ 383)   5.7 W45-W55 S33-S29
     10 1135  356 (S40W34 L=295)   327     63        0 (1109/ 336-1168/ 379)   8.1 W30-W37 S41-S39
  

    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