Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1745 (N62W06 L=243)   626    112        1 ( 708/1705- 803/1784)   3.8 W00-W15 N58-N64
      2 1398 1706 (N34W53 L=290)  1877    316        4 (1263/1659-1532/1747) -37.4 W44-W61 N27-N40
      3 1118 1647 (N40W28 L=266)   277     41        0 (1100/1623-1141/1675)  76.0 W27-W30 N38-N42
      4 1161 1486 (N29W22 L=260)   225     67        1 (1133/1467-1202/1499)   0.0 W20-W25 N28-N30
      5  435 1210 (N32E35 L=202)   706    198        0 ( 363/1188- 519/1235) -15.9 E41-E27 N29-N34
      6 1847  923 (S24W55 L=292)   663     67        0 (1829/ 881-1867/ 958)  45.0 W52-W57 S25-S21
      7 1805  806 (S29W49 L=286)  1026    102        3 (1768/ 768-1859/ 844) -26.6 W44-W56 S32-S25
      8 1131  665 (S16E03 L=234)   273     97        0 (1096/ 651-1184/ 690) -21.8 E05-W00 S17-S13
      9 1086  508 (S25E11 L=227)   709    260        0 (1022/ 467-1143/ 570) -14.0 E14-E06 S27-S21
     10  802  444 (S22E30 L=207)   256     65        0 ( 762/ 426- 836/ 462) -36.9 E32-E28 S23-S19
     11  656  393 (S22E43 L=195)  1441    203        5 ( 584/ 356- 745/ 442) -48.8 E48-E36 S25-S16
     12  720  308 (S29E44 L=193)   367     54        1 ( 689/ 283- 750/ 329) -45.0 E45-E41 S31-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