Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  873 1849 (N41E40 L=177)  1007    292        0 ( 847/1788- 909/1918) -30.5 E49-E32 N36-N47
      2  742 1434 (N09E27 L=190)  1039    409        0 ( 635/1287- 820/1609)  90.0 E29-E23 S01-N21
      3  369 1450 (N03E54 L=162)   320     72        0 ( 357/1410- 386/1487) -68.2 E55-E52 N00-N05
      4 1674 1307 (N27W34 L=250)   582    110        1 (1637/1247-1706/1361) -35.5 W30-W37 N24-N30
      5  406 1252 (S08E44 L=172)   353    114        0 ( 362/1195- 448/1295)  60.3 E46-E41 S11-S04
      6  859 1269 (N02E15 L=201)   946    250        3 ( 721/1222- 983/1291)  14.9 E23-E08 N00-N05
      7 1024 1086 (S04E02 L=215)   798    246        2 ( 987/ 966-1047/1219) -54.0 E06-W02 S10-N02
      8  885  989 (S13E07 L=209)   967    295        5 ( 785/ 948-1015/1028)  17.1 E12-W01 S16-S08
      9 1911  852 (N10W68 L=285)   625    116        0 (1883/ 825-1931/ 876)  -6.3 W63-W72 N08-N11
     10  633  648 (S39E16 L=201)  1671    413        4 ( 460/ 628- 842/ 682)  16.7 E30-W00 S43-S33
     11 1212  623 (S24W23 L=239)  1141    301        3 (1037/ 595-1387/ 665)  25.5 W12-W33 S29-S18
  

    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