Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  976 1550 (N35W10 L=131)   949    257        1 ( 869/1519-1085/1575) -14.9 W02-W17 N32-N36
      2 1610 1494 (N16W51 L=172)   548     81        1 (1576/1467-1647/1517) -53.1 W47-W53 N13-N17
      3 1298 1328 (N14W23 L=144)   363    125        0 (1241/1313-1365/1346) -14.0 W19-W27 N12-N15
      4  301 1262 (N32E47 L= 74)  1923    276        1 ( 213/1165- 406/1340) -33.7 E59-E37 N24-N38
      5  838 1222 (N18E07 L=114)  2238    670        2 ( 565/1131-1114/1318)  -8.1 E25-W10 N13-N21
      6  146 1164 (N29E67 L= 54)   756     73        2 ( 126/1130- 166/1197) -76.0 E71-E64 N26-N30
      7 1692  841 (S23W40 L=161)   970    213        0 (1638/ 752-1719/ 939)  51.3 W34-W43 S27-S17
      8 1394  736 (S22W16 L=137)   733    270        1 (1335/ 626-1452/ 850)  45.0 W10-W21 S27-S16
      9 1866  787 (S32W60 L=181)   879     85        1 (1830/ 755-1899/ 823)  11.3 W55-W65 S33-S29
     10 1521  611 (S33W24 L=145)   287     71        0 (1496/ 581-1537/ 645)  45.0 W21-W24 S34-S31
     11  523  559 (S14E43 L= 78)   622    127        0 ( 504/ 488- 538/ 619)  60.9 E46-E40 S18-S09
     12 1042  546 (S26E10 L=111)   296     58        1 (1011/ 509-1063/ 577)  53.1 E12-E08 S27-S23
  

    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