Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  623 1713 (N48E41 L=344)   545     94        1 ( 601/1684- 643/1747) -39.8 E44-E37 N45-N51
      2  568 1635 (N41E42 L=344)   332     34        0 ( 549/1614- 586/1654)  56.3 E43-E39 N39-N42
      3  885 1489 (N32E12 L= 14)   521    126        0 ( 830/1459- 947/1521) -12.5 E15-E06 N30-N34
      4  469 1366 (N22E41 L=345)   835    212        3 ( 445/1261- 492/1455) -67.4 E42-E36 N15-N27
      5 1382 1333 (N23W23 L= 49)  2078    530        2 (1126/1228-1585/1444) -22.1 W05-W37 N16-N29
      6  941 1340 (N22E06 L= 19)   569    119        3 ( 905/1276- 969/1391) -74.1 E09-E04 N18-N25
      7 1091 1298 (N20W03 L= 29)   250     61        2 (1064/1272-1127/1328) -36.9 W01-W05 N18-N22
      8  595 1303 (N19E30 L=356)   232     58        1 ( 570/1278- 628/1323)  45.0 E31-E27 N16-N20
      9  967 1247 (N17E04 L= 21)   275     64        1 ( 929/1221-1000/1271) -26.6 E06-E02 N14-N18
     10  745 1090 (N06E18 L=  8)  1370    339        3 ( 662/ 983- 878/1225) -45.0 E23-E08 N00-N14
     11  614  940 (S04E26 L=360)   238     46        1 ( 582/ 921- 643/ 955)   0.0 E27-E23 S04-S02
  

    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