Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  474 1554 (N30E50 L=  2)   403     68        0 ( 458/1516- 500/1578)  90.0 E52-E47 N26-N31
      2  759 1441 (N25E22 L= 30)  1272    328        1 ( 662/1367- 810/1557)  81.3 E28-E18 N20-N33
      3 1456 1513 (N37W29 L= 81)   272     50        0 (1425/1501-1488/1526)  18.4 W25-W31 N36-N38
      4 1682 1349 (N27W47 L= 99)  3494    605        1 (1487/1147-1871/1511) -31.7 W31-W66 N15-N37
      5 1693 1192 (N17W46 L= 98)   324    111        0 (1666/1142-1708/1236) -63.4 W43-W47 N14-N20
      6 1557 1158 (N14W34 L= 86)   618    141        1 (1516/1120-1602/1215) -39.8 W31-W37 N11-N17
      7  220  948 (S11E59 L=353)   527     52        0 ( 194/ 929- 253/ 969)   0.0 E62-E55 S12-S09
      8  322  831 (S17E49 L=  3)   323     43        0 ( 304/ 806- 343/ 859) -71.6 E50-E46 S18-S15
      9  641  795 (S16E23 L= 29)   243     68        1 ( 613/ 777- 680/ 812)  36.9 E24-E20 S17-S14
     10  609  725 (S21E25 L= 27)   270     71        0 ( 595/ 686- 622/ 761) -90.0 E25-E24 S23-S18
     11 1631  650 (S17W47 L= 99)  1082    241        0 (1544/ 576-1705/ 734)  50.7 W41-W51 S22-S10
     12  700  639 (S26E18 L= 34)   301     72        0 ( 665/ 619- 742/ 658) -11.3 E20-E15 S26-S24
     13 1439  489 (S29W37 L= 89)   294     62        0 (1408/ 475-1466/ 516)  53.1 W34-W38 S30-S26
  

    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