Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1387 1719 (N46W20 L=259)  1611    214        4 (1295/1661-1490/1755)  -3.2 W11-W29 N42-N49
      2  285 1510 (N18E64 L=175)  1773    159        4 ( 253/1441- 321/1583) -51.3 E69-E60 N13-N23
      3 1216 1412 (N20W07 L=247)  1202    363        1 (1076/1295-1336/1508) -30.5 E02-W15 N14-N24
      4 1804 1334 (N25W55 L=294)  2215    434        2 (1716/1172-1869/1459) -58.0 W47-W64 N15-N31
      5 1681 1310 (N21W41 L=281)  1196    266        2 (1603/1276-1787/1360)   0.0 W34-W51 N18-N23
      6  158 1034 (S12E65 L=174)   750    135        2 ( 129/1004- 193/1072) -20.6 E70-E60 S13-S09
      7  324 1016 (S12E47 L=192)   429     83        0 ( 297/ 975- 357/1063)  51.3 E49-E44 S14-S09
      8  711  875 (S18E18 L=221)   497    240        0 ( 620/ 825- 823/ 909)  24.8 E23-E10 S21-S14
      9  215  844 (S23E60 L=180)   383     76        0 ( 181/ 829- 246/ 864)   7.1 E64-E55 S23-S21
     10  490  700 (S31E34 L=205)   557    118        0 ( 451/ 657- 537/ 748)  36.9 E38-E29 S33-S27
     11 1384  599 (S26W30 L=270)  2033    474        2 (1206/ 466-1556/ 747)  54.0 W21-W39 S37-S15
     12  390  528 (S43E49 L=190)   366     42        0 ( 373/ 498- 408/ 552)  33.7 E52-E46 S45-S41
  

    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