Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1023 1905 (N52E44 L= 90)  1300    154        1 ( 969/1857-1057/1940)  23.2 E54-E36 N46-N57
      2 1088 1832 (N47E29 L=105)   516     61        2 (1075/1797-1106/1869) -41.2 E33-E25 N44-N51
      3  702 1769 (N29E50 L= 84)   651    127        0 ( 658/1714- 749/1818) -20.6 E58-E42 N26-N32
      4  704 1627 (N20E39 L= 95)   361     38        1 ( 674/1606- 732/1646)  45.0 E40-E37 N18-N21
      5  829 1496 (N15E25 L=109)   414    172        0 ( 761/1447- 894/1531)  58.0 E28-E22 N10-N19
      6 1440  910 (S00W27 L=160)  1177    373        0 (1344/ 782-1498/1070) -30.5 W16-W33 S05-N05
      7  544  851 (S29E25 L=109)   235     68        0 ( 508/ 830- 577/ 872)  31.0 E27-E22 S30-S27
      8  382  427 (S65E50 L= 84)   635     58        0 ( 356/ 403- 402/ 463)  20.9 E63-E41 S68-S60
      9 1498  360 (S27W59 L=193)   719    153        0 (1434/ 334-1573/ 386)  45.0 W55-W64 S30-S22
     10 1288  342 (S36W45 L=179)   794     92        2 (1232/ 311-1331/ 376)  66.8 W42-W47 S39-S32
     11  766  272 (S63W15 L=149)   343     38        0 ( 753/ 246- 782/ 299) -26.6 W11-W19 S64-S59
     12  883  207 (S62W34 L=168)   628     76        2 ( 845/ 180- 914/ 235)  80.5 W29-W40 S64-S58
  

    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